土圧

**土圧** (どあつ、)とは、地盤内における土による圧力のことで、状態によって水平土圧の値が変化する。擁壁に裏込めされた土により、擁壁には土圧が作用する。擁壁が転倒しないように設計を行うためには、土圧の算定が重要となる。

## 歴史
### クーロンの土圧理論
1763年パリ条約でイギリスと講和をむすんだフランスは、その際、イギリスに占領されていたマルチニク島が返還され、その破壊された城と町とを再建しなければならなくなった。1764年に再建案がまとまり、一団の軍技術者達がブレスト港からマルチニク島に向けて出港した。その時１人の技術者が病気で倒れ、その代わりとしてクーロンが指名され、急遽赴任させられることになった。

このマルチニク島における（後のFort Desaix）の築城工事は大工事であった。その上、疫病の流行により技術者達が次々と倒れたため、当時27才のクーロンが実質的な責任者として城を完成させなければならなかった。彼は後にこの工事について、以下のように述べている。

1772年にマルチニク島から本国へ帰ったクーロンは、で軍の勤務についた。この翌年の1773年、マルチニク島での研究をもとに、フランス本国へ帰ってきて最初に発表した論文が、有名な土圧論の論文、『建築学に関して、静力学問題における最大・最小原理の応用について（Essai Sur une application des règles de Maximis et Minimis à quelques Problèmes de Statique, relatifs à l'Architecture）』である。この論文は、最大・最小原理という統一的な観点から、さまざまな力学的問題を解いたものであり、「土圧論」もその一つの応用例であった。この最大・最小の原理については、ダニエル・ベルヌーイやレオンハルト・オイラーなどの後に変分法に発展する業績に影響を受けたものと思われる。後にが認めていたことであるが、クーロンの論文は豊富な内容が非常に簡潔な形で書かれていたため、そのほとんどが専門分野で見落とされていた。

現在、我々が一般にクーロンの土圧公式と呼んでいるのは、が1906年に記したものである。
### ランキンの土圧理論
いわゆる「ランキンの土圧論」が書かれたのは、ランキンがグラスゴー大学に着任したその翌年で、さらにその翌年、1857年、王立学士院の紀要（Philosophical Transactions）に、『緩い土の安定について（On the Stability of Loose Earth）』という題で掲載された。このときランキンは37歳であり、これはクーロンが土圧論を発表したのと同じ年齢であった。

この論文の内容は、彼が大学に勤めるようになり、講義の準備の際に土木工事や石工工事の安定について数学的な理論を再検討しているときに考え出したと述べている。ランキンは父が顧問をしていたカレドニアン鉄道会社が進めていたいくつかの鉄道の建設に携わっていたことがあったが、このときに擁壁に作用する土圧に関わりを持ったと考えられる。

## 土圧と水圧の違い
静止した状態にある水において、ある点における水圧はどの方向からも等しい大きさであり、水の単位体積重量にその点よりも上にある水の高さを乗ずることで得られる。土圧の場合、鉛直方向に関しては土の単位体積重量に深さを乗じた値が土圧となり、これは水と同様である。一方、水平方向は、土の単位体積重量に深さを乗じた値の0.4～0.7倍が土圧となる。土の状態によって水平土圧の値が変わり、それぞれ主働土圧、受働土圧、静止土圧と呼ばれる。

## 土圧の種類

土圧の種類は3つあり、以下の通りである。
;主働土圧
:鉛直応力が卓越して土が破壊する時の水平土圧
;受働土圧
:水平応力が卓越して土が破壊する時の水平土圧
;静止土圧
:地盤内で静止している時の水平土圧
水平土圧と変位の関係は右図の通り。
また、地盤の状態はそれぞれ下図の通りである。それぞれK_0\,\!、K_a\,\!、K_p\,\!は静止土圧係数、主働土圧係数、受働土圧係数である。

主働土圧と受働土圧の計算方法は2つ存在し、それぞれランキン土圧、クーロン土圧と呼ばれる。かつてはテルツァギの土圧図表による土圧計算も多く用いられていたが、現在はほぼ用いられていない。

## ランキン土圧
ランキン土圧を算出する時は下記のような仮定を用いている。
#擁壁は考えない。(擁壁の摩擦及び形状は考えない。) 
#塑性平衡状態となる。モール・クーロンの破壊規準に従う。 
#傾斜角を考慮しない。

;主働土圧状態
:モール・クーロンの破壊規準の主応力表示は下記の通りである。
::\sigma_1-\sigma_3=2\mathrm{ccos}\phi+(\sigma_1+\sigma_3)\mathrm{sin}\phi\,\!
:主働土圧の時、最大主応力は\sigma_v\,\!、で最小主応力は\sigma_{ha}\,\!、なので、上式に代入すると、以下の式を得る。
::\sigma_v-\sigma_{ha}=2\mathrm{ccos}\phi+(\sigma_v+\sigma_{ha})\mathrm{sin}\phi\,\!
:上式を\sigma_{ha}\,\!について整理する。
::\sigma_{ha}=\sigma_v\mathrm{tan}^2\left( \frac{\pi}{4}- \frac{\phi}{2}\right)-2\mathrm{ctan}\left( \frac{\pi}{4}- \frac{\phi}{2}\right)\,\!
:また、主働土圧係数K_a\,\!を用いて上式を書く。
::\sigma_{ha}=\sigma_vK_a-2\mathrm{c}\sqrt{K_a}\,\!・・・・①
::K_a=\mathrm{tan}^2\left( \frac{\pi}{4}- \frac{\phi}{2}\right)\,\!
;受働土圧状態
:
:モール・クーロンの破壊規準の主応力表示は下記の通りである。
::\sigma_1-\sigma_3=2\mathrm{ccos}\phi+(\sigma_1+\sigma_3)\mathrm{sin}\phi\,\!
:受働土圧の時、最大主応力は\sigma_{hp}\,\!、で最小主応力は\sigma_v\,\!、なので、上式に代入すると、以下の式を得る。
::\sigma_{hp}-\sigma_v=2\mathrm{ccos}\phi+(\sigma_{hp}+\sigma_v)\mathrm{sin}\phi\,\!
:上式を\sigma_{hp}\,\!について整理する。
::\sigma_{hp}=\sigma_v\mathrm{tan}^2\left( \frac{\pi}{4}+\frac{\phi}{2}\right)-2\mathrm{ctan}\left( \frac{\pi}{4}+\frac{\phi}{2}\right)\,\!
:また、受働土圧係数K_p\,\!を用いて上式を書く。
::\sigma_{hp}=\sigma_vK_p-2\mathrm{c}\sqrt{K_p}\,\!・・・・・②
::K_p=\mathrm{tan}^2\left( \frac{\pi}{4}+\frac{\phi}{2}\right)\,\!

## クーロン土圧
クーロン土圧を算出する時は下記のような仮定を用いている。
#粘着力の無い砂質土を対象とする 
#壁体の背後の土の中に直線状の滑り面が生じ、くさび状の土塊が滑り面に沿って動く
クーロン土圧はランキン土圧よりも適用範囲が広く、壁体との摩擦、壁体の傾斜、背後の地表面の傾斜も考慮している。 
;主働土圧状態
:クーロンの主働土圧計算時の地盤の状態は下図の通りである。
:
: 
:右図を連力図という。土のくさびの重量W\,\!は既知。３つ力のベクトルが閉じた 三角形になるように主働土圧の合力P_a\,\!と滑り面に作用する力F\,\!の大きさを決める。その時の\theta\,\!が滑り面の角度となる。
:(土と壁体の摩擦角\delta\,\!と\phi\,\!内部摩擦角は既知、\alpha\,\!と\beta\,\!は土や擁壁の形を決めるものなので既知) 
:連力図に着目すると、正弦定理より以下の関係式を得る。
::\frac{W}{\mathrm{sin}(\alpha+\delta-\theta+\phi)}=\frac{P_a}{\mathrm{sin}(\theta-\phi)}\!
:したがって、P_a\,\!は以下の通り。
::P_a=\frac{W\mathrm{sin}(\theta-\phi)}{\mathrm{sin}(\alpha+\delta-\theta+\phi)}\!
::W\,\!は土の重さなので別途計算する必要がある。
:つまり、土塊の体積V\,\!を計算し、土の単位体積重量を乗ずればよい。
:この時、土の体積は単位奥行きあたりの体積であるので、実質的には土くさびの面積S\,\!を求めればよい。
:右図より土くさびの面積は以下の通りとなる。
:
::S=\frac{1}{2}AB*AC\mathrm{sin}(\alpha-\theta)
:また、正弦定理より以下の関係式を得る。
::\frac{AB}{\mathrm{sin}(\theta-\beta)}=\frac{AC}{\mathrm{sin}(\pi-\alpha+\beta)}=\frac{AC}{\mathrm{sin}(-\alpha+\beta)}
:辺ABを壁体の高さHを用いて表す。
::AB=\frac{H}{\mathrm{sin}\alpha}
:以上より土くさびの面積は以下の通りである。
::S=V=\frac{H^2\mathrm{sin}(\alpha-\theta)\mathrm{sin}(-\alpha+\beta)}{2\mathrm{sin}^2\alpha \mathrm{sin}(\theta-\beta)}
:したがって、土の重さは以下のように得られる。
::W=\frac{\gamma_tH^2\mathrm{sin}(\alpha-\theta)\mathrm{sin}(-\alpha+\beta)}{2\mathrm{sin}^2\alpha \mathrm{sin}(\theta-\beta)}
:つまり、クーロンの主働土圧は以下のように導かれる。
::P_a=\frac{\gamma_tH^2\mathrm{sin}(\alpha-\theta)\mathrm{sin}(-\alpha+\beta)\mathrm{sin}(\theta-\phi)}{2\mathrm{sin}^2\alpha\mathrm{sin}\theta-\beta)\mathrm{sin}(\alpha+\delta-\theta+\phi)}
::\theta\,\!の最適解は極限定理の上界定理の考え方からP_a\,\!を最大とする時の\theta\,\!となる。
:すなわち、以下の関係を満たす時の\theta\,\!が最適解となり、その時のP_a\,\!がクーロンの主働土圧となる。
::\frac{\partial P_a}{\partial \theta}=0
:クーロンの主働土圧は以下の通りである。
::P_a=\frac{1}{2}\gamma_tK_aH^2\,\!
::K_a=\frac{\mathrm{sin}^2(\phi-\alpha)}{\mathrm{sin}^2\alpha \mathrm{sin}(\alpha+\delta)\left[1+\sqrt{\frac{\mathrm{sin}(\phi-\beta)\mathrm{sin}(\phi+\delta)}{\mathrm{sin}(\alpha-\beta)\mathrm{sin}(\alpha+\delta)}}\right]^2\,}\,\!
;受働土圧状態
:クーロンの受働土圧計算時の地盤の状態は下図の通りである。
:
:主働土圧と同様に連力図を用いて解く。正弦定理より以下の関係式を得る。
:
::\frac{W}{\mathrm{sin}(\alpha-\delta-\theta-\phi)}=\frac{P_p}{\mathrm{sin}(\theta+\phi)}\!
:したがって、P_p\,\!は以下の通り。
::P_p=\frac{W\mathrm{sin}(\theta+\phi)}{\mathrm{sin}(\alpha-\delta-\theta-\phi)}\!
:また、右図より土くさびの面積は以下の通りとなる。
::S=\frac{1}{2}AB*AC\mathrm{sin}(\alpha-\theta)
:また、正弦定理より以下の関係式を得る。
::\frac{AB}{\mathrm{sin}(\theta-\beta)}=\frac{AC}{\mathrm{sin}(\pi-\alpha+\beta)}=\frac{AC}{\mathrm{sin}(-\alpha+\beta)}
:辺ABを壁体の高さHを用いて表す。
::AB=\frac{H}{\mathrm{sin}\alpha}
:以上より土くさびの面積は以下の通りである。
::S=V=\frac{H^2\mathrm{sin}(\alpha-\theta)\mathrm{sin}(-\alpha+\beta)}{2\mathrm{sin}^2\alpha \mathrm{sin}(\theta-\beta)}
:したがって、土の重さは以下のように得られる。
::W=\frac{\gamma_tH^2\mathrm{sin}(\alpha-\theta)\mathrm{sin}(-\alpha+\beta)}{2\mathrm{sin}^2\alpha \mathrm{sin}(\theta-\beta)}
:つまり、クーロンの受働土圧は以下のように導かれる。
::P_p=\frac{\gamma_tH^2\mathrm{sin}(\alpha-\theta)\mathrm{sin}(-\alpha+\beta)\mathrm{sin}(\theta+\phi)}{2\mathrm{sin}^2\alpha \mathrm{sin}(\theta-\beta)\mathrm{sin}(\alpha-\delta-\theta-\phi)}
::\theta\,\!の最適解は極限定理の上界定理の考え方からP_p\,\!を最大とする時の\theta\,\!となる。
:すなわち、以下の関係を満たす時の\theta\,\!が最適解となり、その時のP_p\,\!がクーロンの受働土圧となる。
::\frac{\partial P_p}{\partial \theta}=0
:クーロンの受働土圧は以下の通りである。
::P_p=\frac{1}{2}\gamma_tK_pH^2\,\!
::K_p=\frac{\mathrm{sin}^2(\phi-\alpha)}{\mathrm{sin}^2\alpha \mathrm{sin}(\alpha-\delta)\left[1-\sqrt{\frac{\mathrm{sin}(\phi+\beta)\mathrm{sin}(\phi+\delta)}{\mathrm{sin}(\alpha-\beta)\mathrm{sin}(\alpha-\delta)}}\right]^2\,}\,\!

## ランキン土圧とクーロン土圧の相違
ランキン土圧における仮定は、以下の通りである。
::\alpha=\frac{\pi}{2},\beta=0,\delta=0\,\!
この時主働土圧および、受働土圧は以下のように得られる。
;主働土圧
::P_a=\frac{\gamma_tH^2\mathrm{cos}(\theta)\mathrm{sin}(\theta-\phi)}{2\mathrm{sin}\theta \mathrm{cos}(\theta-\phi)}\,\!
:先ほどと同様に上界定理を用いる。
::\frac{\partial P_a}{\partial \theta}=0
:この方程式を解くと、\theta\,\!は次のように得られる。
::\theta=\frac{\pi}{4}+\frac{\phi}{2}\,\!
:この時の主働土圧は次のように得られる。
::P_a=\frac{1}{2}\gamma_tK_aH^2\,\!
::K_a=\mathrm{tan}^2\left( \frac{\pi}{4}- \frac{\phi}{2}\right)\,\!
;受働土圧
::P_p=\frac{\gamma_tH^2\mathrm{cos}(\theta)\mathrm{sin}(\theta+\phi)}{2\mathrm{sin}\theta \mathrm{cos}(\theta+\phi)}\,\!
:先ほどと同様に上界定理を用いる。
::\frac{\partial P_p}{\partial \theta}=0
:この方程式を解くと、\theta\,\!は次のように得られる。
::\theta=\frac{\pi}{4}-\frac{\phi}{2}\,\!
:この時の受働土圧は次のように得られる。
::P_p=\frac{1}{2}\gamma_tK_pH^2\,\!
::K_p=\mathrm{tan}^2\left( \frac{\pi}{4}+ \frac{\phi}{2}\right)\,\!

ランキン土圧とクーロン土圧は一見違う理論のように見えるが、実は同じである。ただし、適用できる対象がやや異なるので下に表にしてまとめた。

## 参考文献
*
*
*
*
*
*

## 脚注

## 関連項目
*土質力学
*モール・クーロンの破壊規準
*応用力学