芝浦工业大学地震防灾研究室.DOC

7．構造物の地震応答理論（6） 7-1　地震力による応答 　地震力による運動方程式は式(3.2)で与えられた． 　　　(3.2) は，地震加速度である．ここで，をg(t)と表すことにする．したがって，式(3.2)は下記のようになる． 　　　(7.1) 地震力はおもりの（－質量×加速度）であり，加速度は地震計で記録される．現在の地震計は，ほとんどがデジタル記録方式であり，0.01秒や0.005秒間隔と一定の間隔で記録される．たとえば，0.01秒間隔であると，得られるgは， g(0),g(0.01),g(0.02),g(0.03),??? となり，これらはg(tn)と表すことができる．ここで， tn = n⊿tであり，⊿tは，記録間隔である．さらに，これを簡単にgnと表すこともできる．したがって，nは時間に対応する． 　ここでは，もっとも簡単なx(t)の求め方を示す．中央差分法を使えば，は以下のように求まる． (7.2) 　(7.3) ここで，は時刻tnにおける加速度，速度，変位である．式(7.2)，(7.3)がそれぞれ加速度，速度の次元になっていることを確認すること．速度は微小時間における変位の変化率であり，加速度は速度の変化率である．したがって，加速度も，速度も式(7.2)，(7.3)のように変位の変化率で表される．また，これらは，常に式(7.1)を満たしている．式(7.1)の両辺をmで割って，(7.2),(7.3)を代入すれば次のようになる． 　(7.3) n=0,1,2,????? x0=x-1=0 は地震記録であるので既知である．x0 = x-1=0とすると，式(7.3)からが求まる．次のはと今求まったから求まる．これを繰り返すことに，順次が求まる．また，加速度，速度は，式(7.2)，(7.3)を使うことにより求めることができる． 補足１：を求める方法は，中央差分法の他にもたくさんある．よく使われている方法としては，Runge-Kutta法，線形加速度法，Newmarkのβ法， Wilsonのθ法などがある．これらの詳しいことは，他の参考書を勉強すること．なお，上記の中央差分法，Runge-Kutta法は陽解法，線形加速度法，Newmarkのβ法， Wilsonのθ法は陰解法に属する．陽解法は，現時点の物理量は過去の物理量だけで決めるが，陰解法では未来の物理量を仮定して，過去と未来の物理量が矛盾なく連続となるように未来の物理量を修正していく方法である．一般に，陰解法は，陽解法に比べ，計算量は多いが精度はよい．これは，陽解法は過去だけの物理量しか使わないため，誤差が蓄積されてしまうためである．陰解法では，未来の物理量を使うことにより，この誤差を抑えている． 7-2　応答スペクトルの定義 　応答スペクトルとは，１自由度系の最大応答値を縦軸に，横軸に固有周期にとって描いた曲線のことである．減衰比hは，パラメータとなる．したがって，hが変われば応答スペクトルの形状は変わる．一般に，h=0.05が使われることが多い．応答スペクトルには，加速度応答スペクトル， 速度応答スペクトル，変位応答スペクトルの３種類がある．それぞれは，以下の式で定義される． 加速度応答スペクトル: 速度応答スペクトル: 変位応答スペクトル: ここで，は固有周期，hは減衰比である．｜｜は絶対値である．加速度応答スペクトルでは絶対加速度をとることに注意すること． 7-3　応答スペクトルの例 図7-1　2000年鳥取県西部地震のK-NET観測点TTR007（震央距離：約13km）で得られた加速度記録波形のNS成分 （a）=0.1秒 （b）=0.5秒 図7-2　各固有周期における１自由度系の相対速度 (減衰比hは0.05とする) (解法：Wilsonのθ法) （c）=1.0秒 （d）=1.5秒 （e）=2.0秒 図7-2（続き）各固有周期における１自由度系の相対速度 (減衰比hは0.05とする) 　　　　　　　（解法：Wilsonのθ法） 図7-3　速度応答スペクトル（ｈ＝0.05） 図7-4　加速度応答スペクトル（ｈ＝0.05） 図7-5　変位応答スペクトル（ｈ＝0.05） 補足 K-NET：強震ネットワーク（防災科学技術研究所インターネット上で公開されている． 2000年鳥取県西部地震：上記で取り上げた地震は，2000年10月6日に発生した鳥取県西部地震（気象庁マグニチュードMj7.3，震源深さ11km）である．気象庁マグニチュードに関しては，1995年兵庫県南部地震（Mj7.3，震源深さ14km）と同レベルであり，共に浅発の内陸直下型の地震である．兵庫県南部地震に比べると，被害規模は遥かに小さいが，全国のK-NET観測点の約3割にあたる304地点で