特許 6032493 建物の健全性確認方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6032493

(24)【登録日】2016年11月4日

(45)【発行日】2016年11月30日

(54)【発明の名称】建物の健全性確認方法

(51)【国際特許分類】

   G01M 99/00 20110101AFI20161121BHJP

【ＦＩ】

   G01M99/00 Z

【請求項の数】1

【全頁数】16

(21)【出願番号】特願2013-127021(P2013-127021)

(22)【出願日】2013年6月17日

(65)【公開番号】特開2015-1483(P2015-1483A)

(43)【公開日】2015年1月5日

【審査請求日】2015年11月18日

(73)【特許権者】

【識別番号】000002299

【氏名又は名称】清水建設株式会社

(74)【代理人】

【識別番号】100064908

【弁理士】

【氏名又は名称】志賀 正武

(74)【代理人】

【識別番号】100108578

【弁理士】

【氏名又は名称】高橋 詔男

(74)【代理人】

【識別番号】100146835

【弁理士】

【氏名又は名称】佐伯 義文

(74)【代理人】

【識別番号】100161506

【弁理士】

【氏名又は名称】川渕 健一

(72)【発明者】

【氏名】斎藤 知生

【審査官】 萩田 裕介

(56)【参考文献】

【文献】 特開２０１２−０１４６０５（ＪＰ，Ａ）

【文献】 特開２００４−２６４２３５（ＪＰ，Ａ）

【文献】 特開２０１２−１７３００１（ＪＰ，Ａ）

【文献】 特開２００１−２００６４５（ＪＰ，Ａ）

【文献】 M. W. Vanik et al，Bayesian Probabilistic Approach to Structural Health Monitoring，Journal of Engineering Mechanics，２０００年 ７月，Vol.126,No.7，738-745

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｍ ９９／００

Ｇ０１Ｍ ７／００ − ７／０８

(57)【特許請求の範囲】

【請求項1】

建物の観測層に設置したセンサによる観測データから建物特性を表すためのパラメータとなる最尤推定値を求め、該最尤推定値を用いて建物の健全性を確認する方法であって、
建物の設計モデルの質量行列Ｍ、減衰係数行列Ｃ、剛性行列Ｋから、下記の式（１）に示す一般固有値問題を解いてｊ次の固有角振動数ｗ_ｊと刺激関数φ_ｊを得る第１工程と、
下記の式（２）に示す確率変数のモデルパラメータ（ｎ_ｐはパラメータ数）を用い、剛性分布ｋを修正する関数△ｋ（θ）を導入することにより、剛性分布をｋ’＝ｋ＋△ｋ（θ）に修正するとともに、剛性行列ＫをＫ’（θ）に修正する第２工程と、
観測層の建物応答絶対角度ｙ_ｐ（θ）を下記の式（３）で表すとともに、建物応答絶対加速度の確率モデルを下記の式（４）で表す第３工程と、
地震時に、下記の式（５）で表す観測データＤをセンサから得るとともに、ベイズの定理によってθの事後分布を下記の式（６）で求める第４工程と、
下記の式（７）から事前分布のｐ（θ）、下記の式（８）から尤度関数のｐ（D|θ）を求める第５工程と、
尤度関数を推定値θ_Ｌ（上に＾）近傍で下記の式（９）のような正規分布の相似形で近似するとともに、式（９）の両辺の対数をとって下記の式（１０）を求める第６工程と、
各θ_ｊそれぞれに−γ，０，γを代入して対数尤度ｌｏｇｐ（Ｄ｜θ）を計算し、それを最大化するθ＝θ_Ｍを選ぶ第７工程と、
各θ_ｊについて、それ以外のパラメータをθ_Ｍに等しく固定した上で、θ_ｊ＝−γ，０，γに対応する対数尤度ｌｏｇＬ_ｊ（θ_ｊ）を下記の式（１１）の２次式で表す第８工程と、
式（１１）の２次式の係数α＝［α_０α_１α_２］^Ｔを下記の式（１２）、式（１３）、式（１４）によって求める第９工程と、
式（１１）を最大化するθ_ｊ＝θ_ｊ（上に＾）を、下記の式（１５）と式（１６）によって得る第１０工程と、
各θ_ｊについて、第８工程から第１０工程を行い、下記の式（１７）によって、最尤推定値θ_Ｌ（上に＾）を得る第１１工程とを備えていることを特徴とする建物の健全性確認方法。

【数1】

【数2】

【数3】

【数4】

ｎ_ｓは、建物に設置されたセンサの数（地動計測用のものを除く）、ｙ_ｐ（上に^（ハット））（θ）は、Ｍ、Ｃ、Ｋ’（θ）で規定される修正設計モデルに観測された地動ｕを入力したときの各時刻におけるセンサ設置階の応答絶対加速度であり、その値を期待値として等しい分散σ_ｙ^２で独立に正規分布していることを示す。

【数5】

【数6】

【数7】

【数8】

【数9】

ｃ、ｃ’はスケーリング係数、Ｃ_Ｌは推定値近傍でのθの誤差共分散行列である。

【数10】

【数11】

【数12】

【数13】

【数14】

【数15】

【数16】

【数17】

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、建物の健全性を確認するための方法に関する。

【背景技術】

【0002】

建物にセンサを設置し、このセンサからの情報に基づいて建物の損傷、劣化の度合いを把握し、建物の損傷検知や健全性評価を行う構造ヘルスモニタリングが注目されている（例えば、特許文献１、特許文献２参照）。また、特に、オフィスビルやマンション等の多層構造の建物に対し、地震発生時の被災状況を早期に且つ精度よく確認（把握）、判定することが求められている。

【0003】

さらに、特に大地震が発生した後に、建物特性の変化から損傷、劣化の度合い、すなわち、建物の健全性を迅速に評価、判定することが求められており、このため、建物特性の評価に要する計算時間をできるだけ短くすることが必要とされている。

【先行技術文献】

【特許文献】

【0004】

【特許文献1】特開２０１１−１３２６８０号公報

【特許文献2】特開２００１−９９７６０号公報

【発明の概要】

【発明が解決しようとする課題】

【0005】

しかしながら、上記従来の建物の健全性確認方法（構造ヘルスモニタリング）では、建物特性を表現するパラメータがセンサによる観測データに基づく最尤推定値として得られ、一般にその推定に非線形最小二乗法等の繰り返し計算を行う最適化手法を用いている。さらに、各々の繰り返しステップにおいて建物の時刻歴応答解析を行うことになる。このため、上記従来の建物の健全性確認方法では、推定に要する計算時間がどうしても長くなってしまう。

【0006】

本発明は、上記事情に鑑み、センサによる観測データから建物特性を表すためのパラメータの最尤推定値を精度よく短時間で得られるようにすることを可能にした建物の健全性確認方法を提供することを目的とする。

【課題を解決するための手段】

【0007】

上記の目的を達するために、この発明は以下の手段を提供している。

【0008】

本発明の建物の健全性確認方法は、建物の観測層に設置したセンサによる観測データから建物特性を表すためのパラメータとなる最尤推定値を求め、該最尤推定値を用いて建物の健全性を確認する方法であって、建物の設計モデルの質量行列Ｍ、減衰係数行列Ｃ、剛性行列Ｋから、下記の式（１）に示す一般固有値問題を解いてｊ次の固有角振動数ｗ_ｊと刺激関数φ_ｊを得る第１工程と、下記の式（２）に示す確率変数のモデルパラメータ（ｎ_ｐはパラメータ数）を用い、剛性分布ｋを修正する関数△ｋ（θ）を導入することにより、剛性分布をｋ’＝ｋ＋△ｋ（θ）に修正するとともに、剛性行列ＫをＫ’（θ）に修正する第２工程と、観測層の建物応答絶対角度ｙ_ｐ（θ）を下記の式（３）で表すとともに、建物応答絶対加速度の確率モデルを下記の式（４）で表す第３工程と、地震時に、下記の式（５）で表す観測データＤをセンサから得るとともに、ベイズの定理によってθの事後分布を下記の式（６）で求める第４工程と、下記の式（７）から事前分布のｐ（θ）、下記の式（８）から尤度関数のｐ（Ｄ|θ）を求める第５工程と、尤度関数を推定値θ_Ｌ（上に＾）近傍で下記の式（９）のような正規分布の相似形で近似するとともに、式（９）の両辺の対数をとって下記の式（１０）を求める第６工程と、各θ_ｊそれぞれに−γ，０，γを代入して対数尤度ｌｏｇｐ（Ｄ｜θ）を計算し、それを最大化するθ＝θ_Ｍを選ぶ第７工程と、各θ_ｊについて、それ以外のパラメータをθ_Ｍに等しく固定した上で、θ_ｊ＝−γ，０，γに対応する対数尤度ｌｏｇＬ_ｊ（θ_ｊ）を下記の式（１１）の２次式で表す第８工程と、式（１１）の２次式の係数α＝［α_０α_１α_２］^Ｔを下記の式（１２）、式（１３）、式（１４）によって求める第９工程と、式（１１）を最大化するθ_ｊ＝θ_ｊ（上に＾）を、下記の式（１５）と式（１６）によって得る第１０工程と、各θ_ｊについて、第８工程から第１０工程を行い、下記の式（１７）によって、最尤推定値θ_Ｌ（上に＾）を得る第１１工程とを備えていることを特徴とする。

【0009】

【数1】

【0010】

【数2】

【0011】

【数3】

【0012】

【数4】

ｎ_ｓは、建物に設置されたセンサの数（地動計測用のものを除く）、ｙ_ｐ（上に^（ハット））（θ）は、Ｍ、Ｃ、Ｋ’（θ）で規定される修正設計モデルに観測された地動ｕを入力したときの各時刻におけるセンサ設置階の応答絶対加速度であり、その値を期待値として等しい分散σ_ｙ^２で独立に正規分布していることを示す。

【0013】

【数5】

【0014】

【数6】

【0015】

【数7】

【0016】

【数8】

【0017】

【数9】

ｃ、ｃ’はスケーリング係数、Ｃ_Ｌは推定値近傍でのθの誤差共分散行列である。

【0018】

【数10】

【0019】

【数11】

【0020】

【数12】

【0021】

【数13】

【0022】

【数14】

【0023】

【数15】

【0024】

【数16】

【0025】

【数17】

【発明の効果】

【0026】

本発明の建物の健全性確認方法においては、従来の最適化手法による場合と同程度の精度を保ちながら、極めて短い計算時間で最尤推定値の推定が可能になる。よって、建物特性の変化を地震直後に迅速に推定でき、特に大地震後に、建物特性の変化から損傷、劣化の度合い、すなわち、建物の健全性を迅速に評価、判定することができ、有効に利用することが可能になる。

【図面の簡単な説明】

【0027】

【図1】本発明の一実施形態に係る建物の健全性確認方法において、対象の建物を示す図である。

【図2】本発明の一実施形態に係る建物の健全性確認方法において、示す図である。シミュレーション結果を示す図である。

【図3】本発明の一実施形態に係る建物の健全性確認方法と従来手法によって求めた最尤推定値を比較した図である。

【図4】本発明の一実施形態に係る建物の健全性確認方法と従来手法によって最尤推定値を計算する時間を比較した図である。

【発明を実施するための形態】

【0028】

以下、図１から図４を参照し、本発明の一実施形態に係る建物の健全性確認方法について説明する。

【0029】

ここで、本実施形態の建物の健全性確認方法は、構造ヘルスモニタリングシステムを用いてオフィスビルやマンション等の多層構造の建物の健全性を確認、把握するための方法に関するものである。

【0030】

そして、本実施形態に係る建物１は、図１に示すように、複数の振動センサ２、３、４・・・ｎが異なる階（観測層）に設けられ、これらセンサ２〜ｎによって地震に伴う振動が検知される。また、各センサの２〜ｎの検知データ（観測データ）がインターフェイス部５を介してヘルスモニタリングシステムに送信され、ヘルスモニタリングシステム側で、この検知データが地震時の建物の応答としてログデータの形で記憶される。

【0031】

また、センサ２〜ｎ、インターフェイス部５とともに構造ヘルスモニタリング装置を構成するヘルスモニタリングシステムは、例えば、図２に示すように、システムバス６、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）７、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）８、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）９、外部情報機器との送受信を行うための通信制御部１０、キーボードコントローラなどの入力制御部１１、ディスプレイコントローラなどの出力制御部１２、外部記憶装置制御部１３、キーボード、ポインティングデバイス、マウスなどの入力機器からなる入力部１４、ＬＣＤディスプレイなどの表示装置や印刷装置からなる出力部１５、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）等の外部記憶装置１６を備えて構成されている。

【0032】

ＣＰＵ７は、ＲＯＭ８内のプログラム用ＲＯＭ、あるいは大容量の外部記憶装置１６に記憶されたプログラム等に応じ、外部機器との通信を行ってデータの検索や取得、また、図形、イメージ、文字、表等が混在した出力データの処理を実行したり、外部記憶装置１６に予め格納されたデータベースの管理を実行するなど、演算処理を行う。また、ＣＰＵ７は、システムバス１０に接続される各デバイスを統括的に制御する。

【0033】

さらに、本実施形態のヘルスモニタリングシステムでは、建物の健全性を評価するにあたり、建物の観測層に設置したセンサ２〜ｎによる観測データから建物特性を表すためのパラメータとなる最尤推定値を求めるための推定アルゴリズムがＣＰＵ７に搭載されている。

【0034】

ちなみに、ＲＯＭ９内のプログラム用ＲＯＭ、あるいは外部記憶装置１６には、ＣＰＵ７の制御用の基本プログラムであるオペレーティングシステムプログラム（以下ＯＳ）等が記憶されている。また、ＲＯＭ９、あるいは外部記憶装置１６には出力データ処理等を行う際に使用される各種データが記憶されている。ＲＡＭ８は、ＣＰＵ７の主メモリ、ワークエリア等として機能する。

【0035】

また、入力制御部１１は、キーボードや不図示のポインティングデバイスからの入力部１４を制御する。出力制御部１２は、ＬＣＤディスプレイ等の表示装置やプリンタなどの印刷装置の出力部１５の出力制御を行う。

【0036】

外部記憶装置制御部１３は、例えば、ブートプログラム、各種のアプリケーション、フォントデータ、ユーザーファイル、編集ファイル、プリンタドライバ等を記憶するＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）等の外部記憶装置１６へのアクセスを制御する。

【0037】

また、通信制御部１０は、ネットワークを介して外部機器との通信を制御するものであり、これにより、システムが必要とするデータを適宜インターネットやイントラネット上の外部機器が保有するデータベースから取得したり、外部機器に情報を送信したりすることができる。

【0038】

そして、本実施形態の建物の健全性確認方法において、まず、建物１の設計モデルの質量行列Ｍ、減衰係数行列Ｃ、剛性行列Ｋが与えられており、式（１８）に示す一般固有値問題を解いてｊ次の固有角振動数ｗ_ｊと刺激関数φ_ｊを得る（第１工程）。

【0039】

【数18】

【0040】

ここで、剛性分布ｋを修正する関数△ｋ（θ）を導入する。これにより、剛性分布がｋ’＝ｋ＋△ｋ（θ）に修正され、これに対応して剛性行列ＫがＫ’（θ）に修正される（第２工程）。このとき、モデルパラメータは、式（１９）に示すように確率変数である（ｎ_ｐはパラメータ数）。

【0041】

【数19】

【0042】

一方、センサ設置階（観測層）の建物応答絶対角度ｙ_ｐ（θ）は式（２０）で表され、この建物応答絶対加速度の確率モデルは式（２１）で表せる（第３工程）。

【0043】

【数20】

【0044】

【数21】

【0045】

ｎ_ｓは、建物に設置されたセンサの数（地動計測用のものを除く）、ｙ_ｐ（上に^（ハット））（θ）は、Ｍ、Ｃ、Ｋ’（θ）で規定される修正設計モデルに観測された地動ｕを入力したときの各時刻におけるセンサ設置階の応答絶対加速度であり、その値を期待値として等しい分散σ_ｙ^２で独立に正規分布していることを示している。

【0046】

そして、地震時に、式（２２）で表す観測データＤが得られると、ベイズの定理によってθの事後分布が式（２３）で求められる（第４工程）。

【0047】

【数22】

【0048】

【数23】

【0049】

ここで、従来、ｐ（θ）は、事前分布で、式（２４）のような互いに独立で平均が０の正規分布である。また、ｐ（Ｄ|θ）は、尤度関数で、式（２５）で求められる（第５工程）。

【0050】

【数24】

【0051】

【数25】

【0052】

一方、本実施形態の建物の健全性確認方法においては、ＣＰＵ７の推定アルゴリズムによって、尤度関数を推定値θ_Ｌ（上に＾）近傍で次の式（２６）のような正規分布の相似形で近似する。ここで、ｃ、ｃ’はスケーリング係数、Ｃ_Ｌは推定値近傍でのθの誤差共分散行列である。なお、この式（２６）のような正規分布の相似形は確率分布ではない。

【0053】

【数26】

【0054】

そして、式（２６）の両辺の対数をとって式（２７）が得られ（第６工程）、本実施形態の建物の健全性確認方法では、この式（２７）が各θ_ｊの２次式になっていることに注意し、下記の手順１）〜手順６）（第７工程〜第１１工程）の手順にて最尤推定値θ_Ｌ（上に＾）を繰り返し計算なしで推定する。

【0055】

【数27】

【0056】

手順１）：各θ_ｊそれぞれに−γ，０，γを代入して対数尤度ｌｏｇｐ（Ｄ｜θ）を計算し、それを最大化するθ＝θ_Ｍを選ぶ（第７工程）。

【0057】

手順２）：各θ_ｊについて、それ以外のパラメータをθ_Ｍに等しく固定した上で、θ_ｊ＝−γ，０，γに対応する対数尤度ｌｏｇＬ_ｊ（θ_ｊ）を次の式（２８）の２次式で表す（第８工程）。

【0058】

【数28】

【0059】

手順３）：上記の式（２８）の２次式の係数α＝［α_０α_１α_２］^Ｔを次の式（２９）、式（３０）、式（３１）によって求める（第９工程）。

【0060】

【数29】

【0061】

【数30】

【0062】

【数31】

【0063】

手順４）：上記の式（２８）を最大化するθ_ｊ＝θ_ｊ（上に＾）を、次の式（３２）を考慮し式（３３）によって得る（第１０工程）。

【0064】

【数32】

【0065】

【数33】

【0066】

手順５）：各θ_ｊについて、上記の手順２）〜手順４）（第８工程〜第１０工程）を行う。

【0067】

手順６）：そして、次の式（３４）によって、最尤推定値θ_Ｌ（上に＾）を得る（第１１工程）。

【0068】

【数34】

【0069】

このように、本実施形態の建物の健全性確認方法においては、上記の手順１）から手順４）によって、繰り返し計算を用いることなく、非常に少ない計算量でθの最尤推定値θ_Ｌ（上に＾）を推定することができる。

【0070】

ここで、図３、図４は、実際に本実施形態の手法を用い、θの最尤推定値を推定し、通常の最適化手法と比べ、どの程度の精度を保持しているか、また、計算に要する時間がどの程度短縮されるかを確認した結果を示している。

【0071】

なお、このシミュレーションでは、モデルパラメータ数をｎ_ｐ＝２、建物階数を２４階、センサ数を３、時刻歴波形長さを４０００ステップとしている。また、計算機の仕様は、ＯＳ：Ｗｉｎｄｏｗｓ（登録商標） ＸＰ Ｐｒｏｆｅｓｓｉｏｎａｌ ＳＰ３、ＣＰＵ：Ｉｎｔｅｌ Ｐｒｎｔｉｕｍ Ｄ ３．２０ＧＨｚである。

【0072】

実際に本実施形態の手法を用いてθの最尤推定値を推定し、通常の最適化手法と比べた図３では、対数尤度の相対値も等高線表示している。そして、この結果から、本実施形態の手法による推定値は、従来の最適化手法による推定値をよく近似していることが分かる。

【0073】

次に、計算所要時間の比較結果を示した図４から、従来の手法では推定値を求めるために２０秒以上を要するのに対し、本実施形態の手法を用いる約０．４秒しかかからないことが分かる。

【0074】

したがって、本実施形態の建物の健全性確認方法においては、従来の最適化手法による場合と同程度の精度を保ちながら、極めて短い計算時間で最尤推定値の推定が可能になる。よって、建物特性の変化を地震直後に迅速に推定でき、特に大地震後に、建物特性の変化から損傷、劣化の度合い、すなわち、建物の健全性を迅速に評価、判定することができ、有効に利用することが可能になる。

【0075】

以上、本発明に係る建物の健全性確認方法の一実施形態について説明したが、本発明は上記の実施形態に限定されるものではなく、その趣旨を逸脱しない範囲で適宜変更可能である。

【符号の説明】

【0076】

１ 建物
２ センサ（振動センサ）
３ センサ（振動センサ）
４ センサ（振動センサ）
ｎ センサ（振動センサ）
５ インターフェイス部
６ システムバス
７ ＣＰＵ
８ ＲＡＭ
９ ＲＯＭ
１０ 通信制御部
１１ 入力制御部
１２ 出力制御部
１３ 外部記憶装置制御部
１４ 入力部
１５ 出力部
１６ 外部記憶装置

【図1】

【図2】

【図3】

【図4】