マンホールの地震時液状化浮き上がり解析
(地震時せん断応力は「略算」で算定) 目次 (1)基本方針 1,本解析の背景 2 2,構造諸元 2 3,本解析の内容 2 4,本解析の目的 2 5,設計方針及び参考文献 2 6.使用プログラム 3 7,変形解析のフロー 3 8,概要図 3 (2)地盤概要 1,地盤の概説 5 (一部省略) 2,ボーリング調査結果 5 3,設計外力 6 (3)自重変形解析の概要 1,ALIDによる自重変形解析法の概説 10 (ALID)(一部省略) 2,解析モデル 11 (4)マンホールの設計 1,マンホールの構造 12 2,構造計算 12 (5)地震時自重変形解析 1,解析モデル 13 (一部省略) 2,ALID解析モデル 15 3,ALID解析結果 16 4,結果の整理 17 マンホールの地震時液状化による浮き上がりを「ALID」により行った。 (株)ブルドジオテクノ(1)基本方針 1,本解析の概説 本解析は、液状化を考慮した静的FEM解析である。解析法は,液状化に伴う堤防及び盛土の変形 を基礎地盤の砂質土層が液状化し,砂質土層の剛性が低下したことに起因して発生するものと仮定し, 解析における外力として土層の自重のみを考慮する静的自重変形解析である。 2,構造諸元 設置場所 ○○県○○市○○町地内 △△線 設備 雨水排水(マンホール)道路地下に設置 地盤 調査報告書参照 3,本解析の内容 本解析は、道路設置マンホールの地震に対する液状化の影響を考慮した浮き上がり量を概略推定する。 FEM解析 液状化自重沈下解析 解析モデル 二次元 解析条件 静的有効応力 ステップ解析 行わない(現状の状態を解析) 液状化の判定 FLによる 液状化定数の計算 N値及び諸定数によるプログラム内部計算 4,本解析の目的 目的 設計状態 想定する荷重など 解析の目的 沈下量予測 地震時 レベル1、2 液状化を考慮した堤防の自重沈下量概略値を推定
6,使用プログラム ・静的FEM解析コード:ALID(ALID研究会) 二次元・有効応力解析 ・マンホールの設計(フォーラムエイト) マンホール構造、地盤液状化の判定 7,変形解析のフロー 土質定数の設定: 試験値の適用 定数の推定・仮定 設計外力の算定 レベル1,2 マンホールの構造 「マンホールの設計」 解析モデル(二次元) 液状化を考慮した土質定数の決定 (1)地震前の応力解析 有限要素法による自重変形解析 (2)液状化時の変形解析 ALIDにより解析 (3)残留変形解析 安全性の確認 8,概要図 ①地盤断面図
(2)地盤概要 1,地盤の概説 調査では1本のボーリング(No.3)が実施され、土質試験が行われている。土層は同一と思われる 層はまとめて、できるだけ簡単なモデルとなるようにしている。ただし、土層の連続については、サン プル及び観察結果等を参考に連続土層を設定した。(ここでは説明のため大きくまとめていることに注 意) 2,ボーリング調査結果 ①No.3:
耐震設計上の地盤種別(道路橋示方書Ⅴ P32による):2種地盤 PAGE=7 液状化の判定資料(道路橋示方書Ⅴ P134による) PAGE=8~9 *ボーリング調査により、地層構成・標準貫入N値が実施され、結果は上図に示す。 *地盤層③~⑦の砂質層は、土の粒度試験が行われている。 3,設計外力 ALIDの解析には液状化に対する抵抗率FLが必要であるため、道路橋示方書Ⅴ(P134)によ り設計値を算定する。 地域区分 (P27) B地域(0.85) レベル1 (P35) khc=0.15x0.85=0.128 レベル2 (P35) khc=0.7(2種地盤)x0.85=0.595 マンホールの設計では、下水道施設の耐震対策指針と解説(P25) 地域区分 (P27) B地域(0.85) レベル1 (P25) 応答変位法による(P120) レベル2 (P28)
地盤種別
ボーリング番号-No.3 平均せん断 弾性波速度 Hi(m) Vsi(m/s) (sec) 1 BS 砂質土 0.95 8 160 0.0059 2 AS 砂質土 2.40 7 153 0.0157 3 AS 砂質土 1.50 13 188 0.0080 4 AS 砂質土 2.33 7 153 0.0152 5 AC 粘性土 1.55 15 247 0.0063 6 AS 砂質土 3.60 9 166 0.0217 7 AS 砂質土 2.21 4 127 0.0174 8 AS 砂質土 2.50 7 153 0.0163 9 AS 砂質土 3.54 5 137 0.0258 20.58 合 計 0.1323Tg = 4 × Σ(Hi/Vsi) =
0.529
(sec) よって、耐震設計上の地盤種別は Ⅱ種地盤 となる。 解説:道路橋示方書Ⅴ耐震設計編P32 耐震設計上の地盤面:十部堅固な地盤(Vsi≧300m/s) 平均N値 Hi/Vsi Ⅰ種地盤:良好な洪積地盤及び岩盤、硬質砂れき層など Tg≦0.2 層 番 号 層名 層種別 層厚 Ⅱ種地盤:Ⅰ種、Ⅲ種にも属さない洪積地盤及び沖積地盤 0.2<Tg≦0.60 Ⅲ種地盤:沖積地盤のうち軟弱地盤 0.60<Tg粘性土層
V
si=100N
i1/31≦N
i≦25
砂質土層
V
si=80N
i1/31≦N
i≦50
1章 液状化
1.1 液状化の判定
(1)レベル1 計算 深度 (m) N値 Ip FC (%) D50 (mm) D10 (mm) R L FL 判定 1.200 2.200 3.200 4.200 5.200 6.200 7.200 8.200 9.200 10.200 11.200 12.200 13.200 14.200 15.200 16.200 17.200 18.200 19.200 8.0 4.0 10.0 15.0 12.0 5.0 8.0 20.0 12.0 10.0 9.0 2.0 5.0 4.0 11.0 4.0 6.0 4.0 5.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 18.9 2.8 2.8 5.1 5.1 5.1 5.1 9.8 9.8 11.1 11.1 6.8 6.8 6.8 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.355 0.256 0.160 0.228 0.204 0.190 0.088 0.136 0.119 0.196 0.116 0.137 0.109 0.120 0.138 0.150 0.159 0.172 0.175 0.177 0.178 0.178 0.176 0.175 0.173 0.171 0.168 0.166 2.570 1.701 1.007 1.326 1.168 1.075 0.494 0.766 0.675 1.120 0.673 0.801 0.650 0.723 × × × × × × × × × × × ○ ○ ○ × ○ ○ ○ ○ 注)地表面から20mまでを表示 判定:×=液状化しない,○=液状化する (2)レベル2 計算 深度 (m) N値 Ip FC (%) D50 (mm) D10 (mm) R L FL 判定 1.200 2.200 3.200 4.200 5.200 6.200 7.200 8.0 4.0 10.0 15.0 12.0 5.0 8.0 0.0 0.0 0.0 18.9 2.8 2.8 0.000 0.000 0.000 0.000 0.000 0.000 0.653 0.388 0.191 0.637 0.694 0.732 1.025 0.558 0.261 × × × × ○ ○ ×計算 深度 (m) N値 Ip FC (%) D50 (mm) D10 (mm) R L FL 判定 12.200 13.200 14.200 15.200 16.200 17.200 18.200 19.200 2.0 5.0 4.0 11.0 4.0 6.0 4.0 5.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 5.1 9.8 9.8 11.1 11.1 6.8 6.8 6.8 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.088 0.152 0.126 0.257 0.123 0.153 0.113 0.128 0.823 0.820 0.814 0.806 0.797 0.788 0.777 0.766 0.107 0.186 0.155 0.319 0.154 0.194 0.145 0.167 ○ ○ ○ ○ ○ ○ ○ ○ 注)地表面から20mまでを表示 判定:×=液状化しない,○=液状化する
(3)自重変形解析法の概説 1,自重変形解析法の概説 ①解析法の概要 本解析法は,二次元断面を対象とし、各節点の変位量を算出する。すなはち、地震によって液状化し た層は剛性が低下し、これによって堤体が沈下する。また、砂質土層の剛性低下に伴う堤防の変形が 地震動が終了した後に発生するものと仮定するため,慣性力に起因する変形は考慮しない。剛性低下 率は液状化に対する抵抗率FL及び細粒分含有率FCの関数として与えられ、地震前の土の剛性、T LとFCが分かれば沈下量を求めることが出来る比較的簡便な方法である。 液状化層の剛性低下の特徴は,下図に示すように、地震前の解析には液状化層のせん断弾性係数Go としてOA線の傾きを使用する。一方、地震後の解析では、液状化層については剛性低下そたせん断弾 性係数G1としてOC線の傾きを使用する。ポアソン比については、地震前の解析では常時ν0を用い、 地震後の解析では体積弾性係数Kが地震前・後で一定となるよう逆算した地震後のポアソン比ν1を用 いる。比液状化層についても、地震前の解析では常時のG・νを使用し、地震後の解析では液状化層の 過剰間隙水圧の影響等によりある程度のせん断剛性の低下を見込んだ物性値G・νを用いて、変形量を 算出する。 (出典:河川堤防の地震時変形量の解析方法) ②解析のフロー 本解析法では, ①液状化層 ②液状化層上部に位置する非液状化層(堤体の盛土と表土層) ③液状化層の下部に位置する非液状化層 の3 種類の土層で液状化時における剛性を設定する。 液状化時における剛性の設定方法は、 ①液状化層は,液状化時の剛性低下を考慮した変位を「応力解放法」により内部計算する。 ②液状化層の上部に位置する非液状化層(堤体の盛土と表土層)については,下部の液状化層の影響を 受け,液状化の程度が大きいほど剛性が低下しやすいため、経験値による剛性の低下率により算定する。 ③液状化層の下部に位置する非液状化層については,液状化の発生に伴う剛性の低下を考慮しない。 液状化層の剛性低下を考慮した有限要素法による変形解析は,下図に示す通り、自重を用いて地震前 のせん断弾性係数とポアソン比を用いた地震前の変形解析と液状化が発生し剛性の低下が生じた際のせ
液状化の判定 解析モデル作成 パラメータ設定 地震前の変形解析 自重変形解析: 非排水条件 液状化時の変形解析 ALID 変形量を合計=自重沈下量 液状化層の 体積圧縮変位量算定 (出典:河川堤防の地震時変形量の解析方法) 2,解析モデル(二次元FEM) ①FEM解析モデル ・側方の境界の影響を少なくするため、左右へ解析地盤深さを目安とし側方地盤域を設定した。 ・モデル(側方地盤域+解析目標域)の左右部分については、モデル全体長(解析目標長の10倍) を考慮した付加地盤(水位考慮)を設定した。 ・モデルの境界条件は、モデルの両側面は付加地盤・水平固定とし、モデル底面は水平変位と 鉛直変位を固定とした。 解析モデル図 各解析結果に記載 ②要素適用モデル
(4)マンホールの設計 1.マンホールの構造
2.構造計算
使用プログラム 「マンホールの設計」(フォーラムエイト) ファイル名:2110
(5)地震時自重変形解析結果 1.解析モデル 解析モデルは上図に示すように、解析目標である「解析目標域」より左右に側方地盤域を設定した。 境界条件は、両側面は付加地盤・水平固定とし、モデル底面は水平変位と鉛直変位を固定とした。 付加地盤の設定長={(解析目標域x10)-(解析目標域+深さx2)}/2 ={(10mx10=100m)-(10+20.6x2=51.2m)}/2 =24.4m 構成要素は図のように①~⑫まで設定した。各要素の設定は以下の通り。 番号 *土質 液状化層 物性モデル 液状化時条件 ① AS MC/DP 剛性低下 ② AS MC/DP 剛性低下 ③ AS MC/DP 剛性低下 ④ AS ○ MC/DP バイリニアー ⑤ AC MC/DP 剛性低下 ⑥ AS ○ MC/DP バイリニアー ⑦ AS ○ MC/DP バイリニアー
各要素のパラメータ 番号 *土質 γ G ν Dr φ c ψ ① AS 17.0 8421 0.33 67.222 32.7 7.6 ② AS 18.0 7368 0.33 53.525 32.0 7.0 ③ AS 17.0 13684 0.33 65.263 34.4 9.4 ④ AS 16.0 7368 0.33 45.670 31.0 6.0 ⑤ AC 18.0 14483 0.33 64.134 3.75 ⑥ AS 19.0 9474 0.33 46.976 31.3 6.3 ⑦ AS 20.0 4211 0.33 29.391 26.8 1.8 ⑧ AS 20.0 7368 0.33 37.008 29.0 4.0 ⑨ AS 20.0 5263 0.33 29.443 26.8 1.8 ⑩ Con 24.5 1041667 0.2 ・γ:単位体積重量(kN/m3)(一般値による推定) ・G:せん断弾性係数(kN/m2)(N値より推定) ・ν:ポアソン比(推定値) ・Dr:相対密度(Meyerhof式より推定) ・φ:せん断抵抗角(゜)(N値より推定) ・c:粘着力(kN/m2)(N値より推定) ・ψ:ダイレイタンシ角(゜)(N値より推定) 変形係数の応力依存性を考慮するためJanbuの実験式による推定式により算定した。 番号 *土質 k n ① AS 95576 0.5 ② AS 38677 0.5 ③ AS 54909 0.5 ⑤ AC 67988 1.0 液状化層の液状化抵抗率FL:レベル2) 番号 *土質 FL ALID ④ AS 0.329 今回のALID解析では、プログラム内で自動計算する ⑥ AS 0.335 ⑦ AS 0.174 ⑧ AS 0.231 ⑨ AS 0.166
2.ALID解析モデル ①ブロック図
3.ALID解析結果 (ここではレベル2タイプⅡを掲載する) ①残留変形図(x1.0)
③液状時間隙水圧分布図
4.結果の整理
①液状化によるマンホール天端浮き上がり量の確認 浮き上がり量(m)