(57)【要約】【課題】  本発明は、地震動を受けた後の構造物の累積損傷を考慮した耐震性を判定することができる構造物の耐震性判定方法及び構造物の耐震性判定システムを提供する。【解決手段】  本発明に係る構造物の固有周期推定方法は、地震動を受ける構造物の加速度データから周波数ごとの振幅データを算出し（Ｓ３１０）、振幅データから構造物の一番下の階に対する一番上の階の修正スペクトル比Ｇ（ω）を算出するに際し、分母にパワースペクトルの平均値に基づく値を加算して修正スペクトル比Ｇ（ω）を算出し（Ｓ３２０）、修正スペクトル比Ｇ（ω）から増幅率が一番大きな固有周波数を選択し（Ｓ３３０）、固有周波数から一次固有周期を算出する（Ｓ３７０）。【選択図】図７
2019158731.tif 000002
【特許請求の範囲】
【請求項１】
地震動を受ける構造物の加速度データから周波数ごとの振幅データを算出し、
前記振幅データから前記構造物の一番下の階に対する一番上の階のスペクトル比を算出するに際し、分母にパワースペクトルの平均値に基づく値を加算して前記スペクトル比を算出し、
前記スペクトル比から増幅率が一番大きな固有周波数を選択し、
前記固有周波数から一次固有周期を算出することを特徴とする、構造物の固有周期推定方法。
【請求項２】
請求項１において、
前記固有周波数における前記スペクトル比が予め設定されたしきい値以上である場合に、前記固有周波数から前記一次固有周期を算出し、
前記しきい値は、２〜１０に設定することを特徴とする、構造物の固有周期推定方法。
【請求項３】
請求項１または２において、
前記固有周波数が前記スペクトル比から得られたピーク周波数のうち最も低周波数側である場合に、前記固有周波数から前記一次固有周期を算出することを特徴とする、構造物の固有周期推定方法。
【請求項４】
請求項１〜３のいずれか１項において、
前記スペクトル比は、下記式（１）で算出されることを特徴とする、構造物の固有周期推定方法。
【数１】
2019158731.tif 000014
【請求項５】
請求項１〜４のいずれか１項の構造物の固有周期推定方法で得られた前記一次固有周期に対応する伝達関数を作成し、
前記伝達関数に予め設定した設定地震動条件を入力して、地震応答解析をすることで前記設定地震動条件に対する前記構造物の層間変形角を演算し、
前記層間変形角が予め設定した許容値を超えているか否かで前記構造物の耐震性を判定することを特徴とする、構造物の耐震性判定方法。
【請求項６】
請求項５において、
前記伝達関数は、
前記振幅データにおける前記一次固有周期に対応する前記振幅データの前記一番上の階に対する加速度センサの設置階のスペクトル比を用いて前記一番上の階からの前記設置階への地震動の伝播時間を求め、さらにそれを線形補間して求まる、前記一番上の階から前記構造物の各階への地震動の伝播時間である各階伝播時間により定まることを特徴とする、構造物の耐震性判定方法。
【請求項７】
構造物の高さ方向における複数の位置の加速度を取得する加速度センサから取得された加速度に基づいて演算を行う演算部を含み、
前記演算部は、
地震動を受ける構造物の加速度データから周波数ごとの振幅データを算出し、
前記振幅データから前記構造物の一番下の階に対する一番上の階のスペクトル比を算出するに際し、分母にパワースペクトルの平均値に基づく値を加算して前記スペクトル比を算出し、
前記スペクトル比から増幅率が一番大きな固有周波数を選択し、
前記固有周波数から一次固有周期を算出することを特徴とする、構造物の固有周期推定システム。
【請求項８】
請求項７において、
前記固有周波数における前記スペクトル比が予め設定されたしきい値以上である場合に、前記固有周波数から前記一次固有周期を算出し、
前記しきい値は、２〜１０に設定することを特徴とする、構造物の固有周期推定システム。
【請求項９】
請求項７または８において、
前記固有周波数が前記スペクトル比から得られたピーク周波数のうち最も低周波数側である場合に、前記固有周波数から一次固有周期を算出することを特徴とする、構造物の固有周期推定システム。
【請求項１０】
請求項７〜９のいずれか１項において、
前記スペクトル比は、下記式（１）で算出されることを特徴とする、構造物の固有周期推定システム。
【数１】
2019158731.tif 000015
【請求項１１】
前記演算部は、
請求項７〜１０のいずれか１項の構造物の固有周期推定システムで得られた前記一次固有周期に対応する伝達関数を算出し、
前記伝達関数に予め設定した設定地震動条件を入力して、地震応答解析をすることで前記設定地震動条件に対する前記構造物の層間変形角を演算し、
前記層間変形角が予め設定した許容値を超えているか否かで前記構造物の耐震性を判定することを特徴とする、構造物の耐震性判定システム。
【請求項１２】
請求項１１において、
前記伝達関数は、
前記振幅データにおける前記一次固有周期に対応する前記振幅データの前記一番上の階
に対する加速度センサの設置階のスペクトル比を用いて前記一番上の階からの前記設置階への地震動の伝播時間を求め、さらにそれを線形補間して求まる、前記一番上の階から前記構造物の各階への地震動の伝播時間である各階伝播時間により定まることを特徴とする、構造物の耐震性判定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、構造物の固有周期推定方法、構造物の耐震性判定方法、構造物の固有周期推定システム及び構造物の耐震性判定システムに関する。
【背景技術】
【０００２】
建物の地震動に対する損傷状況を判定するための耐震性判定システムとして、複数の加速度センサを建物の複数の階に設置して、構造躯体の損傷、使用安全性、修復可能性などを適切に評価する耐震性判定システムが提案されている（特許文献１）。この耐震性判定システムでは、地震が発生した場合に、加速度センサの出力から層間変形角を算出し、算出された層間変形角を用いて予め設定された複数の基準に従って建物の損傷度を判定している。しかしながら、地震動を受けた後の建物の固有周期を用いた累積的損傷を考慮した評価は行われていなかった。
【０００３】
また、近年の巨大地震に対する対策についての関心の高まりを受け、地震後の建物健全性評価を目的とするモニタリング技術に対する期待が高まっている。そのようなモニタリング技術がこれまでにも複数提案されている（非特許文献２）。しかしながら、いずれの提案も地震動を受けた後の建物の固有周期を効率よく的確に推定し、推定された固有周期を用いることで累積的損傷を考慮する実用的な健全性評価方法ではなかった。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−２５４２３９公報
【非特許文献】
【０００５】
【非特許文献１】白石理人、外２名、「モニタリング技術の現状と将来展望」、２０１６年日本建築学会大会（九州）、日本建築学会、２０１６年８月、「将来の大地震に備える強震観測とモニタリング」ｐ．２５〜３５
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は、地震動を受けた後の構造物についての実用的な固有周期推定方法及びこれにより推定された固有周期を用いた構造物の耐震性判定方法を提供することを目的とする。また、本発明は、地震動を受けた後の構造物についての実用的な固有周期推定システム及びこれにより推定された固有周期を用いた構造物の耐震性判定システムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は上述の課題の少なくとも一部を解決するためになされたものであり、以下の態様または適用例として実現することができる。
【０００８】
［適用例１］
本適用例に係る構造物の固有周期推定方法は、
地震動を受ける構造物の加速度データから周波数ごとの振幅データを算出し、
前記振幅データから前記構造物の一番下の階に対する一番上の階のスペクトル比を算出するに際し、分母にパワースペクトルの平均値に基づく値を加算して前記スペクトル比を算出し、
前記スペクトル比から増幅率が一番大きな固有周波数を選択し、
前記固有周波数から一次固有周期を算出することを特徴とする。
【０００９】
本適用例によれば、パワースペクトルの平均値に基づく値を分母に加算したスペクトル比を用いることにより、実用的な方法で地震動を受けた後の構造物の一次固有周期を推定することができる。
【００１０】
［適用例２］
上記適用例に係る構造物の固有周期推定方法において、
前記固有周波数における前記スペクトル比が予め設定されたしきい値以上である場合に、前記固有周波数から前記一次固有周期を算出し、
前記しきい値は、２〜１０に設定することができる。
【００１１】
本適用例によれば、スペクトル比をしきい値以上とすることで、構造物の共振時のスペクトル比に絞り込むことができ、振幅データにおける誤差を排除して、一次固有周期を算出することができる。
【００１２】
［適用例３］
上記適用例に係る構造物の固有周期推定方法において、
前記固有周波数が前記スペクトル比から得られたピーク周波数のうち最も低周波数側である場合に、前記固有周波数から前記一次固有周期を算出することができる。
【００１３】
本適用例によれば、一次固有周期に対応する固有周波数は最も低周波数側にあるため、選択された固有周期から一次固有周期を確実に算出することができる。
【００１４】
［適用例４］
上記適用例に係る構造物の固有周期推定方法において、
前記スペクトル比は、下記式（１）で算出されることができる。
【００１５】
【数１】
2019158731.tif 000003
【００１６】
本適用例によれば、スペクトル比を用いて適切な一次固有周期を得ることができる。
【００１７】
［適用例５］
本適用例に係る構造物の耐震性判定方法は、
上記適用例に係る構造物の固有周期推定方法で得られた前記一次固有周期に対応する伝達関数を作成し、
前記伝達関数に予め設定した設定地震動条件を入力して、地震応答解析をすることで前記設定地震動条件に対する前記構造物の層間変形角を演算し、
前記層間変形角が予め設定した許容値を超えているか否かで前記構造物の耐震性を判定することを特徴とする。
【００１８】
本適用例によれば、地震動を受けた後の構造物の累積損傷を考慮した耐震性を判定することができる。
【００１９】
［適用例６］
上記適用例に係る構造物の耐震性判定方法において、
前記伝達関数は、
前記振幅データにおける前記一次固有周期に対応する前記振幅データの前記一番上の階に対する加速度センサの設置階のスペクトル比を用いて前記一番上の階から前記設置階への地震動の伝播時間を求め、さらにそれを線形補間して求まる、前記一番上の階から前記構造物の各階への地震動の伝播時間である各階伝播時間により定まることができる。
【００２０】
本適用例によれば、一次固有周期に対応するように求められる各階伝播時間を用いて伝達関数を作成することで、多層階の構造物にも対応することができる。
【００２１】
［適用例７］
本適用例に係る構造物の固有周期推定システムは、
構造物の高さ方向における複数の位置の加速度を取得する加速度センサから取得された加速度に基づいて演算を行う演算部を含み、
前記演算部は、
地震動を受ける構造物の加速度データから周波数ごとの振幅データを算出し、
前記振幅データから前記構造物の一番下の階に対する一番上の階のスペクトル比を算出するに際し、分母にパワースペクトルの平均値に基づく値を加算して前記スペクトル比を算出し、
前記スペクトル比から増幅率が一番大きな固有周波数を選択し、
前記固有周波数から一次固有周期を算出することを特徴とする。
【００２２】
本適用例によれば、パワースペクトルの平均値に基づく値を分母に加算したスペクトル比を用いることにより、実用的な方法で地震動を受けた後の構造物の一次固有周期を推定することができる。
【００２３】
［適用例８］
上記適用例に係る構造物の固有周期推定システムにおいて、
前記固有周波数における前記スペクトル比が予め設定されたしきい値以上である場合に、前記固有周波数から前記一次固有周期を算出し、
前記しきい値は、２〜１０に設定することができる。
【００２４】
本適用例によれば、スペクトル比をしきい値以上とすることで、構造物の共振時のスペクトル比に絞り込むことができ、振幅データにおける誤差を排除して、一次固有周期を算出することができる。
【００２５】
［適用例９］
上記適用例に係る構造物の固有周期推定システムにおいて、
前記固有周波数が前記スペクトル比から得られたピーク周波数のうち最も低周波数側である場合に、前記固有周波数から一次固有周期を算出することができる。
【００２６】
本適用例によれば、一次固有周期に対応する固有周波数は最も低周波数側にあるため、選択された固有周期から一次固有周期を確実に算出することができる。
【００２７】
［適用例１０］
上記適用例に係る構造物の固有周期推定システムにおいて、
前記スペクトル比は、下記式（１）で算出されることができる。
【００２８】
【数１】
2019158731.tif 000004
【００２９】
本適用例によれば、スペクトル比を用いて適切な一次固有周期を得ることができる。
【００３０】
［適用例１１］
本適用例に係る構造物の耐震性判定システムは、
前記演算部が、
上記適用例に係る構造物の固有周期推定システムで得られた前記一次固有周期に対応する伝達関数を作成し、
前記伝達関数に予め設定した設定地震動条件を入力して、地震応答解析をすることで前記設定地震動条件に対する前記構造物の層間変形角を演算し、
前記層間変形角が予め設定した許容値を超えているか否かで前記構造物の耐震性を判定することを特徴とする。
【００３１】
本適用例によれば、地震動を受けた後の構造物の累積損傷を考慮した耐震性を判定することができる。
【００３２】
［適用例１２］
本適用例に係る構造物の耐震性判定システムにおいて、
前記伝達関数は、
前記振幅データにおける前記一次固有周期に対応する前記振幅データの前記一番上の階に対する加速度センサの設置階のスペクトル比を用いて前記一番上の階から前記設置階への地震動の伝播時間を求め、さらにそれを線形補間して求まる、前記一番上の階から前記構造物の各階への地震動の伝播時間である各階伝播時間により定まることができる。
【００３３】
本適用例によれば、一次固有周期に対応するように求められる各階伝播時間を用いて伝達関数を作成することで、多層階の構造物にも対応することができる。
【発明の効果】
【００３４】
本発明によれば、地震動を受けた後の構造物についての実用的な固有周期推定方法及びこれにより推定された固有周期を用いた構造物の耐震性判定方法を提供することができる。また、本発明によれば、地震動を受けた後の構造物についての実用的な固有周期推定システム及びこれにより推定された固有周期を用いた構造物の耐震性判定システムを提供することができる。
【図面の簡単な説明】
【００３５】
【図１】本実施形態に係る構造物の耐震性判定システムの概要を示す図である。
【図２】加速度センサから得られる加速度データの一例を示すグラフである。
【図３】本実施形態に係る構造物の耐震性判定システムの構成を示すブロック図である。
【図４】本実施形態に係る構造物の耐震性判定方法のフローチャートである。
【図５】本実施形態に係る構造物の耐震性判定システムで得られる地震応答解析モデルの一例を示す図である。
【図６】地震応答解析モデルにおける各階伝播時間の線形補間を説明する概念図である。
【図７】一次固有周期の推定方法のフローチャートである。
【図８】従来のスペクトル比と実施例１の一次固有周期の推定方法の修正スペクトル比とを示すグラフである。
【図９】実施例２の一次固有周期の推定方法を用いて一次固有周期を算出した結果を示すグラフである。
【発明を実施するための形態】
【００３６】
以下、本発明の好適な実施形態について、図面を用いて詳細に説明する。なお、以下に説明する実施形態は、特許請求の範囲に記載された本発明の内容を不当に限定するものではない。また、以下で説明される構成の全てが本発明の必須構成要件であるとは限らない。
【００３７】
１．構造物の耐震性判定システム
図１〜図３を用いて、構造物１０の耐震性判定システム２０の概要について説明する。図１は本実施形態に係る構造物１０の耐震性判定システム２０の概要を示す図であり、図２は加速度センサから得られる加速度データの一例を示すグラフであり、図３は本実施形態に係る構造物１０の耐震性判定システム２０の構成を示すブロック図である。なお、以下、構造物１０の耐震性判定システム２０は、単に「システム２０」という。
【００３８】
図１に示すように、システム２０は、構造物１０に設けられる。システム２０は、構造物１０の固有周期推定システム３０を含む。システム２０は、構造物１０に設置された加速度センサ２２から取得した構造物１０の高さ方向における複数の位置の加速度データに基づいて演算を行う演算部２６を含むモニタリングサーバ２４と、を含む。構造物１０の固有周期推定システム３０は、システム２０の演算部２６の一部の処理を行うものである。システム２０は、地震動を受けた後の構造物１０の累積損傷を考慮した耐震性を判定することができる。判定結果は、記憶部２８に保存され、例えば図示しないディスプレイに表示し、又はアラームとして出力等してもよい。
【００３９】
構造物１０は、土地に定着する人工物である工作物であって、ビル等の建築物を含む。図１における構造物１０は、複数階を有する多層建築物である。構造物１０は、１階（１Ｆ）から屋上階（ＲＦ）までの複数階を有し、図１では３階から屋上階から５階下の階までを省略して示す。
【００４０】
１−１．加速度センサ
加速度センサ２２は、地震動を受ける構造物１０の高さ方向における複数の位置に設けられ、例えば、構造物１０の複数の階に設置される。各階の地震時の変位を求めるためには全ての階に加速度センサ２２が設けられることが望ましいが、構造物１０における設置スペースや設置コストとの関係で加速度センサ２２が設置されない階が存在してもよい。例えば、図１では屋上階（ＲＦ）から３階下の階（ＲＦ−３Ｆ）には加速度センサ２２が
設置されていない状態を示す。
【００４１】
加速度センサ２２は、構造物１０の一番下の階（図１では１階であるが、地下階がある場合にはその地下階）と一番上の階（図１では屋上階）に少なくとも配置されることが望ましい。これは、固有周期の推定精度を高めるためである。例えば、屋上階に加速度センサ２２を設置できない場合には、できるだけ屋上階に近い階に加速度センサ２２を設置することが望ましい。
【００４２】
加速度センサ２２は、構造物１０が地震動を受けるときの構造物１０の加速度を測定するものである。加速度データｕ（ｔ）は、図２に示すような時間（ｔ）に対する加速度（ｕ）のデータであり、各階の床面に対して平行な２軸の加速度データを取得する。加速度センサ２２は、構造物１０の設置階のフロア又はフロア付近の加速度を測定することが好ましい。加速度センサ２２としては、例えば、ＭＥＭＳ（Ｍｉｃｒｏ  Ｅｌｅｃｔｒｏ  Ｍｅｃｈａｎｉｃａｌ  Ｓｙｓｔｅｍｓ）技術を用いたものを採用することができる。加速度センサ２２は、例えばフロアに平行な２軸の加速度を測定できる。
【００４３】
１−２．演算部
モニタリングサーバ２４は、構造物１０に設置された複数の加速度センサ２２と接続されており、加速度センサ２２からの電気信号を受信することができる。モニタリングサーバ２４は、図示しないＣＰＵ（中央演算処理装置）、ＲＯＭ、ＲＡＭ等のメモリやハードディスク装置等の記憶装置、外部装置との通信を行う通信インターフェース等を備える。図３に示すように、モニタリングサーバ２４は、少なくとも演算部２６と記憶部２８とを含む。演算部２６はＣＰＵやＲＡＭ等から構成することができ、記憶部２８はハードディスク装置等の記憶装置から構成することができる。
【００４４】
図１においてシステム２０は構造物１０内にモニタリングサーバ２４が設置されているが、図示しない別の構造物にモニタリングサーバ２４を設置して加速度センサ２２等と有線又は無線により通信してもよい。
【００４５】
演算部２６は、地震動における加速度センサ２２の出力に基づいて、その地震動を受けた構造物１０の一次固有周期を算出する。演算部２６は、算出される一次固有周期に対応する伝達関数を作成する。演算部２６は、作成される伝達関数に予め設定した設定地震動条件を入力して、地震応答解析をすることでその設定地震動条件に対する構造物１０の層間変形角を演算する。演算部２６は、演算される層間変形角が予め設定した許容値を超えているか否かで構造物１０の耐震性を判定する。
【００４６】
構造物１０は、新築当初の耐震性が経時的に徐々に変化する。例えば、構造物１０が比較的強い地震動を受けた場合には、構造物１０の一部が塑性変形を生じることもある。また、構造物１０の構造材が経時的に劣化していくこともある。したがって、構造物１０の正確な耐震性の判定は、現状の構造物１０の状況を反映したものでなければならない。そこで、システム２０は、地震発生直後の構造物１０の状況を把握するために、地震時における構造物１０の固有周期から構造物１０の耐震性を推定することとしている。
【００４７】
２．構造物の耐震性判定方法
本実施形態に係る構造物の耐震性判定方法は、構造物の固有周期推定方法で選択された前記一次固有周期に対応する伝達関数を作成し、前記伝達関数に予め設定した設定地震動条件を入力して、地震応答解析をすることで前記設定地震動条件に対する前記構造物の層間変形角を演算し、前記層間変形角が予め設定した許容値を超えているか否かで前記構造物の耐震性を判定することを特徴とする。
【００４８】
図１〜図６を用いて、構造物１０の耐震性判定方法について説明する。図４は、本実施形態に係る構造物１０の耐震性判定方法（以下「耐震性判定方法」という）のフローチャートであり、図５は本実施形態に係る構造物の耐震性判定システムで得られる地震応答解析モデルの一例を示す図であり、図６は地震応答解析モデルにおける各階伝播時間の線形補間を説明する概念図である。
【００４９】
図４に示す耐震性判定方法は、図１〜図３を用いて説明したシステム２０を採用することができる。
【００５０】
システム２０を起動し、耐震性判定を開始する。
【００５１】
Ｓ１０：地震が発生すると、モニタリングサーバ２４は、地震動における構造物１０の高さ方向における複数の位置の加速度センサ２２から加速度データｕ（ｔ）の取得を開始し、地震が終了するまで加速度データｕ（ｔ）を取得し、記憶部２８に記録する。
Ｓ２０：モニタリングサーバ２４の演算部２６は、取得した加速度データｕ（ｔ）に基づいて当該地震動を受けた構造物１０の一次固有周期の推定値を算出する。一次固有周期は、構造物１０の固有周期のうち一番長い周期を有するものである。固有周期は、構造物１０の自由振動のときの周期であり、構造物１０に固有な値を取るものである。構造物１０の一次固有周期は、構造物１０が地震動を受けて損傷を受けた場合には、構造物１０の設計時の固有周期計算により得られたものとは異なることになる。構造物１０の一次固有周期は、地震動を受ける構造物１０の加速度センサ２２からの加速度データに基づいて算出する。一次固有周期の推定は、後述するスペクトル比を用いた一次固有周期の推定方法を用いて行うことができる。
【００５２】
Ｓ３０：一次固有周期に対応する地震応答解析モデル８０（例えば図５）の伝達関数を作成する。図５を用いて、伝達関数の算出の作成について、説明する。図５は、本実施形態に係るシステム２０で得られる地震応答解析モデル８０の一例を示す図である。地震応答解析モデル８０は、演算部２６が算出した一次固有周期に対応する多質点系モデルである。多質点系モデルは、構造物１０の重量分布を各階の床面に重量が集中していると考える。図５では、一次モードにおける共振時（一次固有周期で揺れた時）の構造物１０の変形を表している。地震応答解析モデル８０は，構造物１０がある地震動を受けるとき、各階における応答変位を解析することができる。
【００５３】
地震応答解析モデル８０により求まる構造物１０の一番上の階の応答変位に対する各階の応答変位の比は、周波数領域において後述の式（４）の通りに表される。式（４）を利用することにより、構造物１０がある地震動を受けるとき、各階における応答変位を求めることができる。この式の計算は、例えば、市販されている表計算プログラム「ＭＩＣＲＯＳＯＦＴ  ＥＸＣＥＬ」（マイクロソフト社の登録商標）や、カシオ計算機社の高機能関数電卓「ｆｘ−９９３ＥＳ」などを用いて実施することができる。
【００５４】
伝達関数は、演算部２６で作成する。演算部２６は、一次固有周期に対応するように構造物１０の一番上の階から加速度センサ２２が設置された複数の位置（例えば構造物１０の各階）への地震波の伝播時間を求めると共に、それを線形補間して、構造物１０の一番上の階から構造物１０の各階への地震波の伝播時間である各階伝播時間を求め、作成する。このように、一次固有周期に対応するように求める各階伝播時間を用いて伝達関数を作成することで、多層階の構造物１０に対応することができる。なお、「一番上の階」は、以下の説明のように屋上階であることが望ましいが、屋上階に近い階としてもよい。
【００５５】
まず、演算部２６は、各階伝播時間を求める演算に必要な構造物１０の一次固有周期、および一次共振時（一次固有周期と地震動の周期が一致した時）における各階の応答変位
（水平方向の変位）を求める。
【００５６】
加速度センサ２２を設置した階における変位を屋上階の変位で割った値（各階の変位を屋上階の変位で正規化した値）をα０，α１，・・・，αｎとすると、屋上階から各階への地震波の伝播時間（各階伝播時間σｉ）は、下記式（２）で求めることができる。全ての階に加速度センサ２２が設置されていれば、全ての階の各階伝播時間が得られる。
【００５７】
【数２】
2019158731.tif 000005
【００５８】
この場合、屋上階のα０は、１になる。屋上階に加速度センサ２２が無い場合には、屋上階に最も近い階に設置された加速度センサ２２の値を用いることができる。なお、ｎは、加速度センサ２２の数である。
【００５９】
設置コストや設置場所の確保の問題から図１のように加速度センサ２２が設置されていない階がある場合が考えられる。そのような場合には、加速度センサ２２が設置されていない階は、各階伝播時間を線形補間して求めることができる。
【００６０】
図６は、地震応答解析モデル８０における各階伝播時間の線形補間を説明する概念図である。図６において、加速度センサ２２が設置されている階には質点（○）があり、加速度センサ２２が設置されていない階には質点がない（●）。図６に示すように、一番上の階（Ｎ）Ｆから数えてｍ０番目の階（Ｎ−ｍ０）Ｆ及びｍ２番目の階（Ｎ−ｍ２）Ｆに加速度センサ２２が設置されており、一番上の階（Ｎ）Ｆから数えてｍ１番目の階（Ｎ−ｍ１）Ｆに加速度センサ２２が設置されていない場合を考える。加速度センサ２２が設置されている階の各階伝播時間は、上記式（２）により求めることができるが、加速度センサ２２が設置されていない階では水平方向の変位は上記式（２）だけでは求めることができない。そこで、加速度センサ２２が設置されていない階の各階伝播時間は、上記式（２）から求められた各階伝播時間を図６及び下記式（３）のように線形補間して求める。構造物１０の総階数をＦとすると、各階伝播時間は下記式（３）を用いて例えばＳ１，Ｓ２，・・・，ＳＦとして得られる。
【００６１】
【数３】
2019158731.tif 000006
【００６２】
次に、各階伝播時間から一番上の階（屋上階）の応答に対する各階の応答の比である、伝達関数を求める。伝達関数は、下記式（４）を用いて求めることができる。下記式（４）は、一番上の階（ＮＦ）から数えてｆ番目の階の応答（揺れ）を求めることができる。下記式（４）では、各階伝播時間をＳｆで表すが、加速度センサ２２が設置されている階
について上記式（２）で得られるσｉに等しい。減衰定数（ｈ）は、一般的な建物で０．０２〜０．０５であるが、構造物１０の健全性を判定するために、例えば０．０１に設定することができる。
【００６３】
【数４】
2019158731.tif 000007
【００６４】
Ｓ４０：伝達関数に予め設定した設定地震動条件（例えばＬ１地震動条件）を入力して、地震応答解析をすることで設定地震動条件に対する構造物１０の層間変形角を演算する。
【００６５】
演算部２６は、作成された伝達関数に予め設定した設定地震動条件を入力することでその設定地震動条件に対する構造物１０の層間変形角を演算する。
【００６６】
層間変形角は、地震などの横揺れによって住宅などの建築物が変形する時、各階の床と真上または真下の床との、水平方向における変形の角度である。
【００６７】
設定地震動条件は、構造物１０の耐震性を判定する基準となる地震動として適当な条件を設定することができるが、建物の耐震設計における地震応答解析に用いられる「レベル１の地震動（以下「Ｌ１地震動」という）」を用いることができる。Ｌ１地震動は２５カイン（ｃｍ／ｓ）以上で基準化した地震波であり、建物の供用期間中複数回受けることが想定される規模の地震であるため、Ｌ１地震動に対して、主要構造体は弾性範囲内で応答しなければならない。
【００６８】
地震応答解析は、地震動に対して、建物等の構造物等の各部がどのような力を受けたり変形したりするかを検討するために、構造物等を適切な解析モデルに置き換え、相互作用を考慮した上で、設計用の地震動を入力してコンピューターで計算し、地震によって構造物等の各位置が受ける力と揺れの大きさを算出する解析法である。ここでは、伝達関数にＬ１地震動の条件を入力して解析を行う。
【００６９】
Ｌ１地震動の条件を入力した場合の各階変位のフーリエスペクトルは、下記式（５）で演算することができる。すなわち、構造物１０の１階部分にＬ１地震動が入力された場合に、Ｌ１地震動は予め定められているので、伝達関数ＨＦ（ω）およびＨｆ（ω）が上記式（４）により求まれば、構造物１０の揺れによる各階の変位のフーリエスペクトルの値が下記式（５）により求められる。
【００７０】
【数５】
2019158731.tif 000008
【００７１】
上記式（５）により得られたＹｆ（ω）をフーリエ逆変換すれば、各階の変位の波形ｙｆ（ｔ）が求められる。ここで、ｔは時間である。
【００７２】
そして、ｙｆ（ｔ）について各層で差分をとった最大値が層間変形角になる。層間変形角は、下記式（６）で求めることができる。
【００７３】
【数６】
2019158731.tif 000009
【００７４】
Ｓ５０：演算部２６は、上記式（６）で得られた層間変形角が予め設定した許容値を超えているか否かで構造物１０の耐震性を判定する。予め設定した許容値は、各階ごとに設定することができる。許容値は、層間変形角が構造物１０の弾性変形可能な範囲内に設定されている。構造物１０の弾性変形可能な範囲とは、構造物１０に加えた力を除いたときに生じていた歪みが力を加える以前の状態に戻る範囲である。したがって、今後Ｌ１地震動が発生した状況においても、構造物１０に塑性変形は生じないと推測され、十分な耐震性を備えると判断できる。
【００７５】
耐震性の判定は、構造物１０において算出される層間変形角の最大値が予め設定した許容値を超えている場合には構造物１０の耐震性能が不足していると判定し、当該最大値が許容値を超えていなければ耐震性能に変化がないと判定する。
【００７６】
このように、本実施形態に係る耐震性判定方法によれば、地震動を受けた後の構造物１０の累積損傷を考慮した耐震性を判定することができる。
【００７７】
２−１．一次固有周期の推定方法
図４におけるＳ２０の一次固有周期の推定方法について説明する。地震動を受けた構造物１０は、いずれかの箇所に損傷を受けることで設計時または施工時の固有周期と異なる固有周期を有する可能性がある。本実施形態では、地震動を受けた後の構造物１０の一次固有周期を加速度データから推定することで、地震動を受けた後の構造物１０の累積損傷を考慮した現実に近い耐震性を判定することができる。
【００７８】
図２及び図７を用いて、一次固有周期の推定方法について説明する。図７は一次固有周期の推定方法のフローチャートである。
【００７９】
一次固有周期の推定方法は、地震動を受ける構造物１０の加速度データｕ（ｔ）から周波数ごとの振幅データを算出し（Ｓ３１０）、振幅データから構造物の一番下の階に対する一番上の階のスペクトル比を算出するに際し、分母にパワースペクトルの平均値に基づく値を加算してスペクトル比を算出し（Ｓ３２０）、スペクトル比から増幅率が一番大きな固有周波数を選択し（Ｓ３３０）、固有周波数から一次固有周期を算出する（Ｓ３５０）。一次固有周期の推定方法によれば、パワースペクトルの平均値に基づく値を分母に加算したスペクトル比を用いることにより、実用的な方法で地震動を受けた後の構造物１０の一次固有周期を推定することができる。上述した構造物１０の固有周期推定システム３０の演算部２６は、少なくともＳ３１０，Ｓ３２０，Ｓ３３０及びＳ３５０の処理を行う。
【００８０】
図７に示すように、一次固有周期の推定方法を開始すると、Ｓ３１０〜Ｓ３８０の処理を実行して適切な一次固有周期を得ることができる。
【００８１】
Ｓ３１０：モニタリングサーバ２４の演算部２６は、地震動を受ける構造物１０の加速度データｕ（ｔ）から周波数ごとの振幅データを算出する。周波数ごとの振幅Ｕ（ω）は、地震動における構造物１０の加速度データｕ（ｔ）をフーリエ変換することで得られる振幅のフーリエスペクトルである。加速度データｕ（ｔ）は、図２を用いて上述したので説明を省略する。
【００８２】
Ｓ３２０：演算部２６は、Ｓ３１０で得られた振幅データから構造物１０の一番下の階（例えば１階（１Ｆ））に対する一番上の階（例えば屋上階（ＲＦ））の修正スペクトル比Ｇ＊（ω）を算出する。修正スペクトル比Ｇ＊（ω）は、下記式（１）で算出する。下記式（１）では、従来のようなスペクトル比と異なり、分母にパワースペクトルの平均値Ｐに基づく値（ａＰ）を加算して修正スペクトル比Ｇ＊（ω）を算出する。従来のスペクトル比と区別するために「修正スペクトル比」という。
【００８３】
【数１】
2019158731.tif 000010
【００８４】
修正スペクトル比Ｇ（ω）は、一番下の階と一番上の階とでスペクトル比をとることで、各周波数について構造物１０内でどれだけ振幅が増幅しているかがわかる。振幅の増幅率が最も大きな周波数に対応する周期が一次固有周期となる。上記式（１）において、補正係数ａはゼロ除算を回避するため適宜設定することができ、例えば０．５とすることができる。パワースペクトルの平均値Ｐは、フーリエスペクトルの振幅に共役複素数を掛けて全周波数帯域で合算して平均値をとったものである。上記式（１）を用いることにより、スペクトル比を用いて適切な一次固有周期を得ることができる。修正スペクトル比Ｇ＊（ω）は、従来の単純なスペクトル比（１階に対する屋上階のスペクトル比）に比べると振幅データのノイズによる誤差を排除することができる。
【００８５】
Ｓ３３０：演算部２６は、修正スペクトル比Ｇ＊（ω）から増幅率が一番大きな固有周波数（ピーク周波数ωｐ）を検出し、選択する。演算部２６は、ピーク周波数ωｐを記憶部２８に記憶する。
【００８６】
Ｓ３４０：演算部２６は、修正スペクトル比Ｇ＊（ω）があらかじめ設定されたしきい値ε１よりも大きいか否かを判定する。固有周波数における修正スペクトル比Ｇ＊（ω）が予め設定されたしきい値ε１よりも大きい場合に、固有周波数から一次固有周期を算出する。しきい値ε１は記憶部２８に保存されている。構造物１０が共振するのは、構造物１０の一番上の階と一番下の階とで少なくとも２倍以上に増幅するため、しきい値ε１は例えば２〜１０と設定することができる。修正スペクトル比Ｇ（ω）をしきい値以上とすることで、構造物の共振時以外の修正スペクトル比Ｇ（ω）を排除して、共振時の修正スペクトル比Ｇ（ω）に絞り込むことができる。これにより、振幅データにおける誤差を排除して、Ｓ３５０で一次固有周期Ｔ１を算出することができる。修正スペクトル比Ｇ＊（
ω）がしきい値ε１よりも大きい値であれば、Ｓ３５０の処理に進む。修正スペクトル比Ｇ＊（ω）がしきい値ε１よりも小さい値であれば、Ｓ３８０の処理に進む。
【００８７】
Ｓ３５０：演算部２６は、ピーク周波数ωｐから地震動を受けた後の構造物１０の一次固有周期Ｔ１を下記式（７）で求める。
【００８８】
【数７】
2019158731.tif 000011
【００８９】
Ｓ３６０：演算部２６は、記憶部２８から地震発生前の構造物１０の一次固有周期Ｔ０と、地震動を受けた後の構造物１０の一次固有周期Ｔ１（上記式（７）で算出される）との比をとり、予め設定しているしきい値ε２より大きい場合にＳ３７０の処理に進み、しきい値ε２以下であればＳ３８０の処理に進む。Ｓ３６０の処理によって、固有周波数が修正スペクトル比Ｇ（ω）から得られたピーク周波数のうち最も低周波数側である場合に、固有周波数から一次固有周期を選択する。一次固有周期Ｔ１に対応する固有周波数は最も低周波数側にあるため、選択された固有周期から一次固有周期Ｔ１を確実に算出することができる。
【００９０】
Ｓ３７０：演算部２６は、Ｓ３５０の処理で得られた一次固有周期Ｔ１を出力する。演算部２６から出力された一次固有周期Ｔ１は、例えば記憶部２８に地震動後の構造物１０の一次固有周期Ｔ１として保存され、図４のＳ３０における一次固有周期Ｔ１に対応する地震応答解析モデル８０の作成に用いられる。すなわち、一次固有周期Ｔ１は構造物１０の累積損傷を考慮した現実に近いものとなる。
【００９１】
Ｓ３８０：演算部２６は、地震動後の構造物１０の一次固有周期Ｔ１は、記憶部２８に保存されている構造物１０の一次固有周期Ｔ０に変化が無いと出力する。すなわち、地震動後の構造物１０の一次固有周期Ｔ１は地震動前の構造物１０の一次固有周期Ｔ０と同じあるため、図４のＳ３０以降の処理を行わない。
【００９２】
Ｓ３７０又はＳ３８０の処理の後、一次固有周期の推定方法が終了する。
【００９３】
２−２．地震応答解析モデルの作成
図４のＳ３０において、一次固有周期に対応する地震応答解析モデル８０（図５）を作成する。
【００９４】
地震応答解析モデル８０の作成は、修正スペクトル比Ｇ＊（ω）を用いて一次固有周期Ｔ１を選択すると判定された場合（図７のＳ３２０〜Ｓ３７０）に、一次固有周期Ｔ１に対応するように構造物１０の一番上の階（例えば屋上階（ＲＦ））から加速度センサ２２が設置された複数の位置への地震波の伝播時間を求めると共に、伝播時間に基づいて、構造物１０の一番上の階から構造物の各階への地震波の伝播時間である各階伝播時間を求め、各階伝播時間から一番上の階に対する各階の応答を求めることで行う。一次固有周期Ｔ１に対応するように求める各階伝播時間を用いて地震応答解析モデルを作成することで、多層階の構造物にも対応することができる。
【００９５】
伝達関数の作成については、上述した「１−２−２．伝達関数の作成」で説明したため、重複する説明を省略する。
【００９６】
構造物１０に加速度センサ２２が設置されていない階がある場合には、以下のように処理することができる。加速度データは、構造物１０の複数階に設置された加速度センサ２２から取得され、加速度センサ２２は、少なくとも構造物１０の一番上の階と一番下の階に設置され、構造物１０の一番上の階から加速度センサ２２が設置されている階への各階伝播時間は、下記式（２）で求められ、構造物１０の一番上の階から加速度センサ２２が設置されていない階への各階伝播時間は、下記式（２）で求める各階伝播時間を線形補間して求めることができる。線形補間をすることで、加速度センサが設置されていない階があっても、各階伝播時間を求めることができる。
【００９７】
【数２】
2019158731.tif 000012
【００９８】
上記式（３）におけるαｉとしては、一次固有周期の推定方法における図７のＳ３１０で地震加速度データｕ（ｔ）をフーリエ変換して得られた振幅データにおけるピーク周波数（一次固有周期Ｔ１に対応する周波数）における振幅の比を用いることができる。αｉは、一番上の階（例えば屋上階ＲＦ）の振幅に対する各階の振幅の比である。αｉは、一番上の階の変形量（ここでは振幅）で正規化された一次共振時の変形量として求めることができる。αｉは、下記式（８）で求めることができる。
【００９９】
【数８】
2019158731.tif 000013
【０１００】
加速度センサ２２が設置されていない階の応答計算については、上述した「１−２−２−１．加速度センサが設置されていない階の応答計算」で説明したため、重複する説明を省略する。
【実施例】
【０１０１】
（実施例１）
上記一次固有周期の推定方法について、従来のスペクトル比９２と図７のＳ３２０で求めた修正スペクトル比９４との違いを図８のグラフに示した。図８の縦軸は一番下の階１Ｆと一番上の階（屋上階）ＲＦとのスペクトル比であり、横軸は周波数である。従来のスペクトル比９２は同じ程度の高さのピークが複数あるが、修正スペクトル比９４は０．７Ｈｚ付近のピークを一次モードの周波数として確実に選択することができた。修正スペクトル比９４を用いることで、従来の単純なスペクトル比９２に比べると振幅データのノイ
ズによる誤差を排除することができることがわかった。
【０１０２】
（実施例２）
実際のビルに設置した加速度センサの加速度データを用いて、上記一次固有周期の推定方法によりビルの一次固有周期を推定した結果を図９に示した。図９において、縦軸は一次固有周期の推定方法により推定した一次固有周期であり、横軸は実際に地震が観測された日付である。このビルの初期の一次固有周期Ｔ０が約４秒であったので、図９で１秒付近に表れている推定された一次固有周期Ｔ１は、Ｔ１／Ｔ０が０．２５となり、図７のＳ３６０で設定したしきい値ε２が０．７であったので、一次固有周期を異常値として排除できた。
【０１０３】
本発明は、上述した実施形態に限定されるものではなく、さらに種々の変形が可能である。例えば、本発明は、実施形態で説明した構成と実質的に同一の構成（例えば、機能、方法、及び結果が同一の構成、あるいは目的及び効果が同一の構成）を含む。また、本発明は、実施形態で説明した構成の本質的でない部分を置き換えた構成を含む。また、本発明は、実施形態で説明した構成と同一の作用効果を奏する構成又は同一の目的を達成することができる構成を含む。また、本発明は、実施形態で説明した構成に公知技術を付加した構成を含む。
【符号の説明】
【０１０４】
１０…構造物、２０…システム、２２…加速度センサ、２４…モニタリングサーバ、２６…演算部、２８…記憶部、８０…地震応答解析モデル
【図１】
2019158731.tif 000016
【図２】
2019158731.tif 000017
【図３】
2019158731.tif 000018
【図４】
2019158731.tif 000019
【図５】
2019158731.tif 000020
【図６】
2019158731.tif 000021
【図７】
2019158731.tif 000022
【図８】
2019158731.tif 000023
【図９】
2019158731.tif 000024