ホーム   免震のすすめ / 高山 峯夫    第2回:免震構造の地震時応答の予測
  • 第1回:免震構造の実用化
  • 第2回:免震構造の地震時応答の予測
  • 第3回:やさしい免震層の設計
  • 第4回:熊本地震における建築物の被害
  • 第5回:熊本地震での免震構造の効果
  • 第6回:免震構造のすすめ
高山 峯夫

(たかやま みねお
/ Mineo Takayama)
福岡大学教授

 

< 略 歴 >

1960年福岡県生まれ。1982年福岡大学工学部卒業。1986年東京大学大学院修了後、福岡大学にて免震構造の実用化に取り組む。免震構造は大地震でも建物とその内部に被害は生じさせない。免震構造の実現に欠かせない積層ゴムの実大破壊実験、有限要素解析に取り組み、積層ゴムの優れた荷重支持性能のメカニズムを解明する。その成果により1998年日本建築学会奨励賞受賞。免震構造だけでなく地域やまちの災害リスク低減をする方策の研究にも力をいれている。著書に「4秒免震への道」(理工図書)、「耐震・制震・免震が一番わかる」(技術評論社)がある。専門は建築構造、免震構造。

 

< 著 書 >

4秒免震への道

耐震・制震・免震が一番わかる

第2回:免震構造の地震時応答の予測

はじめに

免震構造は免震層にエネルギー(変形)を集中させる明快な構造である。このため地震により建物に入力された全てのエネルギーを免震層で吸収することになる。入力エネルギーと吸収エネルギーのバランスを考慮することで免震層の応答変位や応答せん断力係数を予測する方法が確立されている。これを包絡解析法1)2)と呼ぶ。

包絡解析法の導出と適用性について紹介するとともに、エネルギーの釣り合いからわかる免震構造の応答特性についても紹介する。

入力エネルギーの計算

【図1】1自由度系振動モデル

建物下部にアイソレータとダンパーを配置した免震層を有する免震建物を対象とする。アイソレータとダンパーのみが地震エネルギーを吸収するものとし、上部構造は地震エネルギーの吸収を行わないものとする。この時、エネルギーの釣り合い式は次式で表される。

We(t)+Wp(t)=E(t)

ここで、We(t)は時刻tまでのアイソレータの弾性歪みエネルギー、Wp(t)はダンパーの消費エネルギー、E(t)は地震による入力エネルギーである。【図1】に示すような1自由度系モデルへ地動加速度¨zが作用する場合の地震入力エネルギーは、次式に示すように時々刻々の慣性力(M¨z)に応答の変位増分(dx)を乗じて積分することで計算できる。

E=M¨zdx

地震入力エネルギーEは建物質量Mで規準化されエネルギーの等価速度VEに変換して用いられることが多い。

VE=2EM

地震入力エネルギーの等価速度VEと周期の関係を示したものがエネルギースペクトルである。【図2】には粘性系や履歴系に対して求められたエネルギースペクトルの例を示す。エネルギースペクトルは一種の応答スペクトルであるが、通常の応答スペクトルでは減衰が大きくなるに従い応答値は減少するのに対し、エネルギースペクトルではある一定値に収束する。すなわち、減衰が0の振動系への入力エネルギーは非常に変動が大きいものの、粘性減衰が増える、あるいは塑性化が大きな系になれば入力エネルギーが平滑化されることがわかっている。

【図2】エネルギースペクトルの例1)

文献1)2)ではこの様なエネルギースペクトルに基づいて設定された設計用エネルギースペクトルも提案されている。設計用エネルギースペクトルは粘性減衰10%の弾性振動系の入力エネルギーとして求めるこができる。ほとんどの観測地震波のエネルギースペクトルは振動系の周期が長くなれば低下する傾向にある。しかし、地盤の種別や地震波の特性によっては、長周期領域での入力エネルギーがどのようになるかは予測できないため、設計用スペクトルとしては長周期領域で一定値、短周期領域で一定の勾配をもつバイリニアタイプのスペクトルが提案されている。このため免震構造の周期領域ではVEは一定値となり、地盤種別ごとにVE=120,150,200,300cm/sが提案されている。免震構造設計指針には様々な観測地震波に対するエネルギースペクトルが示されているので参照されたい。これまでの時刻歴応答解析でよく利用されてきたEL CENTRO波、八戸波を50cm/sに規準化した入力波の入力エネルギーの速度換算値は150cm/s程度、兵庫県南部地震の際の震源地近傍での観測波では、300cm/s程度であったとの報告もある。

入力エネルギー量を時刻歴応答解析により求め、入力エネルギーと時間との関係を描くことができる。弾性振動系への入力エネルギーは最大の変位応答が発生する時刻で最大の入力エネルギーを示す。一方、減衰の大きな系や塑性化が進行する系では入力エネルギーは単調増加となり、入力エネルギーの大部分が減衰あるいは塑性ひずみによる吸収エネルギーとなる。

(1)式は免震層が最大変形を示す時刻t=tmにおいても当然成立する。一般に塑性化の程度が大きな系では最大変形を示す時刻での入力エネルギー量E(tm)は最大値をとらないため、地震終了時(t=t0)の入力エネルギー総量E(t0)E(tm)E(t0)となる。従って、(1)式においてt=tmとして免震層の最大変位を予測する時、右辺をE(t0)で置き換えることは安全側の予測を行うことにつながる。しかし、地震波の特性によってはVE(t0)VE(tm)の差が大きくなる場合があり、応答予測の精度を低下させる原因となる。

免震層の応答予測式

【図3】アイソレータとダンパーの復元力特性

ここでは、【図3】に示すようにアイソレータが弾性、ダンパーが完全弾塑性型の復元力特性を有する場合を対象とする。免震部材には高減衰ゴム系積層ゴムや鉛プラグ型積層ゴムなどのようにダンパーを一体としたデバイスも用いられるが、これらの復元力特性は基本的にバイリニア型でモデル化されており、ここで紹介する方法に適用可能である。また様々な特性を有するデバイスが実用化されているものの、最も単純な復元力モデルを用いて検討しておくことは有効であり適用性も高い。

最大変形を生じるときのアイソレータとダンパーの吸収エネルギー量はそれぞれ次式で表せる。

We(tm)=12Kfδ2max

Wp(tm)=sQysδp

ここで、Kfはアイソレータの水平剛性、δmaxは免震層の最大変形であり、Weはアイソレータの弾性部分の履歴面積を算出している。また、sQyはダンパーの降伏耐力、sδpはダンパーの累積塑性変形量である。累積塑性変形量は塑性化領域の変形量(絶対量)を加算することで求められる。

sδp=δpi

Wpは、降伏体力と累積塑性変形量を乗じることで長方形の履歴面積を求めている。

累積塑性変形量と最大変形量の関係は、変数κを導入することで次式のように表せる3)

sδp=κ(δavesδy)

ここで、δaveは正側と負側の最大変形の平均値、sδyはダンパーの降伏変位、κは累積塑性変形量と最大変形量の比率である。δaveに比べ が十分に小さく、δaveδmaxに大きな差が無ければ(正側と負側に同じ程度の最大変形が発生していれば)、(3)式はsδp=κδmaxと近似できる。

この近似式と(2)式を(1)式に代入すれば、次式が得られる。

4π2(δmaxTf)2+2κgαsδmax=V2E

ここで、VEは地震入力エネルギーの等価速度、Mは建物総質量、gは重力加速度、Tfはアイソレータのみの水平剛性に基づいた周期で次式となる。

Tf=2πMKf

これをαsについて解けば、

αs=sQyMg=12κg(V2Eδmax4π2δmaxT2f)

が得られ、これにアイソレータのせん断力係数αfを加えることで、免震層の最大せん断力係数(以下、ベースシア係数)α1が次式のように求められる。

α1=αf+αs=4π2δmaxgT2f(112κ)+V2E2κgδmax

なお、アイソレータの最大せん断力係数αfは、次式で表せる。

αf=KfδmaxMg=4π2δmaxgT2f

同様に、免震層の最大変形量δmaxは次式で求められる。

δmax=κgαsT2f4π2[1+(2πVEκgαsTf)2+1]

応答予測図

(a)免震周期Tfを4秒とした場合

(b)入力エネルギーVEを150cm/sとした場合

(c)入力エネルギーVEを200cm/sとした場合

【図4】免震層のベースシア係数と最大変位の関係

【図4】にTfを4秒とした時のα1δmaxの関係を(4)式に基づいて示す。なお、κについては平均的に8を使うことができると言われており3)、【図4】でもκ=8を用いている。κ=8は、最大変形が起こるまでにダンパーが最大変形下での履歴ループに換算して2サイクル分のエネルギーを吸収したことに対応している。同図(a)は、Tf=4秒に固定した時のベースシア係数と免震層の変形の関係を描いている。下に凸の曲線は入力エネルギー量VEが一定の時の関係であり、右上がりの直線はダンパーの降伏せん断力係数が一定の時の関係である。VEの大きさによりベースシア係数が極小となる変形量は変化するものの、VE=150cm/sでは20~30cm程度、VE=300cm/sでは40cm程度の変形となる。なお、免震周期が短くなれば、免震層の変形は小さくなるものの、ベースシア係数が上昇することになる。様々なTfに関して同様の図を描いてみると免震構造の特性をより良く理解できると考える。同図(b)(c)ではVE=150cm/sおよび200cm/sと一定にした時の関係図である。免震周期が長くなるに従って、ベースシア係数が小さくなることがわかる。また、ダンパーの降伏せん断力係数、すなわち減衰が小さいと免震層の変形は大きくなり、減衰が大きすぎると変形は小さくなるもののベースシア係数が急激に大きくなる。このように想定される入力地震動の大きさに対して免震構造の応答を最小にするアイソレータの周期とダンパーの減衰量が存在することになる。

ベースシア係数の最小値α1minは、

α1min=15π4gVETf

となり、免震周期Tfと入力エネルギーの速度換算値VEにしか影響されないことがわかる。この時のアイソレータとダンパーのせん断力係数αoptsαoptf、及び変形量δoptmax

αopts=7π415gVETf

αoptf=8π415gVETf

δoptmax=1215πTfVE

となる。上式よりアイソレータとダンパーのせん断力の比率αf/αsが8/7とほぼ同じせん断力を負担している時が最適な状態であることが判る。【図5】には(6)式と(7)式から得られる免震層のベースシア係数と水平変形を示す。免震周期がのびるほど、応答変形は大きくなるものの、ベースシア係数は小さくなり、入力エネルギーの大きさによる変化も小さくなる。

【図5】免震層のベースシアの極小値とその時の変形量の関係

ベースシア係数の極小値をとる変位は周期が伸びる程、VEが大きくなる程増加するが、VE=300cm/sにおいても40cm程度の変形を見込んでおけば良いことがわかる。また免震建物をどの程度の地震入力エネルギー量に対して設計しているのか、又設計されているのかを次式により逆算することも性能を表す1つの指標になると思われる。

VE=g(2καs+αf)δmax

予測式の検証

エネルギーの釣り合いによる応答予測の精度を確認するために1質点系モデルによる応答解析を実施した。解析モデルの周期Tfは1~6sec、降伏変位sδyは1cmとした。入力地震波は9種類で、El Centro(NS)、Taft(EW)、八戸(NS)は最大速度を50cm/sに規準化している。その他のMexico地震やNorthridge地震などの観測波は原波形を用いた。

解析結果より、(3)式で定義されるκは周期の影響をあまり受けていないが、入力波の特性の影響を受け2~15程度まで変動する。大部分の場合κは8以下となっている。一方、VE(tm)/VE(t0)は周期3sec以上では0.5~0.8程度であり、κが小さな値をとる場合ほど低い値を示している。なお、入力エネルギーの等価速度VE(t0)は周期3sec以上ではほぼ一定値を示し、100~300cm/sの範囲にあった。

【図6】応答解析での最大変位と予測値の比較

【図6】には免震層最大変位と(5)式による予測値の関係を示す。(5)式の適用ではκ=8VE=VE(t0)とした。なお、ここでのt0は応答解析終了時としている。同図より、特にSylmar波(Northridge地震)の場合に予測値が小さくなっている以外は、予測値はほぼ安全側の予測になっていることが判る。

【図7】κによる応答予測の変動

※クリックで大きな画像が開きます

エネルギーの釣り合いに基づいた応答予測では、実際の応答がκ=8である場合でもVE(tm)/VE(t0)の値によっては予測値の精度は変化する。Northridge地震のような地震波の場合、κは4程度でVE(tm)/VE(t0)も0.8~1.2と大きかったために、予測値の精度が低下したと考えられる4)κを4,8,16と変化させた場合の応答値の変動を【図7】に示す。同図より、κが8よりも大きければ応答値は小さくなり、κが小さければ応答値は大きくなることが分かる。想定する地震動の特性によっては応答の予測値が変動することを考慮しておく必要がある。

時刻歴応答解析では個々の地震波の特性に応じて応答結果が変化するが、エネルギーの釣り合いによる応答予測では幾つかの仮定に基づく多少の変動はあるものの、包括的に応答を評価することが可能であり、免震建物(通常の建物も含め)の応答特性を誤ることなく評価できる手法であると言える。

ここで紹介したエネルギーの釣り合いに基づいた地震時応答予測手法は、免震建物の最大応答を簡単に求めることができ、さらに免震周期やダンパーの量に応じて応答値がどのように変動するかを把握することを可能とする。こういった手法で時刻歴応答解析結果を検証したり、設計した建物がどの程度までの地震エネルギーを吸収できるのかをチェックすることに活用されることを期待したい。なお、免震構造設計指針(日本建築学会編、2013年)では、ここで紹介した手法を拡張して、オイルダンパーを使うような場合の応答予測についても示されているので参考にしてほしい。

参考文献

1)
秋山宏:建築物の耐震極限設計(第2版)、東京大学出版会、1980
2)
秋山宏:エネルギーの釣合に基づく建築物の耐震設計、技報堂出版、1999
3)
秋山宏:第1層エネルギー集中型柔剛混合構造の地震応答予測、日本建築学会論文報告集、第400号、1989.6
4)
高山、森田:エネルギーの釣り合いに基づいた免震構造の応答予測法に関する研究、日本建築学会九州支部研究報告、第36号、1997.3

>> 「第3回:やさしい免震層の設計」