JP2013152197A - 建物の地震被害損失の評価装置・評価方法・評価プログラム - Google Patents
建物の地震被害損失の評価装置・評価方法・評価プログラム Download PDFInfo
- Publication number
- JP2013152197A JP2013152197A JP2012014046A JP2012014046A JP2013152197A JP 2013152197 A JP2013152197 A JP 2013152197A JP 2012014046 A JP2012014046 A JP 2012014046A JP 2012014046 A JP2012014046 A JP 2012014046A JP 2013152197 A JP2013152197 A JP 2013152197A
- Authority
- JP
- Japan
- Prior art keywords
- ground
- earthquake
- loss
- response
- building
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims description 70
- 230000004044 response Effects 0.000 claims abstract description 248
- 238000004364 calculation method Methods 0.000 claims abstract description 197
- 238000001228 spectrum Methods 0.000 claims abstract description 147
- 230000001133 acceleration Effects 0.000 claims abstract description 100
- 230000003321 amplification Effects 0.000 claims abstract description 96
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 96
- 239000010410 layer Substances 0.000 claims abstract description 92
- 238000010276 construction Methods 0.000 claims abstract description 88
- 238000004458 analytical method Methods 0.000 claims abstract description 83
- 238000006073 displacement reaction Methods 0.000 claims abstract description 51
- 238000000605 extraction Methods 0.000 claims abstract description 49
- 239000002344 surface layer Substances 0.000 claims abstract description 27
- 230000010365 information processing Effects 0.000 claims abstract description 25
- 239000000284 extract Substances 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 165
- 230000008569 process Effects 0.000 claims description 53
- 230000001186 cumulative effect Effects 0.000 claims description 39
- 239000002689 soil Substances 0.000 claims description 21
- 238000005315 distribution function Methods 0.000 claims description 16
- 230000001419 dependent effect Effects 0.000 claims description 12
- 238000012546 transfer Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 description 26
- 101100182248 Caenorhabditis elegans lat-2 gene Proteins 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 238000013461 design Methods 0.000 description 6
- 239000002356 single layer Substances 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000013016 damping Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000008439 repair process Effects 0.000 description 3
- 230000000717 retained effect Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Images
Abstract
【解決手段】 地震抽出手段1により、震源断層およびその地震の情報をデータベース9より抽出し、距離減衰式に従い、建設地の地盤の基盤における加速度応答スペクトルを作成する。地盤増幅率計算手段2は、多層の地盤の情報から一様な表層地盤と基盤との2層にモデル化し、表層地盤固有の地震動増幅率を求める。応答解析手段7は、建物につきモデル化して耐力スペクトルを求め、前記加速度応答スペクトルと表層地盤増幅率とから地表面応答スペクトルを求め、これら耐力スペクトルと地表面応答スペクトルとの関係から建物の応答変位,応答加速度を求める。損失・PML計算手段8は、応答解析結果を用いてPML値を計算する。
【選択図】 図1
Description
また、簡易にPML値を求めるについても、多層の地盤に応じた表層地盤の地震動の増幅率を踏まえ、ある程度は精度が担保されたPML値を算出することが望まれる。しかし、このような多層の地盤を考慮してPML値を算出できる簡易な装置はない。
また、建物各層の荷重変形関係(構造計算結果)を踏まえた建物モデルを用いた応答解析を踏まえ、PML値を算出できる簡易な装置はない。
この発明の他の目的は、建物各層の荷重変形関係を踏まえた建物モデルを用いた応答解析が行えて、より精度の良いPML値が簡易に求められるようにすることである。
前記入力部1は、評価対象となる建物の建設地の位置および地盤の情報、並びに建物の構造計算結果および部位毎の価格の情報の入力を受付ける。
前記地震抽出手段5は、建設地を基準とする対象範囲内に存在する各震源断層、並びにその震源断層で想定される地震のマグニチュード、年発生確率、震源深さ、および位置の情報をデータベースより抽出し、抽出した位置の情報から各震源断層と建設地との距離を計算して、定められた距離減衰式に従い、建設地の地盤の基盤における前記各震源断層で地震が生じた場合の建物周期ごとの最大加速度応答を示す加速度応答スペクトルを作成する。
前記地盤増幅率計算手段5は、建設地の地盤につき、入力された多層の地盤の情報から、表層地盤と基盤との2層に大別してモデル化し、表層地盤は一様な地質の1層または複数層の層構成とし、基盤の地震動に対して表層地盤の地表面に生じる地震動が増幅される率である表層地盤固有の地震動増幅率を定められた理論計算式に従って求める。
前記応答解析手段7は、多層からなる建物を1質点系にモデル化してそのモデルの荷重と変形の関係を示す耐力スペクトルを求め、前記基盤の加速度応答スペクトルと前記表層地盤固有の地震動増幅率とから、前記各震源断層で生じた地震による建設地の地表面における加速応答スペクトルと変位応答スペクトルより設定される地表面応答スペクトル(実施形態のデマンドスペクトルが該当)を求め、これら耐力スペクトルと地表面応答スペクトルとの関係から、応答解析結果として前記各地震による建物の応答変位,応答加速度を求める。
前記損失・PML計算手段8は、前記応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算し、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を例えばβ分布等の分布関数でモデル化して90%非超過損失を算定し、90%非超過損失の大きくなる地震より年発生確率の累積である累積確率を求め、累積確率が定められた値となる90%非超過損失であるPML値を計算する。
前記出力部4は前記損失・PML計算手段8の計算結果となるPML値を表示装置の画面に出力する。
なお、前記「一つの情報処理機器」とは、パーソナルコンピュータ等の機器のハードウェアを言い、機器本体の他に、USBメモリやハードディスクドライブ等の周辺機器が接続されたものを含む。
このように、地盤増幅率計算手段6が多層の地盤に応じた表層地盤増幅率を簡単な演算で求めることができるため、地盤増幅率計算手段6が簡易な構成で済み、そのため、地震抽出手段5、地盤増幅率計算手段6、応答解析手段7、および損失・PML計算手段8を含む評価装置を、ノート型等パーソナルコンピュータ等の一般的な一つの情報処理機器に設けることができる。
また、それぞれ前記処理を行う地震抽出手段5、地盤増幅率計算手段6、応答解析手段7、および損失・PML計算手段8を一つの情報処理機器に備えるため、建設地の位置、多層の地盤の情報、建物構造の情報があれば、地震抽出からPML値の計算までの一連の計算を1か所で纏めて行え、しかも多層の地盤に応じた地震動の増幅特性を考慮した比較的精度の良いPML値の計算が行える。
表層地盤増幅率の計算には、重複反射理論に基づく伝達関数を用いるようにしても良い。
この構成の場合、多層の地盤に応じた表層地盤増幅率を簡単な演算で求めることが実現できる。
土質を粘性土と砂質土との種類し、その割合で1層にモデル化するため、モデル化の演算が簡易に行える。
このように建物の各層の荷重と変形の関係から1質点系に集約したモデルを応答解析に用いることで、建物各層の荷重変形関係を踏まえた応答解析が行え、より精度の良いPML値を簡易に求めることができる。
前記入力過程S1は、評価対象となる建物の建設地の位置および地盤の情報、並びに建物の構造計算結果および部位毎の価格の情報の入力を受付ける。
前記地震抽出過程S2は、建設地を基準とする対象範囲内に存在する各震源断層、並びにその震源断層で想定される地震のマグニチュード、年発生確率、震源深さ、および位置の情報をデータベースより抽出し、抽出した位置の情報から各震源断層と建設地との距離を計算して、定められた距離減衰式に従い、建設地の地盤の基盤における前記各震源断層で地震が生じた場合の建物周期ごとの最大加速度応答を示す加速度応答スペクトルを作成する。
前記地盤増幅率計算過程S3は、建設地の地盤につき、入力された多層の地盤の情報から、表層地盤と基盤との2層に大別してモデル化し、表層地盤は一様な地質の1層または複数層の層構成とし、基盤の地震動に対して表層地盤の地表面に生じる地震動が増幅される率である表層地盤固有の地震動増幅率を定められた理論計算式に従って求める。
前記応答解析過程S4は、多層からなる建物を1質点系にモデル化してそのモデルの荷重と変形の関係を示す耐力スペクトルを求め、前記基盤の加速度応答スペクトルと前表層地盤固有の地震動増幅率とから、前記各震源断層で生じた地震による建設地の地表面における加速応答スペクトルと変位応答スペクトルより設定される地表面応答スペクトルを求め、これら耐力スペクトルと地表面応答スペクトルとの関係から、応答解析結果として前記各地震による建物の応答変位,応答加速度を求める。
前記損失・PML計算過程S5は、前記応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算し、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を例えばβ分布等の分布関数でモデル化して90%非超過損失を算定し、90%非超過損失の大きくなる地震より年発生確率の累積である累積確率を求め、累積確率が定められた値が定められた値となる90%非超過損失であるPML値を計算する。
前記出力過程S6は前記損失・PML計算過程S5の計算結果となるPML値を表示装置の画面に出力する。
なお、前記地震抽出過程S2と地盤増幅率計算過程S3とは、いずれを先に行っても、また並行して行っても良い。
前記入力手順R1は、評価対象となる建物の建設地の位置および地盤の情報、並びに建物の構造計算結果および部位毎の価格の情報の入力を受付ける。
前記地震抽出手順R2は、建設地を基準とする対象範囲内に存在する各震源断層、並びにその震源断層で想定される地震のマグニチュード、年発生確率、震源深さ、および位置の情報をデータベースより抽出し、抽出した位置の情報から各震源断層と建設地との距離を計算して、定められた距離減衰式に従い、建設地の地盤の基盤における前記各震源断層で地震が生じた場合の建物周期ごとの最大加速度応答を示す加速度応答スペクトルを作成する。
前記地盤増幅率計算手順R3は、建設地の地盤につき、入力された多層の地盤の情報から、表層地盤と基盤との2層に大別してモデル化し、表層地盤は一様な地質の1層または複数層の層構成とし、基盤の地震動に対して表層地盤の地表面に生じる地震動が増幅される率である表層地盤固有の地震動増幅率を定められた理論計算式に従って求める。
前記応答解析手順R4は、多層からなる建物を1質点系にモデル化してそのモデルの荷重と変形の関係を示す耐力スペクトルを求め、前記基盤の加速度応答スペクトルと前記表層地盤固有の地震動増幅率とから、前記各震源断層で生じた地震による建設地の地表面における加速応答スペクトルと変位応答スペクトルより設定される地表面応答スペクトルを求め、これら耐力スペクトルと地表面応答スペクトルとの関係から、応答解析結果として前記各地震による建物の応答変位,応答加速度を求める。
前記損失・PML計算手順R5は、前記応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算し、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を例えばβ分布等の分布関数でモデル化して90%非超過損失を算定し、90%非超過損失の大きくなる地震より年発生確率の累積である累積確率を求め、累積確率が定められた値が定められた値となる90%非超過損失であるPML値を計算する。
前記出力手順R6は前記損失・PML計算手順の計算結果となるPML値を表示装置の画面に出力する。
なお、前記地震抽出手順R2と地盤増幅率計算手順R3とは、いずれを先に行っても、また並行して行っても良い。
図2は、図1の演算処理部3の処理を示す流れ図である。図2における丸囲みの符号1〜4は、図1の地震抽出手段5、地盤増幅率計算手段6、応答解析手段7、損失・PML計算手段8に付した丸囲みの符号1〜4で示す手段で行う処理であることを示す。図2は入力情報の種類および各過程の処理結果の利用の関係が明確となるように各過程を図示している。各過程の詳細は図25に示す。
出力部4は、図22に示す画面を表示装置14aに出力し、この画面の中で、損失・PML計算手段8の計算結果となるPML値や、イベントリスク曲線、計算対象範囲、建物損失に影響を与える地震のリスト等を表示する。
(1)建設地と震源断層の緯度経度を3次元直交座標に置換
図6に示すように、震源断層の断層原点Oを座標系の原点とし、X0を東方向、Y0 を北方向、Z0 を法線方向とした座標系で、建設地の位置を(x0 ,y0 ,z0 )で表す。この処理は、角震央距離から震央距離に変換するための係数p(地球の平均半径)を用いると、理科年表より以下の式で、断層原点の経度long、地理緯度lat2、建設地の経度long' 、地理緯度lat2` を座標系に変換できる。
なお、「震央距離」は、地球の大円に沿って測った、震央から観測点までの距離である。
「角震央距離」は角度のことであるため、次の計算式では角度に地球の平均半径を掛けて、長さの単位となる「距離」を求めている。震央距離の具体的は計算方法は、理科年表などに記載されている。
「角震央距離」は、震央及び観測点と地球の中心とを結ぶ直線の間の角度である。
y0 =±arccos(sin(lat2)・sin(lat2')+cos(lat2)・cos(lat2')・cos(long−long))・p
z0 =d
ただし、
p=(a2 b)1/3 (赤道半径(a=6378.137 Km)および、極半径(b=6356.752
Km))
d:断層上端深さ
図7の領域(1)は震源断層面を表している。例えば、観測点が、領域(1)にある場合は、震源深さzが断層最短距離となり、領域(3)の上にある場合は、(x2-W)^2+(y2-D)^2+Z2^2の平方が断層最短距離となり、領域(4)の上にある場合は、(y2-D)^2+(Z2)^2の平方が断層最短距離となる。同図、計算式の違いにより、領域を分けている。なお、符号(1)〜は、図では丸囲みの数字で示している。
図8に示すような工学的基盤における応答スペクトルは、例えば安中らの距離減衰式(非特許文献2)を用いる。周期T、減衰定数h=5%における工学的基盤面の地震動強さSab (T,h=5%)は以下の式より求められる。
log10 Sab (T,h=5%)
=Cm (T)M+Ch (T)H−Cd (T)log10 {R+0.344exp(0.653 M)}+C0 M:気象庁マグニチュード、H:断層中心深さ[km]、
R:断層最短距離[km]、
Cm (T),Ch (T),Cd (T),C0(T):回帰係数
なお、各回帰係数の数値は非特許文献2から読み取る等する。
「工学的基盤」は、建物の基礎として耐え得るような地盤〈N値50以上〉に相当するS波速度が300〜700m/s(およそ400m/s)の層を言う。
上記「基盤」は「工学的基盤」であり、上記工学的基盤の定義に該当する範囲で任意に定めた深さ位置を基盤の表面とした。基盤表面位置の定め方について、具体的には、例えば柱状図等からN値が50〜60以上で、支持層と考えることができる層とし、S波速度を計算すると、概ね工学的基盤(S波速度400m/s)に相当する。
歪依存特性とは、土質材料が振動などを受けた場合に歪が生じ(非線形化)、その歪量γに応じた剛性G及び減衰hの変化を示したもので非特許文献9で挙げられているように複数提案されている。後述の段落〔0044〕のh−γ曲線や、〔数4〕のh(γ)、段落〔0018〕のG(γ)/Goが歪依存特性に該当する。提案されている歪依存特性の多くは、粘性土と砂質土に分けて設定されている。そこでこの実施形態の評価装置で入力した多層の地盤を、粘性土と砂質土に分類し、その割合を用いて1層地盤の歪依存特性と作成している。
(1)表層地盤のモデル化
表層地盤における第i層の層厚・せん断波速度・密度をそれぞれhI ・VsI ・ρI とすると、1層に縮約された表層地盤(等価表層地盤と呼ぶ)の層厚hE ・せん断波速度VsE ・密度ρE は以下の式で求められる。
図9に示すような表層地盤の地震動増幅率Gsは、重複反射理論に基づく伝達関数を用いる。非線形時の複素円振動数p、インピーダンス比α、h−γ曲線より求められるhs(γ)を用いると、以下の式で求められる。
ここで、T:周期、
γ:一次固有周期TS でのせん断ひずみ
HS :表面地盤の層厚
VS (γ):せん断ひずみγに応じた表面地盤のS波速度(m/s)
歪依存特性(せん断ひずみγと剛性低下率G/Goの関係より、弾性時の表層地盤のS波速度Vs0を用いて、次の式で表される。
VS (γ)=Vs0・√G(γ)/Go
hs(γ):せん断ひずみγに応じた減衰係数
VSB:基盤のS波速度(m/S)
複素円振動数p、インピーダンス比αは、〔数4〕のp(T,γ)、α(γ)の式より求める。h-γ曲線は、先の歪依存特性のことで、hs(γ)が得られる。
なお、表層基盤の増幅率については、例えば、柴田明徳氏の「最新耐震構造解析」を参考に計算できる。
heS(γ)=1/(2・GS1(γ))
GS1(γ)=GS (TS ,γ)=1/(1.57・hS (γ)+α0 ・√G(γ)/G0 ) α0 =ρE ・VsE/ρB ・VsB (=インピーダンス比)
非線形時の1次固有周期TS に対応するせん断ひずみγを与えることで、周期TS における基盤の加速度応答スペクトルSaeB (T=TS)を算定することができる。基盤の加速度応答スペクトルをSaB(T=TS )とおき、せん断ひずみγを新たにγ・SaB/SaeB 等と修正して再度、SaeB (T=TS )を算定する。以上を繰り返して、SaB=SaeB となるような、せん断ひずみγを求める。
また、上記耐力スペクトルを求めるについては、構造計算結果、例えば建物各層の荷重変形関係を踏まえ、耐力スペクトルを作成する。
建物各層の荷重変形を作成する際の構造計算結果の入力方法は、例えば次の3種類から選択して入力可能とする。
・ プッシュオーバー(PushOver)解析結果の直接入力(推奨)。
・ 降伏点(Q,δ)(耐力、変位)の入力。
・ 保有耐力の余裕度などの入力
これら3種類の入力方法のうち、プッシュオーバー(PushOver)解析結果の直接入力が最も好ましい。
応答解析の手順の概略は次の通りである。上記3種の入力方法の中から任意の方法を選んで構造計算結果を入力する。多質点系から1質点系モデルに変換する。応答スペクトル法による地震応答解析を実施し、一質点系モデルの応答変位,応答加速度を求める。
(1)建物各層の荷重変形関係を作成
以下の3種類の中の任意の方法で各層の荷重変形関係を作成可能としている。
・ プッシュオーバー(PushOver)解析結果の直接入力(図11)をした場合の曲線の延長方法。
・ 降伏点(耐力・変位)の入力による曲線の作成方法〈図12〉。
・ 保有水平耐力の余裕度等の入力による曲線の作成方法〈図12〉。
図11は、解析結果の直接入力による建物各層の荷重変形関係を表しており、横軸は建物変位、縦軸は建物耐力を表している。解析結果の直接入力とは、SS2やBUS(商品名)など市販の構造計算ソフトウェアを用い、二次設計(PushOver解析)まで実施した際に出力可能な、CSV形式などで作成される外力を漸増させたステップごとの荷重(耐力)・変位の値を、この実施系他の評価装置に直接コピー・ペーストする操作を指す。
建物の各層の荷重変形関係から、1質点系に縮約した場合の建物の、荷重変形関係(=耐力スペクトル)を求める。なお、この処理は非特許文献6を参考にした処理である。
デマンドスペクトルは入力地震動の周期特性を示すもので、加速度応答スペクトルと変位応答スペクトルから求める。加速度応答スペクトルは、地震抽出手段5で求めた基盤の応答スペクトルに、地盤増幅率計算手段6で求めた表層地盤の増幅率を乗じたものとする。変位応答スペクトルは擬似とし、次の式で求める。
SD =SA/ω2
ここで、SD :変位スペクトル:SA :加速度応答スペクトル:ω:固有円振動数:である。
耐力スペクトルとデマンドスペクトルを同一のグラフに表示し、2つの曲線の交点を一質点系の応答変位・応答加速度とする。図16において、実線は耐力スペクトルを、各破線は各地震のデマンドスペクトルを示す。同図の交点が、一質点系の応答変位・応答加速度である。
ただし、交点が耐力スペクトルの塑性域となっている場合には、塑性化による建物減衰の増大を考慮し、デマンドスペクトルの低減を行う。低減したスペクトルとの交点により求まる塑性率μと設定した塑性率μが等しくなるまで収斂計算を行う。
SA (T,h=5%):表層地盤考慮後の加速度スペクトル、
なお、等価粘性減衰定数heqは、等価−自由度系の塑性率μを用いて、次の式で表される。
heq=γ(1−1/√μ)+0.05
多質点系の応答変位は、耐力スペクトル法での応答解析結果のステップ数を用いて、構造計算における増分解析(PushOver解析と同じ)の同一ステップ数の数値を各階での応答変位とする。各階の応答加速度Acci についてはAi分布を用いて次式で求める。
損傷確率計算部8aは、応答解析手段7の応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算する。
90%非超過損失算定部8bは、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を例えばβ分布等の分布関数でモデル化して90%非超過損失を算定する。
PML値計算部8cは、90%非超過損失の大きくなる地震より年発生確率の累積である累積確率を求め、累積確率が定められた値となる90%非超過損失であるPML値を計算する。
累積確率は、90%非超過損失の年発生確率の累積である。 PML値は、予想最大損失率(額)であり、累積確率が定められた値(1/475)となるような90%非超過損失を言う。累積確率が1/475となるような90%非超過損失であるとは、475年に1回生じる90%非超過損失確率である。
建物の構造形式に応じた既存のフラジリティ曲線を用いて、構造部材、非構造部材の損傷確率を計算する。フラジリティ曲線は、建物応答(もしくは地震動強さ)に対する損傷確率である。
イベントツリー解析を用い、地震毎に損失期待値と、部位毎の損失の割合を計算する。 β分布を用い、損失期待値を中央値として損失分布をモデル化し、90%非超過損失を算定する。
損失に影響の大きい地震から順に年超過確率の累積確率が1/475となる90%非超過損失をPML値とする。
(1)フラジリティ曲線を用いた損傷確率の計算
建物の地震応答に対する損傷確率を表した構造形式に応じたフラジリティ曲線を設定する。フラジリティ曲線は、図18のように一般的に対数正規分布関数で表され、応答変位・加速度に応じた関数を用いる。なお、このフラジリティ曲線は、非特許文献7を参考にした曲線である。
フラジリティ曲線や、中央値、ばらつき等のデータは、例えばデータベース10に登録しておき、損失・PML計算手段8によりデータベース10にアクセスして計算に必要なデータを得るようにしている。
各地震に対して図19に示すようなイベントツリー解析を行い、損失(損傷費用)の期待値と標準偏差を求める。なお、この解析は非特許文献8に示された方法である。応答変位に依存する損失としては、構造部材や、外壁などの非構造部材がある。一方、応答加速度に依存する損失としては、設備機器などの非構造部材がある。
修理可能な場合: DCI =CSt×DCRst,a,α+Cd ×DCRd,β+Ca ×DCRa ,γ
修理不可能な場合 DCi =E[CD |Fail ]
ここで、
α:構造部材の損傷モード、
β:変形依存非構造部材の損傷モード、
γ:加速度依存非構造部材の損傷モード、
DCR:損傷モードに対する部位ごとの損失
Cst:建物全体の構造部材の再調査費用
Cd:変位依存型非構造部材の費用
Ca:加速度依存型非構造部材の費用
E[CD |Fail ]:期待損傷費用
ある地震に対する損失を例えばβ分布等の分布関数でモデル化する。図20に示すように分布の確率密度関数の面積、つまり分布関数が0.9になるような損失(=90%非超過損失)を地震ごとに求める。
90%非超過損失を、対象とする全ての地震について計算する。図7に示すように、損失の大きい順に並べて、地震の年発生確率の累積を求める。累積確率はpi を地震の年発生確率として次式のように求める。
P=1−Πi(1-Pi)
90%非超過損失と累積確率をプロットしたものをイベントリスク曲線と呼び、累積確率が1/475となるような90%非超過損失〈率〉をPML値とする。
すなわち、抽出手順R2aは、建設地を基準とする対象範囲内に存在する各震源断層、並びにその各震源断層で生じる地震のマグニチュード、年発生確率、震源深さ、および震源の位置、大きさ、角度の情報をデータベース9より抽出する。基盤加速度応答スペクトル作成手順R2bは、抽出した位置の情報から各震源断層と建設地との距離を計算して、定められた距離減衰式に従い、建設地の地盤の基盤における前記各震源断層で地震が生じた場合の加速度応答スペクトルを作成する。
すなわち、地盤モデル生成手順R3aは、建設地の地盤につき、入力された多層の地盤の情報から、表層地盤と基盤との2層にモデル化する。表面地盤増幅率計算部6bは、基盤の地震動に対して表層地盤の地表面に生じる地震動が増幅される率である表層地盤固有
の地震動増幅率を定められた理論計算式に従って求める。
すなわち、建物モデル・耐力スペクトル生成手順R4aは、建物を1質点系にモデル化してそのモデルの荷重と変形の関係を示す耐力スペクトルを求める。地表面応答スペクトル生成手順R4bは、前記基盤の加速度応答スペクトルと前記表層地盤増幅率とから、前記各震源断層で生じる地震による建設地の地表面における加速応答スペクトルと変位応答スペクトルより設定される地表面応答スペクトルを求める。応答変位・応答加速度計算手順R4cは、これら耐力スペクトルと地表面応答スペクトルとの関係から、応答解析結果として前記各地震による建物の応答変位,応答加速度を求める。
すなわち、損失計算手順R5aは、前記応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算する。90%非超過損失算定手順R5bは、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を例えばβ分布等の分布関数でモデル化して90%非超過損失を算定する。PML値計算手順R5cは、90%非超過損失の大きくなる年発生確率の累積である累積確率を求め、累積確率が定められた値となる90%非超過損失であるPML値を計算する。
なお、このプログラム100の各手順の詳細は、評価装置の対応する各手段に説明した処理を行う手順である。
地震抽出過程S2は、抽出手順S2aと基盤加速度応答スペクトル作成手順S2bからなる。地盤増幅率計算過程S3は、地盤モデル生成過程S3aと表面地盤増幅率計算過程S3bからなる。前記応答解析過程S4は、建物モデル・耐力スペクトル生成過程S4a、地表面応答スペクトル生成過程S4b、および応答変位・応答加速度計算過程S4cからなる。前記損失・PML計算過程S5は、損失計算過程S5a、90%非超過損失算定過程S5b、およびPML値計算過程S5cからなる。
これらの各過程S1〜S6,S2a,S2b,S3a,S3b,S4a,S4b、S4c,S5a,S5b,およびS5cは、評価プログラム100の各手順につき対応する名称を付した各手順の実行により行われる過程である。
・建物の地盤モデルに応じた増幅率特性、かつ、建物各層の荷重変形関係(構造計算結果)を踏まえた建物モデルを用いた応答解析を実施可能。
・一連の計算を一つのシステムで実施することが可能。短時間で計算が完了。
・今後の新たな知見により、入力部や演算処理部、データベースが更新されても、入力データ処理部2による更新により、プログラム全体への影響を軽減できる。
・構造設計実施前にPML値を推定し、構造設計に役立てることが可能。
・一般的なPC(パーソナルコンピュータ)の性能で十分に動作可能。
・ノートPCで営業先でプレゼンツールとして使うことができる。
2…入力データ処理部
4…出力部
5…地震抽出手段
5a…地震抽出部
5b…基盤加速度応答スペクトル作成部
6…地盤増幅率計算手段
6a…地盤モデル生成部
6b…表面地盤増幅率計算部
7…応答解析手段
7a…建物モデル・耐力スペクトル生成部
7b…地表面応答スペクトル生成部
7c…応答変位・応答加速度計算部
8…損失・PML計算手段
8a…損失計算部
8b…90%非超過損失算定部
8c…PML値計算部
9,10…データベース
14…出力機器
14a…表示装置
100…評価プログラム
Claims (7)
- 入力部と、地震抽出手段と、地盤増幅率計算手段と、応答解析手段と、損失・PML計算手段と、出力部を一つの情報処理機器に備え、
前記入力部は、評価対象となる建物の建設地の位置および地盤の情報、並びに建物の構造計算結果および部位毎の価格の情報の入力を受付け、
前記地震抽出手段は、
建設地を基準とする対象範囲内に存在する各震源断層、並びにその各震源断層で想定される地震のマグニチュード、年発生確率、震源深さ、および位置の情報をデータベースより抽出し、抽出した位置の情報から各震源断層と建設地との距離を計算して、定められた距離減衰式に従い、建設地の地盤の基盤における前記各震源断層で地震が生じた場合の建物周期ごとの最大加速度応答を示す加速度応答スペクトルを作成し、
前記地盤増幅率計算手段は、
建設地の地盤につき、入力された多層の地盤の情報から、表層地盤と基盤との2層に大別してモデル化し、表層地盤は一様な地質の1層または複数層の層構成とし、基盤の地震動に対して表層地盤の地表面に生じる地震動が増幅される率である表層地盤固有の地震動増幅率を定められた理論計算式に従って求め、
前記応答解析手段は、
建物を1質点系にモデル化してそのモデルの荷重と変形の関係を示す耐力スペクトルを求め、前記基盤の加速度応答スペクトルと前記表層地盤固有の地震動増幅率とから、前記各震源断層で生じた地震による建設地の地表面における加速応答スペクトルと変位応答スペクトルより設定される地表面応答スペクトルを求め、これら耐力スペクトルと地表面応答スペクトルとの関係から、応答解析結果として前記各地震による建物の応答変位,応答加速度を求め、
前記損失・PML計算手段は、
前記応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算し、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を分布関数でモデル化して90%非超過損失を算定し、90%非超過損失の大きくなる地震より年発生確率の累積である累積確率を求め、累積確率が定められた値となる90%非超過損失であるPML値を計算し、
前記出力部は前記損失・PML計算手段の計算結果となるPML値を表示装置の画面に出力する、
ことを特徴とする建物の地震被害損失の評価装置。 - 請求項1において、前記地盤増幅率計算手段は、表層地盤を全体が一様な1層としてモデル化する建物の地震被害損失の評価装置。
- 請求項2または請求項3において、前記地盤増幅率計算手段は、入力された多層の地盤を、粘性土と砂質土とに分類し、粘性土と砂質土の割合で、1層にモデル化された表層地盤の歪依存特性を定める建物の地震被害損失の評価装置。
- 請求項1ないし請求項3のいずれか1項において、前記応答解析手段は、前記モデル化による建物のモデル化を、建物の各層の荷重と変形の関係から1質点系に集約したモデルとする建物の地震被害損失の評価装置。
- 入力過程と、地震抽出過程と、地盤増幅率計算過程と、応答解析過程と、損失・PML計算過程と、出力過程とを1台の情報処理機器を用いて実行し、
前記入力過程は、評価対象となる建物の建設地の位置および地盤の情報、並びに建物の構造計算結果および部位毎の価格の情報の入力を受付け、
前記地震抽出過程は、
建設地を基準とする対象範囲内に存在する各震源断層、並びにその各震源断層で想定される地震のマグニチュード、年発生確率、震源深さ、および位置の情報をデータベースより抽出し、抽出した位置の情報から各震源断層と建設地との距離を計算して、定められた距離減衰式に従い、建設地の地盤の基盤における前記各震源断層で地震が生じた場合の建物周期ごとの最大加速度応答を示す加速度応答スペクトルを作成し、
前記地盤増幅率計算過程は、
建設地の地盤につき、入力された多層の地盤の情報から、表層地盤と基盤との2層に大別してモデル化し、表層地盤は一様な地質の1層または複数層の層構成とし、基盤の地震動に対して表層地盤の地表面に生じる地震動が増幅される率である表層地盤固有の地震動増幅率を定められた理論計算式に従って求め、
前記応答解析過程は、
建物を1質点系にモデル化してそのモデルの荷重と変形の関係を示す耐力スペクトルを求め、前記基盤の加速度応答スペクトルと前記表層地盤固有の地震動増幅率とから、前記各震源断層で生じた地震による建設地の地表面における加速応答スペクトルと変位応答スペクトルより設定される地表面応答スペクトルを求め、これら耐力スペクトルと地表面応答スペクトルとの関係から、応答解析結果として前記各地震による建物の応答変位,応答加速度を求め、
前記損失・PML計算過程は、
前記応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算し、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を分布関数でモデル化して90%非超過損失を算定し、90%非超過損失の大きくなる地震より年発生確率の累積である累積確率を求め、累積確率が定められた値となる90%非超過損失であるPML値を計算し、
前記出力過程は前記損失・PML計算過程の計算結果となるPML値を表示装置の画面に出力する、
ことを特徴とする建物の地震被害損失の評価方法。 - コンピュータに実行されるプログラムであって、
入力手順と、地震抽出手順と、地盤増幅率計算手順と、応答解析手順と、損失・PML計算手順とを含み、
前記入力手順は、評価対象となる建物の建設地の位置および地盤の情報、並びに建物の構造計算結果および部位毎の価格の情報の入力を受付け、
前記地震抽出手順は、
建設地を基準とする対象範囲内に存在する各震源断層、並びにその各震源断層で想定される地震のマグニチュード、年発生確率、震源深さ、および位置の情報をデータベースより抽出し、抽出した位置の情報から各震源断層と建設地との距離を計算して、定められた距離減衰式に従い、建設地の地盤の基盤における前記各震源断層で地震が生じた場合の建物周期ごとの最大加速度応答を示す加速度応答スペクトルを作成し、
前記地盤増幅率計算手順は、
建設地の地盤につき、入力された多層の地盤の情報から、表層地盤と基盤との2層に大別してモデル化し、表層地盤は一様な地質の1層または複数層の層構成とし、基盤の地震動に対して表層地盤の地表面に生じる地震動が増幅される率である表層地盤固有の地震動増幅率を定められた理論計算式に従って求め、
前記応答解析手順は、
建物につき1質点系にモデル化してそのモデルの荷重と変形の関係を示す耐力スペクトルを求め、前記基盤の加速度応答スペクトルと前記表層地盤固有の地震動増幅率とから、前記各震源断層で生じた地震による建設地の地表面における加速応答スペクトルと変位応答スペクトルより設定される地表面応答スペクトルを求め、これら耐力スペクトルと地表面応答スペクトルとの関係から、応答解析結果として前記各地震による建物の応答変位,応答加速度を求め、
前記損失・PML計算手順は、
前記応答解析結果を用いて、前記各地震が生じた場合の建物の構造部材、非構造部材の損傷確率を計算し、地震毎に損失期待値と部位毎の損失の割合を計算し、地震に対する損失を分布関数でモデル化して90%非超過損失を算定し、90%非超過損失の大きな地震より年発生確率の累積である累積確率を求め、累積確率が定められた値となる90%非超過損失であるPML値を計算し、
前記出力手順は前記損失・PML計算手順の計算結果となるPML値を表示装置の画面に出力する、
ことを特徴とする建物の地震被害損失の評価プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012014046A JP5840513B2 (ja) | 2012-01-26 | 2012-01-26 | 建物の地震被害損失の評価装置・評価方法・評価プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012014046A JP5840513B2 (ja) | 2012-01-26 | 2012-01-26 | 建物の地震被害損失の評価装置・評価方法・評価プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013152197A true JP2013152197A (ja) | 2013-08-08 |
JP5840513B2 JP5840513B2 (ja) | 2016-01-06 |
Family
ID=49048641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012014046A Active JP5840513B2 (ja) | 2012-01-26 | 2012-01-26 | 建物の地震被害損失の評価装置・評価方法・評価プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5840513B2 (ja) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106354897A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种基于gpu的卷积最佳匹配层边界条件实现方法 |
JP2017048583A (ja) * | 2015-09-01 | 2017-03-09 | 東日本旅客鉄道株式会社 | 地盤振動予測方法 |
JP2017531799A (ja) * | 2014-10-25 | 2017-10-26 | ビルディング マグニチュード テクノロジー アナリシス アンド リサーチ サービス カンパニー リミテッドBuilding Magnitude Technology Analysis And Research Service Co.,Ltd. | 既存構造物の耐震規模の算定方法及びそのシステム |
JP2019027984A (ja) * | 2017-08-02 | 2019-02-21 | ビイック株式会社 | 地震による建物変位量の判定方法 |
JP2020125911A (ja) * | 2019-02-01 | 2020-08-20 | 株式会社益田建設 | 家屋の耐震性評価方法 |
CN112613189A (zh) * | 2020-12-31 | 2021-04-06 | 中国再保险(集团)股份有限公司 | 一种考虑地震动幅值不确定性的建筑物易损性曲线建立方法及装置 |
JP2022148889A (ja) * | 2021-03-24 | 2022-10-06 | パナソニックホームズ株式会社 | 建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラム |
CN115963557A (zh) * | 2022-11-25 | 2023-04-14 | 哈尔滨工业大学(深圳) | 基于地震危险一致性的地震动记录选取方法及装置 |
CN115983062A (zh) * | 2022-12-08 | 2023-04-18 | 武汉大学 | 一种基于有限元模型修正的高拱坝地震损伤评估方法及系统 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102032304B1 (ko) * | 2018-05-17 | 2019-10-14 | 부산대학교 산학협력단 | 지진세기에 따른 복수 기기의 동시손상확률 산출 장치 및 방법 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003147970A (ja) * | 2001-11-12 | 2003-05-21 | Takenaka Komuten Co Ltd | 建物の耐震性能評価方法及びその耐震性能評価値に基く改修費用評価方法 |
JP2003296396A (ja) * | 2002-03-29 | 2003-10-17 | Yuji Takahashi | 建物の期待ライフサイクルコスト評価システムおよび期待ライフサイクルコスト評価プログラムを記録した記録媒体 |
JP2005208832A (ja) * | 2004-01-21 | 2005-08-04 | Ohbayashi Corp | 施工中建物の地震損失評価システム、及びそのシステムを利用した施工中建物の地震リスク評価システム、これらのシステムを実行させるためのプログラム、及び、これらのプログラムを記録したコンピュータ読取可能な記録媒体 |
JP2006078272A (ja) * | 2004-09-08 | 2006-03-23 | Sekisui Chem Co Ltd | 地震危険度評価システムおよび建築物選定方法 |
JP2006133871A (ja) * | 2004-11-02 | 2006-05-25 | Shimizu Corp | 建築物評価システムおよび建築物評価方法 |
JP2006318399A (ja) * | 2005-05-16 | 2006-11-24 | Shimizu Corp | 建造物の構造選択支援システムおよび構造選択支援方法 |
JP2007093619A (ja) * | 2006-12-04 | 2007-04-12 | Takenaka Komuten Co Ltd | 地震被害予測装置、地震被害予測方法及び地震被害予測プログラム |
JP2009271684A (ja) * | 2008-05-03 | 2009-11-19 | Kajima Corp | 構造物の定量的耐震性能評価プログラム |
JP2011027481A (ja) * | 2009-07-22 | 2011-02-10 | Ohbayashi Corp | 地震リスク評価方法 |
-
2012
- 2012-01-26 JP JP2012014046A patent/JP5840513B2/ja active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003147970A (ja) * | 2001-11-12 | 2003-05-21 | Takenaka Komuten Co Ltd | 建物の耐震性能評価方法及びその耐震性能評価値に基く改修費用評価方法 |
JP2003296396A (ja) * | 2002-03-29 | 2003-10-17 | Yuji Takahashi | 建物の期待ライフサイクルコスト評価システムおよび期待ライフサイクルコスト評価プログラムを記録した記録媒体 |
JP2005208832A (ja) * | 2004-01-21 | 2005-08-04 | Ohbayashi Corp | 施工中建物の地震損失評価システム、及びそのシステムを利用した施工中建物の地震リスク評価システム、これらのシステムを実行させるためのプログラム、及び、これらのプログラムを記録したコンピュータ読取可能な記録媒体 |
JP2006078272A (ja) * | 2004-09-08 | 2006-03-23 | Sekisui Chem Co Ltd | 地震危険度評価システムおよび建築物選定方法 |
JP2006133871A (ja) * | 2004-11-02 | 2006-05-25 | Shimizu Corp | 建築物評価システムおよび建築物評価方法 |
JP2006318399A (ja) * | 2005-05-16 | 2006-11-24 | Shimizu Corp | 建造物の構造選択支援システムおよび構造選択支援方法 |
JP2007093619A (ja) * | 2006-12-04 | 2007-04-12 | Takenaka Komuten Co Ltd | 地震被害予測装置、地震被害予測方法及び地震被害予測プログラム |
JP2009271684A (ja) * | 2008-05-03 | 2009-11-19 | Kajima Corp | 構造物の定量的耐震性能評価プログラム |
JP2011027481A (ja) * | 2009-07-22 | 2011-02-10 | Ohbayashi Corp | 地震リスク評価方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017531799A (ja) * | 2014-10-25 | 2017-10-26 | ビルディング マグニチュード テクノロジー アナリシス アンド リサーチ サービス カンパニー リミテッドBuilding Magnitude Technology Analysis And Research Service Co.,Ltd. | 既存構造物の耐震規模の算定方法及びそのシステム |
CN106354897A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种基于gpu的卷积最佳匹配层边界条件实现方法 |
JP2017048583A (ja) * | 2015-09-01 | 2017-03-09 | 東日本旅客鉄道株式会社 | 地盤振動予測方法 |
JP2019027984A (ja) * | 2017-08-02 | 2019-02-21 | ビイック株式会社 | 地震による建物変位量の判定方法 |
JP2020125911A (ja) * | 2019-02-01 | 2020-08-20 | 株式会社益田建設 | 家屋の耐震性評価方法 |
JP7253231B2 (ja) | 2019-02-01 | 2023-04-06 | 株式会社益田建設 | 家屋の耐震性評価システム |
CN112613189B (zh) * | 2020-12-31 | 2023-12-26 | 中国再保险(集团)股份有限公司 | 一种考虑地震动幅值不确定性的建筑物易损性曲线建立方法及装置 |
CN112613189A (zh) * | 2020-12-31 | 2021-04-06 | 中国再保险(集团)股份有限公司 | 一种考虑地震动幅值不确定性的建筑物易损性曲线建立方法及装置 |
JP2022148889A (ja) * | 2021-03-24 | 2022-10-06 | パナソニックホームズ株式会社 | 建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラム |
JP7444811B2 (ja) | 2021-03-24 | 2024-03-06 | パナソニックホームズ株式会社 | 建物の被災度推定方法、建物の被災度推定システム、建物の被災度学習方法、建物の被災度学習システム及びプログラム |
CN115963557A (zh) * | 2022-11-25 | 2023-04-14 | 哈尔滨工业大学(深圳) | 基于地震危险一致性的地震动记录选取方法及装置 |
CN115963557B (zh) * | 2022-11-25 | 2024-03-12 | 哈尔滨工业大学(深圳) | 基于地震危险一致性的地震动记录选取方法及装置 |
CN115983062A (zh) * | 2022-12-08 | 2023-04-18 | 武汉大学 | 一种基于有限元模型修正的高拱坝地震损伤评估方法及系统 |
CN115983062B (zh) * | 2022-12-08 | 2023-09-12 | 武汉大学 | 一种基于有限元模型修正的高拱坝地震损伤评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
JP5840513B2 (ja) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5840513B2 (ja) | 建物の地震被害損失の評価装置・評価方法・評価プログラム | |
Lu et al. | A numerical coupling scheme for nonlinear time history analysis of buildings on a regional scale considering site‐city interaction effects | |
van Elk et al. | A probabilistic model to evaluate options for mitigating induced seismic risk | |
Ancheta et al. | NGA-West2 database | |
Rodgers et al. | Broadband (0–4 Hz) ground motions for a magnitude 7.0 Hayward fault earthquake with three‐dimensional structure and topography | |
Kenda et al. | Subsurface structure at the InSight landing site from compliance measurements by seismic and meteorological experiments | |
Karapetrou et al. | “Time-building specific” seismic vulnerability assessment of a hospital RC building using field monitoring data | |
US9383465B2 (en) | Quantitative analysis of time-lapse seismic data | |
KR101663679B1 (ko) | 기하학적 증폭계수를 고려한 지진시 산사태 재해도 도출방법 및 도출시스템 | |
Huang et al. | Fragility assessment of tunnels in soft soils using artificial neural networks | |
Trifunac et al. | Evolution of accelerographs, data processing, strong motion arrays and amplitude and spatial resolution in recording strong earthquake motion | |
Erdogan et al. | Investigation of the seismic behavior of a historical masonry minaret considering the interaction with surrounding structures | |
JP7373616B2 (ja) | 構造物の常時微動に基づく構造物の診断評価方法 | |
Gribovszki et al. | Estimating the upper limit of prehistoric peak ground acceleration using an in situ, intact and vulnerable stalagmite from Plavecká priepast cave (Detrekői-zsomboly), Little Carpathians, Slovakia—first results | |
JP2012168008A (ja) | 地震被害判定システム、地震被害判定システムを備えた構造物、及び地震被害判定プログラム | |
Sabia et al. | Model updating to forecast the dynamic behavior of the Ghirlandina Tower in Modena, Italy | |
Boyd et al. | Crustal deformation in the New Madrid seismic zone and the role of postseismic processes | |
Choudhury et al. | A review of seismic hazard assessment of Gujarat: a highly active intra-plate region | |
Boukri et al. | Seismic risk and damage prediction: case of the buildings in Constantine city (Algeria) | |
Calcina et al. | Ambient vibration tests of an arch dam with different reservoir water levels: experimental results and comparison with finite element modelling | |
JP4506625B2 (ja) | リアルタイム地震情報を利用した地震動の予測システム | |
Johari et al. | Stochastic nonlinear ground response analysis: A case study site in Shiraz, Iran | |
Gu et al. | A record-based method for the generation of tridirectional uniform hazard-response spectra and ground motions using the Hilbert-Huang transform | |
Dragos et al. | Characteristics and results of the near real-time system for estimating the seismic damage in Romania | |
Chen et al. | Numerical damage localisation for building systems including dynamic soil-structure interaction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20141226 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20151015 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20151027 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20151111 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5840513 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |