6. 圧密理論 6. 圧密理論 6.. 圧密方程式の誘導 粘土層の圧密原因とメカニズム 地下水位の低下 盛土建設 最終圧縮量と圧縮速度
6. 圧密理論 記号の統一間隙水圧 ( 絶対圧 ): u 間隙水圧 (gauge 圧 ): u u p a ( 大気圧 ) 過剰間隙水圧 : Δu ( 教科書は これを u と記している 初期状態が u p a で u の時で uδu の状態を対象にしている ) 微小の増分 : du, d( u), dt の様に表現 水頭 h の定義両者が異なるから 間隙水が移動する 地表 静水時の全水頭 パイプ内水位 砂層 ( 静水時のパイプ内の 地下水位 ΔhΔu/γ からの水面の高さ ) 厚さ の粘土層 - h : 位置の水頭 ( 静水圧水頭 ) ( 全水頭 ) u/γ h: 全水頭 -: 静水圧水頭 ΔhΔu/γ ( 位置水頭 ) ( 静水時の の位置からの h の分布 パイプ内の水面の高さ ) 岩盤 head この位置を ead (h) と定義する 静水時の h の分布 全水頭 : パイプ内の の位置からの水面の高さ 間隙水圧 u の定義 地表 パイプ内水位 砂層 地下水位 Δu/γ 厚さ の粘土層 - h ( 静水圧 : u /γ ) u/γ u: 間隙水圧 u : 静水圧 Δu: 過剰間隙水圧 ( 位置水頭 ) (- ) γ 岩盤 γ 間隙水圧 u この位置を ead (h) と定義する
6. 圧密理論 3 ( 全 ) 水頭 h: パイプ内の水面の の面からの高さ 全水頭 静水圧水頭 + 位置の水頭 +( 過剰間隙水圧 Δu)/γ h ( ) + + Δu/γ (6.3) ( 静水圧時の全水頭 )+ Δu/γ 全水頭 - 位置の水頭 静水圧水頭 (u /γ )+( 過剰間隙水圧 Δu)/γ ( γ ) 間隙水圧 u 静水圧 u + 過剰間隙水圧 Δu pore ater pressure hydro-static pressure + Excess pore ater pressure 過剰間隙水圧 Δu 間隙水圧 u 静水圧 u Excess pore ater pressure pore ater pressure hydro-static pressure 有効応力 全応力 - 間隙水圧 σ σ u σ γ total ( ) d( depth) ; 浸透圧が無い静水時 σ γ '( ) d( depth) 水の流れを支配する Darcy 則 i k ( 今 i の正は上方への流れと定義する ) での動水勾配 : i h ( + u / γ ) (6.3) 式で は定数だから ( u / γ ) ( u) γ ( マイナス記号に注意 ) 平均流速 i k ( u) k γ (i 正は 上方への流れを意味する ) (6.4) 要素 a ( 飽和した土の要素 ) この場合 ( u) は負 ( u) は正 a 流出速度 + d d 水平面断面積 A 流入速度 3
6. 圧密理論 4 要素 a から 時間増分 dt の間に排出される水の量 : dq d A dt (6.) (6.) + (6.4) から ; ( u) dq { k } d A dt (6.5) γ Darcy 式 : これまでに用いた 唯一の物性 (k) を含む式 連続の式要素 a は 飽和しているから 流出した間隙水の量 dq 要素の体積減少量 dv 体積ひずみ 間隙の減少量/ 体積 : 時間 dt の間に生じた間隙比の減少量は e dq dv dε d A d A dt + e dε d A dv de Vs de dε V ( + e) V + e e de dt なので t (6.6) s (6.5) (6.6) から γ k ( u) e { } + e γ k ( u) { } e + e (6.7a) (6.7) 次に u と e の関連をつける必要がある 粘土の圧縮性を表す式 ( 有効応力 ~ひずみ関係 ) いま ある深さで全応力 σが一定であり 静水圧を u とすると σσ + uσ + Δu + u 一定 また u 一定 σ 従って d(δu) - dσ また dσ ' dε (m は体積圧縮係数 )[ 右図 ] m 従って d( u) d m ε de m + e (6.) ( 両辺 正 ) m が有効応力 σ の変化に対して一定で σ に対しても一定で m m である場合は 時間的にも一定となり となる また (6.) 式の両辺を時間 σ ' e t で偏分すると d( u) dt de dt である 以上から 次式が得られる t t e dt dt t m + e t ; u e ; 従って e t m + e t + e t m ε (6.) t m 4
6. 圧密理論 5 (6.) を (6.7) に代入して 次式を得る γ m k ( u) { } k が深さ 方向に一定であるとすると 次の圧密方程式を得る k m γ ( u) 圧密方程式 u c ( u) : c k γ m (6.) c : 圧密係数 (coefficient of consolidation): 圧密速度を表す係数単位 : cm /sec) この式は m と k が空間的時間的に定数であると言う重大な仮定に基づいて導かれた式である 圧密方程式の意味 : 深さ での過剰間隙水圧 Δu の時間的変化率 ( 実際に負 ) は ; ( u) a) 過剰間隙水圧の深さ方向の変化の変化率に比例 k b) 圧密係数 c γ m に比例 c が大きいほど 過剰間隙水圧の減少率 が大きくなり 圧密が早く終わる t b) c は 透水係数 k に比例 b) c は 体積圧縮係数 m に逆比例 5
6. 圧密理論 6 例 地表面 Δp の瞬間載荷 i γ 動水勾配 i は正で 方向に一定ではない * γ 水の流れ ただし t 大 t ( u) <, < Δp 過剰間隙水圧 Δu 不透水面 * 上方ほど i は大きい 上方ほど流速 は大きい 下方から絞り出された水を集めて 流れているから ( u) < から 圧密方程式により u < である すなわち 時間とともに過剰間隙水圧 Δu は減少する また 下端 では水の流れはない 従って u ここは 直角 例 地表面 i γ 動水勾配 i は正で 方向に一定 *. 水の流れただし ( u), u < で に対して一定 過剰間隙水圧 Δu 透水面 * 流速 は 深さ方向で一定 地盤の間隙比 ( 密度 ) は どの深さでも時間的に一定 圧密現象は 生じていない ( u) から, である 従って 時間に対して過剰間隙水圧 Δu は一定 下端 でも水の流れはある 従って u < ここは 直角ではない これは 下端から常に時間的に一定量の水が湧いてくる定常流 6
6. 圧密理論 7 例 3 地表面 あるいは 急速埋め立て直後 水の流れ 砂地盤の完全液状化 σ ' 定常流 圧密 静水圧 u 有効応力 σ u 過剰間隙水圧 Δu 不透水面この範囲 t は時間とともに減少地表面 tt の時 この範囲では勾配 i は正で 方向に一定 u は t に対して一定であり 水の流れ ( u), u < であり 密度が t 対して一定 t の増加 tt 過剰間隙水圧 Δu 不透水面 t tt の時 この範囲では 動水勾配 i は正で 方向に一定ではなく ( u) <, u < であり 従って u < であり 沈殿が生じていて 密度が経時的に増加している範囲 下端から水は湧きだしておらず 下端から沈殿が開始される 7
6. 圧密理論 8 6.. 圧密方程式の解 例 地表面 ( 排水面 ) i γ 地表面でΔp が瞬間載荷されたとする. 水の流れ t t 大 過剰間隙水圧 Δu 不透水面 Δp c ( u) : c k γ m (6.) 距離に関する二回の偏微分 二つの境界条件が必要 Δu( ) ; ( ) 時間に関する一回の偏微分 初期条件が必要 Δu() (t) を与える ; 一般的には 任意の初期分布 ( 上図の例では Δu() (t) Δp) 変数分離法 ( u c ) k : c (6.) γ m Δu Z() X(t) (6.9) を式 (6.) に代入すると 次式を得る Z X t c Z X 両辺を Z X c で除して X c X t Z β ( ) ( 定数 ) (6.3) Z ここで 左辺は t だけの関数であり 右辺は だけの関数であるので 両辺が等しいと言うことは 両辺は定数であることを意味する 後での便利さを考えて定数を-β と置く この式の解は 次の二つの式である Z( ) + β Z( ) (6.3) Z( ) Acos β + Bsin β (6.33) 8
6. 圧密理論 9 X ( t) + cβ X ( t) X ( t) C exp( β c t) (6.3) X ( t) C exp( β c t) (6.34) C c が大きいほど X が ( 過剰間隙水圧 Δu が ) 早く減少する time t 式 (6.33) に二つの境界条件を適用して 定数 β を求める まず で u ( ) である X(t) なので この境界条件は Z ( ) を意味する Z 一方 式 (6.33) から A βsin β+ B βcos β 従って Z ( ) を用いると Z ( ) + β B 従って B 従って Z A cos β が得られる 次に 境界条件 Δu( ) は X(t) なので Z( ) を意味する これを用いると 次式 が得られる Z( ) A cos(β ) A だと (6.33) から常に Z 従って常に Δu である これは有り得ない A なので Z( ) A cos(β ) が得られる この式から 次式が得られる cos(β ) これの一般解は 従って Z n β n A n n の時は β n π (6.35) n cos( π ) (6.36) Z π, Z A cos π n Z A 9
6. 圧密理論 従って (6.35) を (6.34) に代入すると n n に関する X の解が以下のように求まる n π Xn( t) Cexp( β c t) Cnexp c t n Cn exp π T (6.37) 式は 物性 c 寸法 を含まない一般的な式であることに注意 c t ここで T 正規化された無次元の時間 ; 時間係数である n Xn( t) Cexp( β ct ) Cnexp π T T c t (6.37) C c と によらない関係 時間係数 T 従って (6.9) 式の 特定の n の値にに対する解 は n n un Xn( t) Zn( ) An Cnexp π T cos π 一般解は全ての n に対する解の集合であるので 未定係数 A n C n A n とおくと次式を得る n n u un An exp π T cos π n n (6.39) (6.39) 式に 既知の初期条件 ;Δu() (t) すなわち T の時 Δu() (T ) を適用すると n u ( )( T ) An cos π n (6.4) m 次に 未定係数 A n を求める このために 両辺に cos π を掛けて の間で積分する m u( )( T ) cos π d A n n n m cos( π ) cos( π ) d (6.4) cos( x) + 右辺は n m の項以外は全てゼロであり (cos x) を用いると 次式を得る
6. 圧密理論 m u( )( t ) cos π d A cos{(m ) π } + d Am (6.4) m ここで m を n として A m を A n とすると次式を得る A n n u t ( )( ) cos( π ) d (n,, ) (6.43) したがって (6.9) 式の一般解は n n n [ exp{ ( ) T} cos( ) u ( )( t ) cos( )] u π π n π d 時間に関する項 初期条件が Δu() (t) Δp の時 (6.44) の中で 距離に関する項 (6.44) n u ( )( t ) cos( π ) d p n cos( π ) d n p [ sin( π )] (n ) π n p ( p [ sin( π ) ] ( n ) π (n ) π したがって 一般解は ) n + n+ 4 p n n u T n [( ( ) exp{ ( π ) } cos( π )] (6.45) π n n ここで exp{ ( π ) T} は 過剰間隙水圧 Δu の時間的変化を n cos( π ) は 空間的変化を表す t ; Δu Δp Time 過剰間隙水圧 Δu Δp
6. 圧密理論 6..3 圧密度 圧密完了時の総地盤沈下量 S f に対する 現在までに終了した沈下量 S の比は? 沈下量 S t ; Δu(t) Δp S f 時間増加 Δp 三つの圧密度 ; 過剰間隙水圧 Δu に関する圧密度 U(T). 過剰間隙水圧 Δu より浅い () 各深さ での local な圧密度 U l (l は local の l) time t () 層全体に対する平均圧密度 U(T); U は T の関数 ; T () 過剰間隙水圧に関する各深さ での local な圧密度 U Δσ u の変化量 Δu(t)- 圧密進行中のΔu Δu(t) Δu(t) -Δu l u( t ) u u( t ) ct 圧密進行中のΔu time t () 過剰間隙水圧に関する全層に対する圧密度 u( t ) u u u( ) d U ( T ) u( t o) u( t o) u( )( t ) d
6. 圧密理論 3 (3) 全層に対する沈下に関する圧密度 U S S/S f (S f は最終沈下量 ). time t Us と U(T) の関係を m がσ に依らない ( すなわち深さ に依らない ) と仮定して求める 沈下計算は 次のようになる S ε d m σ ' d ( 深さに対する積分 ) σ σ ' + u は一定であり σ に等しいので それぞれの深度 において次式がなりたつ Δσ u の変化量 Δu(t)- 圧密進行中のΔu Δu(t) Δu(t) ΔuΔσ 圧密進行中の Δu time t 沈下量 S t ; Δu(t)Δp 時間経過 圧密進行中のΔu Δp 従って t t Δσ Δu(t) - 圧密進行中の Δu 過剰間隙水圧 Δu S ε d m σ ' d m { u( )( t ) d u( ) d} 一方 Δσ ( 圧密終了時 ) Δu(t) なので S f (S henδu) m u( )( t ) d ( m u ( )( t ) ) 3
6. 圧密理論 4 従って, 両式から S S m u( ) d が得られる 従って 次式が得られる f S m U s S u( ) d m u( ) d S f m u( )( t ) d f u( ) d U(T) u( )( t ) d 従って U s U(T); すなわち U s は T の関数となる これは m がσ に依らない ( すなわち深さ に依らない ) と言う仮定が厳密には正しくないので 厳密には正しくない しかし 多くの場合は 十分な近似 この仮定を用いないと Us を T の関数として示すことは 非常に難しい ct 次に U s U(T) とT の関係を求める この関係は 初期過剰間隙水圧の深さ方向の分布型 Δu()(t) によって異なる 今 Δu()(t)Δp ( 一定 ) とする すると n+ 4 p ( ) n n u ( ) exp{ ( π ) T} cos( π ) π n n (6.45) から次式が得られる o n+ 4 p ( ) n u( ) d exp{ ( π ) T} π n n n cos( π ) d n+ 4 p ( ) n n ( exp{ ( π ) T} [sin( π )] π n n (n ) π 従って 次式が得られる U s U(T) u( ) d u( ) d u( )( t ) d p 8 n [( exp{ ( π ) T} π ( n ) n ct U s U(T) は T だけの関数であることに注意 (-) n+ (6.55) 4
6. 圧密理論 5 U(T) ( 教科書 4 頁 図 6.8)...4 n の解.6 厳密解 (6.55.) 式.8 U(T).9 U(T)....4.6.8. T ct hen t 近似解 U( T ) π T (6.63) T.5 程度までは 厳密解と殆ど同じ 圧密試験結果の解析に利用 式 (6.55) は複雑すぎて そのまま実務に用いるのは難しい 多くの場合 この図を用いる ( 図解法 ) n だけの解は π u ( )( T ) u cos( ) (6.4-a) π π u ( ) exp( T) u c os( ) (6.45-a) 4 U s U(T) exp( π T ) (6.55-a) 4 沈下量 S t ; Δu(t) (n) Time Δu 過剰間隙水圧 Δu 5
6. 圧密理論 6 図 6.8 の用い方 : C [ 土質試験で求めた値 ] が既知として 粘土地盤の排水距離 ( 既知 ) 所定の時間 t TC t/ を計算する ( 上図 ) U(T) を求める 時間 t の時の沈下量 S S f U(T) を求める [S f は別途計算 ] ct 9 a) U 9 % の時 T. この時殆ど圧密は終了したと見なせる 実際の時間は t 9 γ m T ( この形に注目 ) c c k 圧密終了時間は a) に比例 b) m に比例 c) k に逆比例 実際は m と k を別々に測定しない c の値を報告する ct b) U 99 % の時 T 3. T. 以降は 圧密の進行が極端に遅くなる 圧密試験 : 供試体初期厚さ cm: 両面排水 : cm (t が小さくなるから ) 砂層と sand drain の機能 水平粘土地盤に対しては Sand drain 排水距離が圧倒的に小さくなる 粘土層 ( 厚さ ) (/3) (/3) ここに排水砂層があると圧密終了時間は /9 に短縮 岩盤 軟弱粘性土を用いた盛土では 水平に排水層を入れる 6
6. 圧密理論 7 正規化した圧密方程式 c ( u) : c k γ m (6.) T ct ( ), U*Δu/Δp ( ); Z / ( ~) を用いると ( u / p) ( u / p) ( c t / ) ( / ) U * U * T Z (6.-a) (6.-b) 非常に簡単な式になる 7
6. 圧密理論 8 6..4 圧密試験と整理法 圧密試験の目的 : 次の二つの物性を求める a) 土の圧縮しやすさを示す圧縮指数 Cc(e~logp 関係の勾配 ) と b) 土の圧密速度を示す圧密係数 C(U~T 関係での TC t/ ) 標準圧密試験 (Q) 供試体の厚さは 何故 cm のように薄いのか? (A) 圧密終了時間は 排水距離の二乗に比例する 試験が実施出来る程度に 圧密終了時間を短くするためは 供試体の厚さは cm 程度にする必要がある しかし 非常に硬い地盤材料 ( 硬質粘性土 堆積軟岩等 ) では 測定されるひずみ量が非常に小さくなるので ひずみの測定精度が落ちる (Q) なぜ 4 時間毎に荷重を 倍にするのか? (A)4 時間は人間の生活のリズム また 倍の率で増加していないと 週間程度の範囲で広い圧力範囲の試験が実施できない しかし 圧力が高くなると測定する圧力の幅が大きくなり過ぎて 圧密降伏応力を正確に測定できなくなる 実験を自動化すれば この制限は無くなる 特に 定ひずみ速度試験が良い 圧密係数 C の求め方 具体的には 二つの方法 a) 圧密初期のデ-タを利用する t 法 ある荷重段階 ttでの沈下量 ΔS の測定 最終圧密量 ΔS の推定 ( これは 面倒 ) 圧密度 U(T)ΔS/ΔS を求める U(T)~T 曲線 ( 図 6.8) から U(T) に対応したTを求める 供試体排水距離 T t から C (T )/t を求める b) 圧密終了時までの全結果を用いる lot(t) 法 t 法 ) 沈下量 S と t の図を用意 ) 初期沈下曲線を直線近似する 初期の点は 供試体内に混入した空気の圧縮が原因で 直線から外れることがある このデ-タポイントは無視する 原点の沈下量 S を ds とする ( 理論的根拠 ) T ct 圧密初期の近似解 U ( T ) T π π c t (6.63) 従って (S-dS) S U と t 関係は 直線になるはず 8
6. 圧密理論 9 3) 角度 α を求める ds. α B /.9 ( S ) 厳密解 α.. +.55 t (t 9 ) 近似解 : c U T T ( ) 沈下量 S π π t t t π 理論的には tanα S ds S U( T ) S c f f.886 (a) S f c U ( T ) π c t ここで まだ最終沈下量 S f は求まっていないなので この式を用いて C を求めることは出来ない 3) tanα (.55) tanα の直線を引く この直線と デ-タにフィットさせた曲線との交点 B を求める この B 点は U 9 % の状態を表す この時 厳密解 ( 図 6.8) から T.848.55 の理論的根拠 : T ct t T tanα ( ) t t ( ) S U S c U 9 t t 9. 848 S c. 9 3..55.886.55 tanα S c S c (b) 式 (a) と式 (b) から tanα (.55) tanα ct9 4) B 点の時間座標 (t 9 ) を読みとる この t 9 を Tt. 848 から求まる c 9 入して C を求める 5) S f は 図のようして求めることが出来るが ここでは使わない. 848 t 9 に代 9