(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024086098
(43)【公開日】2024-06-27
(54)【発明の名称】所定層の層剛性推定方法及び層剛性推定装置
(51)【国際特許分類】
G01M 99/00 20110101AFI20240620BHJP
G01H 17/00 20060101ALI20240620BHJP
【ＦＩ】
G01M99/00 Z
G01H17/00 Z
【審査請求】有
【請求項の数】2
【出願形態】ＯＬ
(21)【出願番号】P 2022201033
(22)【出願日】2022-12-16
(71)【出願人】
【識別番号】000166432
【氏名又は名称】戸田建設株式会社
(74)【代理人】
【識別番号】100090387
【弁理士】
【氏名又は名称】布施 行夫
(74)【代理人】
【識別番号】100090398
【弁理士】
【氏名又は名称】大渕 美千栄
(74)【代理人】
【識別番号】100213388
【弁理士】
【氏名又は名称】竹内 康司
(72)【発明者】
【氏名】成田 修英
(72)【発明者】
【氏名】保井 美敏
(72)【発明者】
【氏名】山本 健史
(72)【発明者】
【氏名】小阪 宏之
【テーマコード（参考）】
2G024
2G064
【Ｆターム（参考）】
2G024AD34
2G024CA13
2G024FA04
2G024FA06
2G024FA11
2G064AA05
2G064AB01
2G064AB02
2G064BA02
2G064BD02
2G064CC02
2G064CC41
2G064CC42
2G064CC43
2G064DD02
(57)【要約】
【課題】本発明は、加速度の測定点数を従来に比べて大幅に削減することができる所定層の層剛性推定方法及び層剛性推定装置を提供する。
【解決手段】本発明に係る所定層の層剛性推定方法は、多層の建築物における所定層の層剛性を多質点系モデルを用いて推定する方法であって、建築物の屋上に対応する質点の下から数えてｉ番目の質点における加速度及びｉ番目の質点の上下の質点の加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する工程と、得られた各加速度のフーリエ・スペクトルと、ｉ番目の質点の質量と、を用いてｉ番目の質点とｉ－１番目の質点の間にある所定層の層剛性（ｋｉ）を算出する工程と、所定層の層剛性（ｋｉ）の中からＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する工程と、を含む。
【選択図】図３
【解決手段】本発明に係る所定層の層剛性推定方法は、多層の建築物における所定層の層剛性を多質点系モデルを用いて推定する方法であって、建築物の屋上に対応する質点の下から数えてｉ番目の質点における加速度及びｉ番目の質点の上下の質点の加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する工程と、得られた各加速度のフーリエ・スペクトルと、ｉ番目の質点の質量と、を用いてｉ番目の質点とｉ－１番目の質点の間にある所定層の層剛性（ｋｉ）を算出する工程と、所定層の層剛性（ｋｉ）の中からＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する工程と、を含む。
【選択図】図３
【特許請求の範囲】
【請求項1】
多層の建築物における所定層の層剛性を多質点系モデルを用いて推定する方法であって、
前記建築物の屋上に対応する質点の下から数えてｉ番目の質点における加速度及び前記ｉ番目の質点の上下の質点の加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する工程と、
前記工程で得られた前記各加速度のフーリエ・スペクトルと、前記ｉ番目の質点の質量と、を用いて下記式（１）及び下記式（２）により前記ｉ番目の質点とｉ－１番目の質点の間にある前記所定層の層剛性（ｋｉ）を算出する工程と、
前記所定層の層剛性（ｋｉ）の中から下記式（３）で得られるＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する工程と、
を含むことを特徴とする、所定層の層剛性推定方法。
前記建築物の屋上に対応する質点の下から数えてｉ番目の質点における加速度及び前記ｉ番目の質点の上下の質点の加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する工程と、
前記工程で得られた前記各加速度のフーリエ・スペクトルと、前記ｉ番目の質点の質量と、を用いて下記式（１）及び下記式（２）により前記ｉ番目の質点とｉ－１番目の質点の間にある前記所定層の層剛性（ｋｉ）を算出する工程と、
前記所定層の層剛性（ｋｉ）の中から下記式（３）で得られるＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する工程と、
を含むことを特徴とする、所定層の層剛性推定方法。
【数1】
【数2】
【数3】
【請求項2】
請求項１に記載の前記所定層の層剛性推定方法を実行する層剛性推定装置であって、
演算部及び記憶部を備え、
前記演算部は、
前記ｉ番目の質点における加速度及び前記上下の質点の加速度を取得する処理と、
取得した加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する処理と、
算出された前記各加速度のフーリエ・スペクトルと、前記記憶部に保存された前記ｉ番目の質点の質量と、を用いて前記所定層の層剛性（ｋｉ）を算出する処理と、
算出された前記所定層の層剛性（ｋｉ）の中からＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する処理と、
を実行することを特徴とする、層剛性推定装置。
演算部及び記憶部を備え、
前記演算部は、
前記ｉ番目の質点における加速度及び前記上下の質点の加速度を取得する処理と、
取得した加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する処理と、
算出された前記各加速度のフーリエ・スペクトルと、前記記憶部に保存された前記ｉ番目の質点の質量と、を用いて前記所定層の層剛性（ｋｉ）を算出する処理と、
算出された前記所定層の層剛性（ｋｉ）の中からＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する処理と、
を実行することを特徴とする、層剛性推定装置。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、建築物における所定層の層剛性推定方法及び層剛性推定装置に関する。
【背景技術】
【0002】
近年、建築物の安全性を確保するために、構造ヘルスモニタリング技術が重要視されている。例えば、地震等による建築物の損傷評価のために、建築物の各層における層剛性を同定する手法が提案されている（特許文献１、非特許文献１，２）。建築物のある層の部材が地震等で損傷すると当該層の剛性が低下するため、地震の前後で層剛性を比較することで損傷の程度を評価することができる。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２０１６－１９４５１４号公報
【非特許文献】
【0004】
【非特許文献1】中村充、安井譲、「微動測定に基づく地震被災鉄骨建物の層損傷評価」、日本建築学会構造系論文集、第５１７号、ｐｐ．６１－６８、１９９９
【非特許文献2】吉村怜毅、三田彰、「多入力多出力モデルに基づく建築構造パラメタのオンライン同定」、日本建築学会構造系論文集、第５７４号ｐｐ．３９－４４、２００３
【発明の概要】
【発明が解決しようとする課題】
【0005】
一般に、高層建築物の場合には地震による構造部材の損傷は下層において発生しやすい。しかしながら、特許文献１及び非特許文献１，２の方法は、いずれも評価対象の層について層剛性を計算するために、少なくとも評価対象の層よりも上層の全ての層における加速度に基づいて層剛性を計算する必要がある。
【0006】
そこで、本発明は、加速度の測定点数を従来に比べて大幅に削減することができる所定層の層剛性推定方法及び層剛性推定装置を提供することを目的とする。
【課題を解決するための手段】
【0007】
本発明は上述の課題の少なくとも一部を解決するためになされたものであり、以下の態様または適用例として実現することができる。
【0008】
［１］本発明に係る所定層の層剛性推定方法の一態様は、
多層の建築物における所定層の層剛性を多質点系モデルを用いて推定する方法であって、
前記建築物の屋上に対応する質点の下から数えてｉ番目の質点における加速度及び前記ｉ番目の質点の上下の質点の加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する工程と、
前記工程で得られた前記各加速度のフーリエ・スペクトルと、前記ｉ番目の質点の質量と、を用いて下記式（１）及び下記式（２）により前記ｉ番目の質点とｉ－１番目の質点の間にある前記所定層の層剛性（ｋｉ）を算出する工程と、
前記所定層の層剛性（ｋｉ）の中から下記式（３）で得られるＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する工程と、
を含むことを特徴とする。
多層の建築物における所定層の層剛性を多質点系モデルを用いて推定する方法であって、
前記建築物の屋上に対応する質点の下から数えてｉ番目の質点における加速度及び前記ｉ番目の質点の上下の質点の加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する工程と、
前記工程で得られた前記各加速度のフーリエ・スペクトルと、前記ｉ番目の質点の質量と、を用いて下記式（１）及び下記式（２）により前記ｉ番目の質点とｉ－１番目の質点の間にある前記所定層の層剛性（ｋｉ）を算出する工程と、
前記所定層の層剛性（ｋｉ）の中から下記式（３）で得られるＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する工程と、
を含むことを特徴とする。
【0009】
【数1】
【0010】
【数2】
【0011】
【数3】
【0012】
［２］本発明に係る層剛性推定装置の一態様は、
上記所定層の層剛性推定方法の一態様を実行する層剛性推定装置であって、
演算部及び記憶部を備え、
前記演算部は、
前記ｉ番目の質点における加速度及び前記上下の質点の加速度を取得する処理と、
取得した加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する処理と、
算出された前記各加速度のフーリエ・スペクトルと、前記記憶部に保存された前記ｉ番目の質点の質量と、を用いて前記所定層の層剛性（ｋｉ）を算出する処理と、
算出された前記所定層の層剛性（ｋｉ）の中からＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する処理と、
を実行することを特徴とする。
上記所定層の層剛性推定方法の一態様を実行する層剛性推定装置であって、
演算部及び記憶部を備え、
前記演算部は、
前記ｉ番目の質点における加速度及び前記上下の質点の加速度を取得する処理と、
取得した加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルを算出する処理と、
算出された前記各加速度のフーリエ・スペクトルと、前記記憶部に保存された前記ｉ番目の質点の質量と、を用いて前記所定層の層剛性（ｋｉ）を算出する処理と、
算出された前記所定層の層剛性（ｋｉ）の中からＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する処理と、
を実行することを特徴とする。
【発明の効果】
【0013】
本発明に係る所定層の層剛性推定方法の一態様及び層剛性推定装置の一態様によれば、加速度の測定点数を従来に比べて大幅に削減することができる。
【図面の簡単な説明】
【0014】
【図1】本実施形態に係る層剛性推定方法に用いる建築物の模式図とその多質点系モデルを概念的に示す図である。
【図2】本実施形態に係る所定層の層剛性推定装置の構成例を示すブロック図である。
【図3】本実施形態に係る所定層の層剛性推定方法の一例を示すフローチャートである。
【図6】実施例における入力（模擬地動ｕｇ・・）である。
【図7】実施例１の層剛性推定方法の検証計算の結果を示す図である。
【図8】実施例２の層剛性推定方法の検証計算の結果を示す図である。
【発明を実施するための形態】
【0015】
以下、本発明の好適な実施形態について、図面を用いて詳細に説明する。なお、以下に説明する実施形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また、以下で説明される構成の全てが本発明の必須構成要件であるとは限らない。
【0016】
１．層剛性推定装置
図１及び図２を用いて、本発明の一実施形態に係る層剛性推定装置２を説明する。図１は、本実施形態に係る層剛性推定方法に用いる建築物１の模式図（ａ）とその多質点系モデル１ａを概念的に示す図（ｂ）であり、図２は、本実施形態に係る所定層の層剛性推定装置２の構成例を示すブロック図である。
図１及び図２を用いて、本発明の一実施形態に係る層剛性推定装置２を説明する。図１は、本実施形態に係る層剛性推定方法に用いる建築物１の模式図（ａ）とその多質点系モデル１ａを概念的に示す図（ｂ）であり、図２は、本実施形態に係る所定層の層剛性推定装置２の構成例を示すブロック図である。
【0017】
図１の（ａ）に示す多層の建築物１は、同図の（ｂ）に示す多質点系モデル１ａに置き換えることができる。各階の床に対応する質点の質量は、屋上に対応する質点から１階に向かってｍ０，ｍ１，ｍ２，…，ｍｉ－１，ｍｉ，ｍｉ＋１，…ｍｎ－１で示すことができ、初期の系の位置で地盤（ＧＬ）から鉛直方向に間隔を空けてｎ個の質点が配置される。質量ｍ１は屋上の質点の下の質点から数えて１番目の質点（１０階建てであれば１０階の床に対応する）の質量、質量ｍｉは屋上の質点の下の質点から数えてｉ番目の質点（１０階建てであれば１０－（ｉ－１）階の床に対応する）の質量を示す。上下に隣接する質点の間には、それぞれ各層の剛性ｋ１…ｎと減衰（または粘性）ｃ１…ｎを有するばねモデルが配置され、両端が上下の質点に固定される系とする。層剛性ｋ１は最上階から数えて１番目の層の剛性であり、層剛性ｋｉは最上階から数えてｉ番目の層の剛性である。本実施形態では、このｉ番目の層の層剛性（ｋｉ）を推定するものであって、言い換えれば、屋上に対応する質点（０番目）の下の質点（１番目）から数えてｉ番目の質点（以下、「ｉ番目の質点」）とｉ－１番目の質点（以下、「ｉ－１番目の質点」）の間にある所定層の層剛性（ｋｉ）を推定するものである。
【0018】
また、図１の（ｂ）において、ｕｉはｉ番目の質点の水平方向の変位の時間分布ｕｉ（ｔ）であり、ｕｇは建築物１の地盤（ＧＬ）の変位の時間分布ｕｇ（ｔ）である。さらに、変位の時間分布ｕｉ（ｔ）の時間に関する微分（速度）はドットを１つ付けてｕｉ・（ｔ）、ｕｇ・（ｔ）で表し、同様に時間に関する２階微分（加速度）はドットを２つ付けてｕｉ・・（ｔ）、ｕｇ・・（ｔ）で表す。
【0019】
建築物１は、少なくともｉ番目の質点とその上下の質点に対応する３階分の床に加速度センサ１４が取り付けられる。加速度センサ１４は、各階の床の少なくとも水平方向の加速度を計測でき、すなわち、ｉ－１番目、ｉ番目及びｉ＋１番目の質点における加速度を計測できる。建築物１は、全ての階の床に加速度センサ１４を備えてもよい。加速度センサ１４は、３軸加速度計であってもよい。
【0020】
多質点系モデル１ａは、建築物１の例えば設計緒元から各質点の質量ｍ０…ｍｎ－１があらかじめ取得でき、加速度センサ１４の出力からｉ－１番目、ｉ番目及びｉ＋１番目の質点における絶対加速度（以下単に加速度という）ｕｉ・・（ｔ）＋ｕｇ・・（ｔ）、ｕｉ－１・・（ｔ）＋ｕｇ・・（ｔ）、ｕｉ＋１・・（ｔ）＋ｕｇ・・（ｔ）を取得できる。
【0021】
このとき、多質点系モデル１ａの運動方程式（力のつり合い式）は、下記式（４）の通りである。
【0022】
【数4】
【0023】
この方程式に基づいて、下記式（１）及び下記式（２）を後述するように導き出すことができる。
【0024】
【数1】
【0025】
【数2】
【0026】
図２に示す本実施形態に係る層剛性推定装置２は、後述する所定層の層剛性推定方法の実施形態を実行する装置である。層剛性推定装置２は、上記式（１）及び上記式（２）を用いて所定層の層剛性（ｋｉ）を推定することができる。
【0027】
層剛性推定装置２は、演算部２１及び記憶部２２を備える。層剛性推定装置２は、例えばパソコンやサーバであり、図示しないＣＰＵ（中央演算処理装置）、ＲＯＭ、ＲＡＭ等のメモリやハードディスク装置等の記憶装置、外部装置との通信を行う通信インターフェース等を備える。演算部２１は、ＣＰＵやＲＡＭ等から構成することができ、記憶部２２に保存されているプログラムを実行する。記憶部２２は、ハードディスク装置等の記憶装置から構成することができ、多質点系モデル１ａの各質点の質量や上記式（１）及び上記式（２）等を記憶する。記憶部２２は、加速度センサ１４から取得した計測データを記憶してもよい。
【0028】
層剛性推定装置２は、通信ネットワークを介して建築物１に設置された複数の加速度センサ１４と接続し、計測データを受信することができる。加速度センサ１４の計測データを取得することができれば、通信ネットワーク以外の方法、例えばＵＳＢメモリによる計測データの取り込みなどでもよい。また、層剛性推定装置２は、層剛性推定装置２を操作
するためのキーボードやマウス等の公知の入力装置と、推定結果を出力するディスプレイやプリンタ等の公知の出力装置を備えてもよい。
するためのキーボードやマウス等の公知の入力装置と、推定結果を出力するディスプレイやプリンタ等の公知の出力装置を備えてもよい。
【0029】
演算部２１は、少なくともｉ番目の質点における加速度ｕｉ・・（ｔ）＋ｕｇ・・（ｔ）及び上下の質点の加速度ｕｉ－１・・（ｔ）＋ｕｇ・・（ｔ）、ｕｉ＋１・・（ｔ）＋ｕｇ・・（ｔ）を取得する処理と、取得した加速度をそれぞれフーリエ変換して各加速度のフーリエ・スペクトルｕｉ・・（ω）＋ｕｇ・・（ω）、ｕｉ－１・・（ω）＋ｕｇ・・（ω）、ｕｉ＋１・・（ω）＋ｕｇ・・（ω）を算出する処理と、算出された前記各加速度のフーリエ・スペクトルと、記憶部２２に保存されたｉ番目の質点の質量ｍｉと、を用いて所定層の層剛性（ｋｉ）を算出する処理と、を実行する。各処理については、後述の層剛性推定方法において説明する。
【0030】
演算部２１によって算出された所定層の層剛性（ｋｉ）は、建築物１を多質点系モデル１ａに置き換えたときの層剛性（ｋｉ）であるため、実際の建築物１における所定層の層剛性を推定した値となる。
【0031】
層剛性推定装置２が建築物１の地震観測結果もしくは振動測定の結果から、所定層の層剛性（ｋｉ）を推定することができれば、建築物１の損傷や劣化の度合いが把握できる。しかも、層剛性推定装置２は、従来のように多くの層の加速度データを必要とせず、少なくとも３つの質点（３つの階の床）における加速度データを取得できればよいので、超高層ビルであれば加速度の測定点数を大幅に削減することが可能となる。層剛性が推定できれば、例えば、地震の前後で比較して層剛性が低下している層があれば、その層に損傷が出ていることが分かる。また、地震前の記録が無くても、各層で層剛性を比較して不自然に層剛性の小さい層があれば、その層に損傷が生じていることが分かる。さらに、地震に限らず定期的に振動測定を行って層剛性を推定すれば、建築物１の経年劣化の進行度合いを把握できる。
【0032】
次に、上記式（４）から上記式（１）及び上記式（２）を導出できることを説明する。以下、必要に応じて上記式（４）における質量行列（ｍｉの並んでいる行列）をＭ、減衰行列（ｃｉの並んでいる行列）をＣ、剛性行列（ｋｉの並んでいる行列）をＫと書くことがある。
【0033】
まず、質点の変位、速度、加速度の時間分布をフーリエ変換すると、変位、速度、加速度の周波数分布が得られる。これをｕｉ（ω）、ｕｉ・（ω）、ｕｉ・・（ω）、ｕｇ・・（ω）のように、時間ｔを角周波数ωで置き換えた記号で書くことができる。ｕｉ（ω）、ｕｉ・（ω）、ｕｉ・・（ω）、ｕｇ・・（ω）についても上記式（４）はそのまま成り立つが、フーリエ変換の性質よりｕｉ（ω）、ｕｉ・（ω）、ｕｉ・・（ω）、ｕｇ・・（ω）の間には下記式（５）及び下記式（６）が成り立つ。
【0034】
【数5】
【0035】
【数6】
【0036】
そして、上記式（５）及び上記式（６）を上記式（４）に代入すれば、下記式（７）及び下記式（８）が成り立つ。
【0037】
【数7】
【0038】
【数8】
【0039】
また、図１の多質点系モデル１ａにおけるｊ次の固有角周波数をλｊとし、ｊ次の固有ベクトルを｛φ０，φ１ｊ，…，φｎ－１，ｊ｝とする。このとき固有角周波数と固有値ベクトルの性質より、固有ベクトル｛φ０，φ１ｊ，…，φｎ－１，ｊ｝は、角周波数がλｊで地盤の変位がゼロ（ｕｇ（ω＝λｊ）＝０）の条件下での変位ベクトルに対応するので、上記式（８）より、固有空間の運動方程式として下記式（９）が成り立つ。
【0040】
【数9】
【0041】
ここで、全質点の質量は例えば設計諸元より既知とし、所定層を含む３つの層の質点（実際は建築物１なので連続する３つの階の床）における加速度が加速度センサ１４から取得できる。現実の測定における加速度ｕｉ・・（ｔ）の記録は一定時間間隔で行われる。この時間間隔をｄｔとして、ｔｐ＝ｄｔ（ｐ－１）（ｐ＝１，２，３，…）とすると、現実の加速度記録の時間分布は、連続的なｕｉ・・（ｔ）ではなく離散的なｕｉ・・（ｔｐ）となり、そのフーリエ変換も離散的な値となる（ｕｉ・・（ωｑ）（ｑ＝１，２，３，…）と書く）。
【0042】
上記式（４）、（７）、（９）の運動方程式における既知の項を右辺に集めると下記式（１０）～（１２）となる。
【0043】
【数10】
【0044】
【数11】
【0045】
【数12】
【0046】
また、表記の簡略化のため下記式（２）及び下記式（１３）のように記号を定義することができる。
【0047】
【数2】
【0048】
【数13】
【0049】
このとき、下記式（１４）及び下記式（１５）の関係が成り立つ。
【0050】
【数14】
【0051】
【数15】
【0052】
上記式（１４）及び上記式（１５）の右辺を下記式（１６）及び下記式（１７）のように書くことができる。
【0053】
【数16】
【0054】
【数17】
【0055】
従来技術である非特許文献１では上記式（１０）が用いられ、非特許文献２では上記式（１２）が用いられ、特許文献１では上記式（１１）が用いられる。
【0056】
本発明では、フーリエ変換の積分規則（上記式（５））より、上記式（１１）を下記のように変形して用いる。
【0057】
【数18】
【0058】
上記式（１８）からｉ＋１行目だけ抜き出すと、下記式（１９）のようになる。
【0059】
【数19】
【0060】
上記式（１９）を整理すると、下記式（２０）が成り立つ。
【0061】
【数20】
【0062】
ｓｉにはその定義式（上記式（２））にωｑが入っているので周波数依存性があるし、そもそもｋｉ，ｃｉが時間変化、周波数変化しないのは、図１の多質点系モデル１ａの性質であって、実際の建築物１においては必ずしもそうでないが、ここでは角周波数がωｑ－ｂからωｑ＋ｂの範囲でｓｉおよびｓｉ＋１の変化が十分に小さいものとすると、上記
式（２０）より、下記式（２１）が成り立つ。
式（２０）より、下記式（２１）が成り立つ。
【0063】
【数21】
【0064】
さらに、上記式（２１）の左辺の変位差分の並んだ行列をＤｉ（ωｑ）と定義すれば、上記式（２１）より、ｓｉ，ｓｉ＋１が下記式（１）の通りに定まる。
【0065】
【数1】
【0066】
以上説明した通り、上記式（４）から上記式（１）及び上記式（２）が成り立つことがわかる。そして、上記式（１）または上記式（２１）より、ｓｉ，ｓｉ＋１が定まるので、非特許文献１，２及び特許文献１と同様に、ｓｉの定義（上記式（２））より、所定層の層剛性（ｋｉ）の他、ｋｉ＋１，ｃｉ，ｃｉ＋１が求まることになる。このとき、計算に用いる周波数幅ｂは最低１以上であれば、解を求めることができる。しかし、周波数幅
ｂが小さすぎると解が安定せず、大きすぎてもｓｉ，ｓｉ＋１の周波数変化を無視した平均的な結果が出てきてしまったり、後述の理由により、計算精度にも悪影響がでたりするため、例えばｂ＝１０～１００の間程度で適宜設定することが好ましい。
ｂが小さすぎると解が安定せず、大きすぎてもｓｉ，ｓｉ＋１の周波数変化を無視した平均的な結果が出てきてしまったり、後述の理由により、計算精度にも悪影響がでたりするため、例えばｂ＝１０～１００の間程度で適宜設定することが好ましい。
【0067】
次に、演算部２１は、算出された前記所定層の層剛性（ｋｉ）の中からＳ／（Ｓ＋Ｎ）比の推定値が１に近い所定基準値以上の層剛性のみを抽出する処理を実行する。当該処理の詳細は、「２．層剛性推定方法」で説明する。
【0068】
２．層剛性推定方法
図１～図５を用いて、本発明の一実施形態に係る所定層の層剛性推定方法を説明する。図３は、本実施形態に係る所定層の層剛性推定方法の一例を示すフローチャートであり、図４は、図３における区間ごとの処理の一例を示すフローチャートであり、図５は、図３におけるＳ／（Ｓ＋Ｎ）推定処理の一例を示すフローチャートである。以下の説明では、図１の（ａ）の建築物１の所定層の層剛性を、図１の（ｂ）に示す多質点系モデル１ａを用いて、図２の層剛性推定装置２で推定する。層剛性推定装置２の説明と重複する部分については省略する。
図１～図５を用いて、本発明の一実施形態に係る所定層の層剛性推定方法を説明する。図３は、本実施形態に係る所定層の層剛性推定方法の一例を示すフローチャートであり、図４は、図３における区間ごとの処理の一例を示すフローチャートであり、図５は、図３におけるＳ／（Ｓ＋Ｎ）推定処理の一例を示すフローチャートである。以下の説明では、図１の（ａ）の建築物１の所定層の層剛性を、図１の（ｂ）に示す多質点系モデル１ａを用いて、図２の層剛性推定装置２で推定する。層剛性推定装置２の説明と重複する部分については省略する。
【0069】
図３に示す所定層（ｉ番目の層）の層剛性推定方法は、多層の建築物１における所定層の層剛性（ｋｉ）を多質点系モデル１ａを用いて推定する方法であって、少なくとも加速度のフーリエ・スペクトルを算出する工程（Ｓ２１）と、所定層の層剛性（ｋｉ）を算出する工程（Ｓ３１）と、所定基準値以上の層剛性のみを抽出する工程（Ｓ４１）と、を含む。図３の例では、Ｓ１０，Ｓ１１，Ｓ１２，Ｓ２１，Ｓ２３，Ｓ２７，Ｓ２９，Ｓ３１，Ｓ３２，Ｓ４１の各処理を実行する。
【0070】
Ｓ１０：加速度データを取得する工程（Ｓ１０）は、所定層とこれに上下で隣接する２つの層の加速度センサ１４から加速度の測定データ（以下、「加速度データ」）を層剛性推定装置２が取得する。加速度データは、現在発生している地震や振動による加速度センサ１４から出力されたデータを直接取得してもよいし、過去の加速度データを記憶部２２に保存することで取得してもよい。加速度データは、建築物１の屋上に対応する質点の下から数えてｉ番目の質点における加速度ｕｉ・・（ｔ）＋ｕｇ・・（ｔ）及びｉ番目の質点の上下の質点の加速度ｕｉ－１・・（ｔ）＋ｕｇ・・（ｔ）、ｕｉ＋１・・（ｔ）＋ｕｇ・・（ｔ）である。
【0071】
次に、各工程を説明する前に、加速度データから算出された層剛性の値から信頼性の高い値を抽出するために、上記式（１）を用いて信頼性の評価指標となるＨｉ（ωｑ）を定義する。まず、上記式（１）は、変形すると下記式（２２）を得ることができる。
【0072】
【数22】
【0073】
上記式（２２）の両辺の各行に加速度フーリエ・スペクトルの共役複素数をかけると、下記式（２３）が成り立つ。
【0074】
【数23】
【0075】
上記式（２３）の左辺の行列をＨｉ（ωｑ）と定義すると下記式（２４）となる。
【0076】
【数24】
【0077】
本実施形態では、Ｈｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比を用いて測定誤差による影響が少ない層剛性の値を抽出する。ＳＮ（シグナル／ノイズ）比（本実施形態では、Ｓ／（Ｓ＋Ｎ）比を使う）は、一般的にはスカラー量に対して定義される値であり、評価対象となるスカラー量の絶対値の２乗値に含まれるシグナル成分Ｓ（真値）とノイズ成分Ｎ（測定誤差、測定と真値とのずれ）の比を表す。本実施形態であれば、Ｈｉ（ωｑ）のＳＮ比よりもＨｉ（ωｑ）を構成する各要素のＳＮ比を考えるのが望ましいとも考えられるが、それだと（層剛性の推定精度の低さに対して）過剰にＳＮ比が高くなる帯域が出てきてしまって、真値に近い推定値だけを抽出することが難しい。本実施形態においては、行列の各要素のＳＮ比を個別に考えるだけでなく、要素同士の関係性も考えることとする。そのため、本実施形態においては、行列Ｈｉ（ωｑ）における絶対値の２乗に相当する量として、ｄｅｔ（Ｈｉ
＊（ωｑ）Ｈｉ（ωｑ））（ここに、ｄｅｔＡは行列Ａの行列式を表す）を「行列Ｈｉ（ωｑ）のパワー（以下、ｐｏｗＨｉ（ωｑ）と書く）」と定義し、ｐｏｗＨｉ（ωｑ）に含まれるシグナル成分Ｓとノイズ成分Ｎの割合を推定することにより、信頼性の高い層剛性の推定値を抽出することができる。
【0078】
Ｓ１１：区間長の設定値（Ｓ１１）は、あらかじめ設定した区間長の設定値のデータをＳ１２の工程に提供する。区間長は、検討対象となる建築物１の１次固有周期に対して、例えばその１０倍～１００倍の間で設定する。
【0079】
Ｓ１２：データの区間分割（Ｓ１２）は、演算部２１がＳ１０で取得した加速度データをＳ１１で呼び出した区間長の設定値で割って分割する。例えば１００秒の加速度データがあった場合に、２０秒の区間長で５つに分割する。
【0080】
Ｓ２１：区間ごとの処理（Ｓ２１）は、加速度のフーリエ・スペクトルを算出する工程を含む。区間ごとの処理（Ｓ２１）は、演算部２１がＳ１２の工程で得られた各加速度のフーリエ・スペクトルを算出し、データのＳＮ（シグナル／ノイズ）比を求めるために上記式（２３）の左辺のＨｉ（ωｑ）を作成する。詳細な処理については、図５を用いて後述する。
【0081】
Ｓ２３：全区間処理済みか否かを確認する処理（Ｓ２３）は、Ｓ２１が全区間で処理済み（Ｙｅｓ）であればＳ２５へ進み、Ｓ２１が全区間で処理済みでない場合（Ｎｏ）にはＳ２１へ戻る。
【0082】
Ｓ２７：Ｈｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定処理（Ｓ２７）は、演算部２１がＳ２５でフィルタリングされたｍ個の区間のＨｉ（ωｑ）を用いたｐｏｗＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比である推定値Ｒｍハットを算出する。なお、詳細な処理については、図７を用いて後述する。
【0083】
Ｓ２９：Ｈｉ（ωｑ）のアンサンブル平均処理（Ｓ２９）は、演算部２１がＳ２５でフィルタリングされたＨｉ（ωｑ）のアンサンブル平均を求める。Ｈｉ（ωｑ）のアンサンブル平均を求めることにより、測定誤差（ノイズ成分）の影響を低減することができる。
【0084】
Ｓ３１：層剛性（ｋｉ）の推定処理（Ｓ３１）は、所定層の層剛性を算出する工程である。層剛性（ｋｉ）の推定処理（Ｓ３１）は、演算部２１がＳ２９の工程で得られたＨｉ（ωｑ）のアンサンブル平均と、ｉ番目の質点の質量ｍｉと、を用いて下記式（２５）及び下記式（２）により推定したｉ番目の質点とｉ－１番目の質点の間にある所定層の層剛性（ｋｉ）を算出する。下記式（２５）は、上記式（１）を変形した基本の方程式であり、上述の工程（Ｓ２９）を実行することで得られる。そして、下記式（２５）によりｓｉを求め、下記式（２）によりｓｉの実部より所定層の層剛性（ｋｉ）を算出し、虚部より減衰ｃｉを算出することができる。
【0085】
【数25】
【0086】
【数2】
【0087】
Ｓ３２：Ｓ／（Ｓ＋Ｎ）の基準値（Ｓ３２）は、あらかじめ設定した所定基準値のデータをＳ４１の工程に提供する。所定基準値は、１に近い値、例えば０．９７とすることにより、Ｓ４１の工程において層剛性の推定値における真値に近い値を抽出できる。
【0088】
Ｓ４１：Ｓ／（Ｓ＋Ｎ）が基準値以上となる周波数の層剛性の抽出処理（Ｓ４１）は、
演算部２１がＳ３１の工程で算出された所定層の層剛性（ｋｉ）の中から下記式（３）で得られるＳ／（Ｓ＋Ｎ）比の推定値がＳ３２の工程で設定した１に近い所定基準値以上の層剛性のみを抽出する。Ｓ４１の工程は、測定誤差による影響が少ない層剛性の値を抽出することがでできる。
演算部２１がＳ３１の工程で算出された所定層の層剛性（ｋｉ）の中から下記式（３）で得られるＳ／（Ｓ＋Ｎ）比の推定値がＳ３２の工程で設定した１に近い所定基準値以上の層剛性のみを抽出する。Ｓ４１の工程は、測定誤差による影響が少ない層剛性の値を抽出することがでできる。
【0089】
【数3】
【0090】
本発明に係る所定層の層剛性推定方法の一態様及び層剛性推定装置の一態様によれば、加速度の測定点数を従来に比べて大幅に削減することができる。
【0091】
２．１．区間ごとの処理
図４を用いて、上記区間ごとの処理（Ｓ２１）について説明する。図４に示す区間ごとの処理（Ｓ２１）は、Ｓ２１０，Ｓ２１２，Ｓ２１４，Ｓ２１６，Ｓ２１８，Ｓ２２０の各処理を実行する。
図４を用いて、上記区間ごとの処理（Ｓ２１）について説明する。図４に示す区間ごとの処理（Ｓ２１）は、Ｓ２１０，Ｓ２１２，Ｓ２１４，Ｓ２１６，Ｓ２１８，Ｓ２２０の各処理を実行する。
【0092】
Ｓ２１０：基線補正の処理（Ｓ２１０）は、演算部２１が図３のＳ１２で分割された各区間の加速度データに対し基線補正を実行する。
【0093】
Ｓ２１２：窓関数適用の処理（Ｓ２１２）は、演算部２１がＳ２１０で基線補正された各区間の加速度データに対し窓関数を適用する。
【0094】
Ｓ２１４：フーリエ変換の処理（Ｓ２１４）は、演算部２１がＳ２１２で窓関数が適用された各区間の加速度データ（ｕｉ・・（ｔ）＋ｕｇ・・（ｔ））に対しフーリエ変換を実行（ｕｉ・・（ω）＋ｕｇ・・（ω））する。
【0095】
Ｓ２１６：層間変位計算の処理（Ｓ２１６）は、演算部２１がＳ２１４でフーリエ変換された加速度（ｕｉ・・（ω）＋ｕｇ・・（ω））を用いて、検証対象帯域内の全ωｑに対してＨｉ（ωｑ）を作成する。
【0096】
２．２．Ｈｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定処理
図５を用いて、上記Ｈｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定処理（Ｓ２７）について説明する。図５に示すＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定処理（Ｓ２７）は、Ｓ２６６でフィルタリングされた周波数で再構築されたＨｉ（ωｑ）について演算部２１がＳ２９０，Ｓ２９２，Ｓ２９４，Ｓ２９６，Ｓ２９８，Ｓ３００の処理を実行する。
図５を用いて、上記Ｈｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定処理（Ｓ２７）について説明する。図５に示すＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定処理（Ｓ２７）は、Ｓ２６６でフィルタリングされた周波数で再構築されたＨｉ（ωｑ）について演算部２１がＳ２９０，Ｓ２９２，Ｓ２９４，Ｓ２９６，Ｓ２９８，Ｓ３００の処理を実行する。
【0097】
Ｓ２９０の工程は、演算部２１がｍ個の区間のＨｉ（ωｑ）を記憶部２２から読み込む。Ｓ２９２の工程は、演算部２１がｍ個のＨｉ
＊（ωｑ）Ｈｉ（ωｑ）のアンサンブル平均を計算する。Ｓ２９４の工程は、演算部２１がｍ個のＨｉ
＊（ωｑ）のアンサンブル平均とｍ個のＨｉ（ωｑ）のアンサンブル平均との掛け算をする。Ｓ２９６の工程は、下記式（２６）でＮを計算し、Ｓ２９８の工程は、下記式（２７）でＳを計算し、Ｓ３００の工程は、下記式（２８）でＲｍハットを計算する。ここで、Ｒｍハットは、ｍ個の区間のＨｉ（ωｑ）を用いたｐｏｗＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定値である。
【0098】
【数26】
【0099】
【数27】
【0100】
【数28】
【0101】
２．３．ｐｏｗＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定式の導出
上記式（２８）を導出した過程について以下説明する。
上記式（２８）を導出した過程について以下説明する。
【0102】
（１）記号の定義
以下、測定値（ノイズ混じりの値）と真値を区別して、測定値に「～」を付ける。ここでまず、絶対加速度は十分にＳ／（Ｓ＋Ｎ）比が高いものとして、近似的に測定値と真値が等しいものと仮定する。
以下、測定値（ノイズ混じりの値）と真値を区別して、測定値に「～」を付ける。ここでまず、絶対加速度は十分にＳ／（Ｓ＋Ｎ）比が高いものとして、近似的に測定値と真値が等しいものと仮定する。
【0103】
【数29】
【0104】
一方、相対変位については、差分をとることで地動成分（ｕｇ（ωｑ））が消去されることによってシグナル（真値）のパワーが低下するため、ノイズの影響は無視できないと考えて、シグナルとノイズの関係を次のようにおく。
【0105】
【数30】
【0106】
表記の簡略化のため、Ｈｉ（ωｑ）の各要素を以下のように書く。
【0107】
【数31】
【0108】
さらに、ｈｋ１～（ｋ＝１，２，…，ｎ），ｌ＝１，２）をシグナル成分ｈｋｌとノイズ成分εｋｌに分けて、以下のように書ける。
【0109】
【数32】
【0110】
ノイズについては、その定義より以下の性質を仮定できる。
【0111】
【数33】
【0112】
【数34】
【0113】
【数35】
【0114】
また、＜Ｈｉ（ωｑ）＞ｍの各要素を、＜ｈｋｌ～＞ｍ＋＜εｋｌ～＞ｍと書くことにすると、期待値ゼロの確率変数の性質より、以下の関係が成り立つ。
【0115】
【数36】
【0116】
さらに、ｎは十分に大きいものとし、Ｈｉ（ωｑ）に関連する任意の変数ｘｌについて、次の近似が成り立つものとする。
【0117】
【数37】
【0118】
ここで、上記式（３１）よりＨｉ
＊（ωｑ）Ｈｉ（ωｑ）は次の通りになる。
【0119】
【数38】
【0120】
上記式（３２）～上記式（３７）の性質を組わせると＜Ｈｉ
＊（ωｑ）Ｈｉ（ωｑ）＞ｍと＜Ｈｉ
＊（ωｑ）＞ｍ＜Ｈｉ（ωｑ）＞ｍは次式のようになる。
【0121】
【数39】
【0122】
【数40】
【0123】
（２）＜Ｈｉ（ωｑ）＞のＳ／（Ｓ＋Ｎ）比の推定式
まず、上記式（３９）と上記式（４０）との差分をとると次式が導ける。
まず、上記式（３９）と上記式（４０）との差分をとると次式が導ける。
【0124】
【数41】
【0125】
上記式（４１）を用いて、上記式（２６）のノイズ成分Ｎの行列について、次式の近似が成り立つ。
【0126】
【数42】
【0127】
よって、上記式（２７）のシグナル成分Ｓの行列について、次式の近似が成り立つ。
【0128】
【数43】
【0129】
上記式（４０）及び上記式（４３）より、Ｓは＜Ｈｉ
＊（ωｑ）＞ｍ＜Ｈｉ（ωｑ）＞ｍからノイズ成分を除去した形になっているので、下記式（２８）により＜Ｈｉ（ωｑ）＞ｍのＳ／（Ｓ＋Ｎ）比が推定できることが導ける。
【0130】
【数28】
【実施例0131】
（検証対象と入力データ）
本実施例は、本発明に係る層剛性推定方法を用いて算出した層剛性の推定値を検証した。検証対象とする多質点系モデルは、５質点（ｍ０＝ｍ１＝…＝ｍ４＝１ｔ）であり、各層の層剛性（ｋＮ／ｍ）は表１の通りであり、減衰はｈ１＝ｈ２＝…＝ｈ４＝０．０５とした。なお、添え字の番号の若い方が系の上部である。一次固有周波数は３．３Ｈｚとした。
本実施例は、本発明に係る層剛性推定方法を用いて算出した層剛性の推定値を検証した。検証対象とする多質点系モデルは、５質点（ｍ０＝ｍ１＝…＝ｍ４＝１ｔ）であり、各層の層剛性（ｋＮ／ｍ）は表１の通りであり、減衰はｈ１＝ｈ２＝…＝ｈ４＝０．０５とした。なお、添え字の番号の若い方が系の上部である。一次固有周波数は３．３Ｈｚとした。
【0132】
【表1】
【0133】
また、入力は、図６のようなホワイトノイズとした。サンプングレートは、１００Ｈｚ、継続時間は１６３．８４秒とした。
【0134】
（実施例１）周波数応答計算で図６の入力に対する多質点系モデルの周波数応答を求め、それを時間領域に戻して質点系応答波形の真値とした。真値をそのまま模擬地震観測記録として用いて実施例１の加速度データとした。そして、模擬地震観測記録をフーリエ変換して、周波数幅０．５Ｈｚ（中心周波数に対し±０．２５Ｈｚ）で下記式（１）によりｓｉを求め、下記式（２）によりｓｉの実部より層剛性ｋｉの推定値（ｋｉハット）を算出し、虚部より減衰ｃｉの推定値（ｃｉハット）を算出した。さらに、減衰ｃｉについて
は、粘性係数と減衰定数の関係より、ｓｉハットの虚部をｋｉハットの２倍で割った値を減衰定数の推定値（ｈｉハット）を求めた。算出された層剛性の推定値（ｋｉハット）と減衰定数の推定値（ｈｉハット）を図７に示した。
は、粘性係数と減衰定数の関係より、ｓｉハットの虚部をｋｉハットの２倍で割った値を減衰定数の推定値（ｈｉハット）を求めた。算出された層剛性の推定値（ｋｉハット）と減衰定数の推定値（ｈｉハット）を図７に示した。
【0135】
【数1】
【0136】
【数2】
【0137】
図７において、真値は破線で示し、下点基準として推定した層剛性の推定値は薄い色で幅の広い線で示し、上点基準として推定した層剛性の推定値は濃い色の幅の狭い線で示した。図７において、下点基準及び上点基準は、推定対象層の上部に位置する質点を基準点（上記式（１）及び上記式（２）の質点ｉ）として推定を行ったか、下部の質点を基準としたかを示している。ｎ層目の層剛性は、ｉ＝ｎ－１のときのｓｉ＋１と、ｉ＝ｎのときのs_iの両方から定まるので、ここではｉ＝ｎ－１で定まる層剛性を上点基準とし、ｉ＝ｎで定まる層剛性を下点基準とした。そして、図７に示すように、真値と推定結果は正確に一致しており、理論の妥当性が確認できた。
【0138】
（実施例２）次に、実施例１における質点系応答波形の真値に地動のＲＭＳ値（約０．７９ｍ／ｓ２）の１％の標準偏差を持つガウスノイズを地動および各点応答に加えて模擬地震観測記録として用いて実施例２の加速度データとした（Ｓ１０）。そして図３～図５を用いて説明したフローチャートに従って、区間ごとの処理（Ｓ２１）を全区間で実行（
Ｓ２３）した。Ｈｉ（ωｑ）のアンサンブル平均を求めたら（Ｓ２９）、周波数幅０．５Ｈｚ（中心周波数に対し±０．２５Ｈｚ）で下記式（２５）によりｓｉを求め、上記式（２）よりｓｉの実部より層剛性ｋｉ（ωｑ）の推定値（ｋｉハット）を算出した（Ｓ３１）。また、下記式（３）よりＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定値（Ｒｍハット）を算出した（Ｓ２７）。Ｓ／（Ｓ＋Ｎ）の基準値を０．９７として（Ｓ３２）、層剛性ｋｉの中からＳ／（Ｓ＋Ｎ）比の推定値（Ｒｍハット）が０．９７以上の層剛性のみを抽出した（Ｓ４１）。算出されたＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定値を図８の右側に示し、その内、推定値が０．９７以上の点を左側の層剛性の推定値（ｋｉハット）として図８の左側に示した。
Ｓ２３）した。Ｈｉ（ωｑ）のアンサンブル平均を求めたら（Ｓ２９）、周波数幅０．５Ｈｚ（中心周波数に対し±０．２５Ｈｚ）で下記式（２５）によりｓｉを求め、上記式（２）よりｓｉの実部より層剛性ｋｉ（ωｑ）の推定値（ｋｉハット）を算出した（Ｓ３１）。また、下記式（３）よりＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定値（Ｒｍハット）を算出した（Ｓ２７）。Ｓ／（Ｓ＋Ｎ）の基準値を０．９７として（Ｓ３２）、層剛性ｋｉの中からＳ／（Ｓ＋Ｎ）比の推定値（Ｒｍハット）が０．９７以上の層剛性のみを抽出した（Ｓ４１）。算出されたＨｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定値を図８の右側に示し、その内、推定値が０．９７以上の点を左側の層剛性の推定値（ｋｉハット）として図８の左側に示した。
【0139】
【数25】
【0140】
【数3】
【0141】
図８において、真値は破線で示し、下点基準として推定した層剛性の推定値は薄い色で幅の広い線（または点）で示し、上点基準として推定した層剛性の推定値は濃い色の幅の狭い線（または点）で示した。層剛性の推定値（ｋｉハット）は、真値（破線）の近傍に分布しており、Ｈｉ（ωｑ）のＳ／（Ｓ＋Ｎ）比の推定値（Ｒｍハット）を用いて、信頼性が高い層剛性の推定値（ｋｉハット）だけを抽出できた。
【0142】
本発明は、上述した実施形態に限定されるものではなく、さらに種々の変形が可能である。例えば、本発明は、実施形態で説明した構成と実質的に同一の構成（例えば、機能、方法、及び結果が同一の構成、あるいは目的及び効果が同一の構成）を含む。また、本発明は、実施形態で説明した構成の本質的でない部分を置き換えた構成を含む。また、本発
明は、実施形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施形態で説明した構成に公知技術を付加した構成を含む。
明は、実施形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施形態で説明した構成に公知技術を付加した構成を含む。
【符号の説明】
【0143】
１…建築物、１ａ…多質点系モデル、１１…層、１４…加速度センサ、２…層剛性推定装置、２１…演算部、２２…記憶部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
