著者

垂水 祐二 伯野 元彦

出版者

東京大学地震研究所

雑誌

東京大学地震研究所彙報 (ISSN:00408972)

巻号頁・発行日

vol.62, no.4, pp.p535-577, 1988-03

被引用文献数

2

---

土とか岩は,砂粒とか土粒子でできていたり,多くの割れ目があったりで,本質的には連続体ではない.勿論,通常の問題では,大変形が生じる訳ではないので,連続体としての解析で何ら支障はないのであるが,土砂崩壊とか液状化とか,大変形時には無理がある.非連続体の解析手法としては,現在までのところ,有限要素法に滑り要素を導入したものや,有限要素法を多少変えた川井のRBSM法(Rigid Body Spring Method),CundallのDEM法(Distinct Element Method)がある.どの手法にも一長一短があるが,有限要素法の流れを汲む前二者では,要素の辺に沿って滑る変形しかできないし,隣り合う要素と要素が離れてしまうなどということもできない.それは,要素の格子点は,隣り合う多くの要素に共有されているので,要素と要素が離れてしまい,今まで1つであった格子点が複数にならなければならないが,それが,全く不可能ではないにせよ,かなり面倒なのである.ところが,実際の現象では,岩石の崩壊の際には,岩塊と岩塊は離れるばかりでなく空中を飛ぶことだってある.一方CundallのDEM法は,要素は元々,バラバラのものとしての取り扱いである.そのバラバラの独立した要素一つ一つについて運動方程式をたてる.したがって,岩石が崩壊する時,岩塊同志が衝突したり,跳んだり,転がったりすることも可能である.現在まで,このDEM法による粒状体シミューショシは,粒体間の間隙に水が飽和している場合については,殆んど行われていなかったので,筆者らは,二次元粒体に水圧を考慮する事によって,地震時の砂の液状化と,大雨の場合の間隙水圧上昇によるQuick Sand現象をこの粒状体シミュレーションで解析することを試みた.そして,或る程度の成功を得た.しかしながら,粒径が砂のように小さくなり,間隙水の体積弾性率の大きい影響で,P波の伝播速度が大きくなり,その結果,計算の安定のため,計算時間刻みを100万分の1秒のように小さくとらねばならず,解析粒数を200個強しかとれなかった.実際の液状化を詳細に模擬するには,3次元解析か必要であり,これは,将来の計算機の高速化をまたねばならない.Various kinds of liquefaction analysis have been carried out in laboratory experiments and using the Finite Element Method (FEM). But, no numerical liquefaction analysis in which sand has been considered a non continuous material has yet been reported. In 1971, the Distinct Element Method (DEM) was introduced by Cundall; a numerical simulation used to analize the behavior of rock, based on the assumption that each individual rock element satisfies the equation of motion. We have developed a modified DEM using Darcy's law that takes into account the pore water pressure. We analyzed the liquefaction of saturated sand under seismic excitation and also the quick sand phenomenon due to water pressure. The assembly model consists of circular elements with log-normal distributed radii, and it is packed by dropping. We used the nonlinear spring in the normal direction between particles in order to express the "Dilatancy" in two dimensional treatment. The excessive pore pressure of the numerical result rose gradually due to the shaking effect. This result agreed with the results of the past laboratory tests. However, the contact of the neighboring particles was not easy to hide even under the high pore pressure, and resultantly, a complete liquefaction phenomenon could not be realized. The complete liquefaction can be analyzed only by the three dimensional analysis or a forced strain input. However, it needs much more computing time, and we have to wait for an advanced and faster computors. We also made the quick sand phenomenon occurr in the model assembly by applying the uplifting water pressure at the bottom of the assembly. A time history of the excessive pore water peressure was also obtained at the point near the bottom or the surface of the assembly.