JP3153385B2 - 土の挙動を捉える地盤分析方法 - Google Patents
土の挙動を捉える地盤分析方法Info
- Publication number
- JP3153385B2 JP3153385B2 JP12021693A JP12021693A JP3153385B2 JP 3153385 B2 JP3153385 B2 JP 3153385B2 JP 12021693 A JP12021693 A JP 12021693A JP 12021693 A JP12021693 A JP 12021693A JP 3153385 B2 JP3153385 B2 JP 3153385B2
- Authority
- JP
- Japan
- Prior art keywords
- fluid
- granular
- granular material
- determining
- mesh
- 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.)
- Expired - Fee Related
Links
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、土の挙動を捉える地盤
分析方法に関する。
分析方法に関する。
【０００２】
【従来の技術】土は、固相、液相、及び気相の三相から
構成されており、地盤の液状化、クイックサンド（浸透
水などにより、地盤土砂が浮遊状態になること）、及び
土石流等は、粒状体（固相）と流体（液相）との二相系
によって起こる現象である。
構成されており、地盤の液状化、クイックサンド（浸透
水などにより、地盤土砂が浮遊状態になること）、及び
土石流等は、粒状体（固相）と流体（液相）との二相系
によって起こる現象である。
【０００３】このような２相解析方法としては、個別要
素法（個々の粒状体相互に作用する力を計算し、時間毎
の粒状体の挙動を追う方法）を粒状体と流体の２相系に
適用し、２次元での粒状体の間隙面積の変化から粒状体
に作用する間隙水圧を計算して行うものや、パイピング
（土中に生じた管状の水路から水が吹き出す現象）試験
等をシュミレーションするもの等がある。そして、この
結果より、土の挙動を捉えて、上述したような地盤の液
状化、クイックサンド、及び土石流を分析するのであ
る。
素法（個々の粒状体相互に作用する力を計算し、時間毎
の粒状体の挙動を追う方法）を粒状体と流体の２相系に
適用し、２次元での粒状体の間隙面積の変化から粒状体
に作用する間隙水圧を計算して行うものや、パイピング
（土中に生じた管状の水路から水が吹き出す現象）試験
等をシュミレーションするもの等がある。そして、この
結果より、土の挙動を捉えて、上述したような地盤の液
状化、クイックサンド、及び土石流を分析するのであ
る。
【０００４】これらの解析方法では、隣合う粒状体の間
隙同士の圧力差を平衡させるように流体が流動すると仮
定して、流体と粒状体との相互作用を評価している。
隙同士の圧力差を平衡させるように流体が流動すると仮
定して、流体と粒状体との相互作用を評価している。
【０００５】しかしながら、実際の粒状体の間隙は、２
次元方向だけでなく３次元方向にも間隙が生じているの
で、間隙中を流動する流体の挙動を２次元方向の間隙間
に圧力平衡だけでは、土の挙動を正確に分析できない。
次元方向だけでなく３次元方向にも間隙が生じているの
で、間隙中を流動する流体の挙動を２次元方向の間隙間
に圧力平衡だけでは、土の挙動を正確に分析できない。
【０００６】
【発明が解決しようとする課題】本発明は上記事実を考
慮して、間隙中の流体の挙動を正確に表現し、流体から
粒状体が受ける力を正しく評価して、土の挙動を捉える
地盤分析方法を提供することを目的とする。
慮して、間隙中の流体の挙動を正確に表現し、流体から
粒状体が受ける力を正しく評価して、土の挙動を捉える
地盤分析方法を提供することを目的とする。
【０００７】
【課題を解決するための手段】本発明は、水の流体と土
粒子の粒状体が共に存在する領域に外力が作用したとき
に、流体と粒状体が相互に干渉しながら挙動する状態を
時間毎に捉えることで地盤の液状化、クイックサンド、
又は土石流を分析する地盤分析方法において、各々粒状
体の初期配列と各々粒状体を囲む境界線を初期条件とし
て設定した後、粒状体同士の接触判定を行ない、粒状体
が相互に受ける作用力と、前記境界線で粒状体が受ける
作用力を求める第１ステップと、粒状体に作用する２つ
の作用力によって粒状体が並進或は回転する運動を表す
運動方程式を前進的に数値解析して、時間毎の粒状体の
移動量や回転角を求める第２ステップと、時間毎の粒状
体の移動量や回転角から各々粒状体の次の配列及び粒状
体同士の接触判定を行う第３ステップと、境界線で囲ま
れた領域を所定サイズのメッシュで区分けして、各メッ
シュ内での流体の初期流量、抵抗係数、密度を設定した
後、各メッシュ内に流入する流体の全流量を求める第４
ステップと、前記メッシュ内の粒状体部分の間隙の面積
を計算し、各メッシュの中心を結ぶ線を管路として管路
網を設定したときの、管路面積とする第５ステップと、
管路網の各節点の流量の連続性と流体の運動方程式より
各管路の流体の流速を決定する第６ステップと、各管路
の流速から管路内の粒状体位置での流速を求め、粒状体
が流体から受ける作用力を求める第７ステップと、前記
第１ステップから前記第７ステップにより、流体から粒
状体が受ける作用力により粒状体の配列状態を時間毎に
求めて粒状体の挙動を時間毎に捉えることで地盤の液状
化、クイックサンド、又は土石流を分析することを特徴
としている。
粒子の粒状体が共に存在する領域に外力が作用したとき
に、流体と粒状体が相互に干渉しながら挙動する状態を
時間毎に捉えることで地盤の液状化、クイックサンド、
又は土石流を分析する地盤分析方法において、各々粒状
体の初期配列と各々粒状体を囲む境界線を初期条件とし
て設定した後、粒状体同士の接触判定を行ない、粒状体
が相互に受ける作用力と、前記境界線で粒状体が受ける
作用力を求める第１ステップと、粒状体に作用する２つ
の作用力によって粒状体が並進或は回転する運動を表す
運動方程式を前進的に数値解析して、時間毎の粒状体の
移動量や回転角を求める第２ステップと、時間毎の粒状
体の移動量や回転角から各々粒状体の次の配列及び粒状
体同士の接触判定を行う第３ステップと、境界線で囲ま
れた領域を所定サイズのメッシュで区分けして、各メッ
シュ内での流体の初期流量、抵抗係数、密度を設定した
後、各メッシュ内に流入する流体の全流量を求める第４
ステップと、前記メッシュ内の粒状体部分の間隙の面積
を計算し、各メッシュの中心を結ぶ線を管路として管路
網を設定したときの、管路面積とする第５ステップと、
管路網の各節点の流量の連続性と流体の運動方程式より
各管路の流体の流速を決定する第６ステップと、各管路
の流速から管路内の粒状体位置での流速を求め、粒状体
が流体から受ける作用力を求める第７ステップと、前記
第１ステップから前記第７ステップにより、流体から粒
状体が受ける作用力により粒状体の配列状態を時間毎に
求めて粒状体の挙動を時間毎に捉えることで地盤の液状
化、クイックサンド、又は土石流を分析することを特徴
としている。
【０００８】
【作用】本発明は、水の流体と土粒子の粒状体が共に存
在する領域に外力が作用したときに、流体と粒状体が相
互に干渉しながら挙動する状態を時間毎に捉えることで
地盤の液状化等を分析する。 先ず、境界線で囲まれた粒
状体が相互に受ける作用力、境界線で粒状体が受ける作
用力を求め、粒状体に作用する２つの作用力によって、
粒状体が並進或は回転する運動を表す運動方程式から、
時間毎の粒状体の移動量や回転角を求め、各々粒状体の
次の配列及び粒状体同士の接触判定を行う。一方、境界
線で囲まれた領域を所定サイズのメッシュで区分けし、
各メッシュ内での流体の初期流量、抵抗係数、密度を設
定した後、各メッシュ内に流入する流体の全流量を求
め、メッシュ内の粒状体部分の間隙の面積を計算する。
そして、この面積を、各メッシュの中心を結ぶ線を管路
として管路網を設定したときの、管路面積とする。ま
た、管路網の各節点の流量の連続性と流体の運動方程式
より各管路の流体の流速を決定し、各管路の流速から管
路内の粒状体位置での流速を求め、粒状体が流体から受
ける作用力を求める。このような構成により、流体から
粒状体が受ける作用力により粒状体の配列状態を時間毎
に求めることができる。また、粒状体部分の間隙を管路
面積とすることで、２次元方向だけでなく３次元方向に
も流動する流体及び粒状体の挙動を正確につかむことが
できる。これにより、水と土粒子が共に存在する土の挙
動が正確に把握でき、地盤の液状化、クイックサンド、
又は土石流の地盤の動きを分析することができる。
在する領域に外力が作用したときに、流体と粒状体が相
互に干渉しながら挙動する状態を時間毎に捉えることで
地盤の液状化等を分析する。 先ず、境界線で囲まれた粒
状体が相互に受ける作用力、境界線で粒状体が受ける作
用力を求め、粒状体に作用する２つの作用力によって、
粒状体が並進或は回転する運動を表す運動方程式から、
時間毎の粒状体の移動量や回転角を求め、各々粒状体の
次の配列及び粒状体同士の接触判定を行う。一方、境界
線で囲まれた領域を所定サイズのメッシュで区分けし、
各メッシュ内での流体の初期流量、抵抗係数、密度を設
定した後、各メッシュ内に流入する流体の全流量を求
め、メッシュ内の粒状体部分の間隙の面積を計算する。
そして、この面積を、各メッシュの中心を結ぶ線を管路
として管路網を設定したときの、管路面積とする。ま
た、管路網の各節点の流量の連続性と流体の運動方程式
より各管路の流体の流速を決定し、各管路の流速から管
路内の粒状体位置での流速を求め、粒状体が流体から受
ける作用力を求める。このような構成により、流体から
粒状体が受ける作用力により粒状体の配列状態を時間毎
に求めることができる。また、粒状体部分の間隙を管路
面積とすることで、２次元方向だけでなく３次元方向に
も流動する流体及び粒状体の挙動を正確につかむことが
できる。これにより、水と土粒子が共に存在する土の挙
動が正確に把握でき、地盤の液状化、クイックサンド、
又は土石流の地盤の動きを分析することができる。
【０００９】
【実施例】図１に示すフローチャートに基づいて、本実
施例に係る流体と粒状体の２相系解析法を説明する。
施例に係る流体と粒状体の２相系解析法を説明する。
【００１０】この２相系解析方法は、個別要素法（ＤＥ
Ｍ）によって粒状体の挙動を計算する第１のプログラム
と、流体の挙動を計算する第２のプログラムとに分かれ
ている。
Ｍ）によって粒状体の挙動を計算する第１のプログラム
と、流体の挙動を計算する第２のプログラムとに分かれ
ている。
【００１１】すなわち、第１のプログラムでは、２相系
解析を行なうかどうかを選択し、２相系解析を行なう場
合は、各時間毎の第２のプログラムの流体解析から、流
体によって粒状体が受ける力で粒状体の次のステップで
の配列位置を決定し、第２プログラムでは、第１のプロ
グラムによって得られる粒状体の配列情報を受け取っ
て、粒状体部分での流速分布及び粒状体に及ぼす力を計
算するものである。
解析を行なうかどうかを選択し、２相系解析を行なう場
合は、各時間毎の第２のプログラムの流体解析から、流
体によって粒状体が受ける力で粒状体の次のステップで
の配列位置を決定し、第２プログラムでは、第１のプロ
グラムによって得られる粒状体の配列情報を受け取っ
て、粒状体部分での流速分布及び粒状体に及ぼす力を計
算するものである。
【００１２】以下詳細に説明すると、ステップ１０で、
図２に示すような個別要素法の要素である各粒状体４２
の初期配列及び粒状体４２の境界条件を決定する。ここ
で、境界条件を決定するとは、粒状体４２で満たされた
領域（境界線Ｂで囲まれた範囲）を設定することをい
い、粒状体４２がない範囲でも流体が満たされていれ
ば、その範囲まで行うこととする。
図２に示すような個別要素法の要素である各粒状体４２
の初期配列及び粒状体４２の境界条件を決定する。ここ
で、境界条件を決定するとは、粒状体４２で満たされた
領域（境界線Ｂで囲まれた範囲）を設定することをい
い、粒状体４２がない範囲でも流体が満たされていれ
ば、その範囲まで行うこととする。
【００１３】次に、ステップ１２で、境界線Ｂで囲まれ
た領域を接触判定用のメッシュ（Δｘ×Δｙ）で分割す
る。なお、このメッシュの分割線は、粒状体４２が移動
しても固定されたものとする。
た領域を接触判定用のメッシュ（Δｘ×Δｙ）で分割す
る。なお、このメッシュの分割線は、粒状体４２が移動
しても固定されたものとする。
【００１４】次に、ステップ１４でメッシュに囲まれた
粒状体４２同士、及び境界線Ｂにおける粒状体４２の接
触判定を行い、接触力の計算を行う。ここで、接触力の
計算は、粒状体４２の形状を２次元的に円と考え、メッ
シュ内にある粒状体４２と接触する粒状体４２の間に作
用する力は、図４及び図５に示すように、接触点の法線
方向と接線方向にばね（定数Ｋ）とダッシュポット（定
数ｎ）があると仮定すると、２要素中心間の距離の変化
によって求められる。すなわち、隣接する粒状体の半径
をＲi 、ＲJ とすると、２要素中心間の距離Ｌが、Ｌ＜
Ｒi ＋ＲJ の場合は接触状態にあり、また、Ｌ＞Ｒi ＋
ＲJ の場合は非接触状態にある。
粒状体４２同士、及び境界線Ｂにおける粒状体４２の接
触判定を行い、接触力の計算を行う。ここで、接触力の
計算は、粒状体４２の形状を２次元的に円と考え、メッ
シュ内にある粒状体４２と接触する粒状体４２の間に作
用する力は、図４及び図５に示すように、接触点の法線
方向と接線方向にばね（定数Ｋ）とダッシュポット（定
数ｎ）があると仮定すると、２要素中心間の距離の変化
によって求められる。すなわち、隣接する粒状体の半径
をＲi 、ＲJ とすると、２要素中心間の距離Ｌが、Ｌ＜
Ｒi ＋ＲJ の場合は接触状態にあり、また、Ｌ＞Ｒi ＋
ＲJ の場合は非接触状態にある。
【００１５】次に、ステップ１６で２相系の解析を行う
が否かを判断する。２相系の解析を行なわない場合、ス
テップ１８で、各粒状体４２及び境界線Ｂにおける粒状
体への作用力の決定を行う。この作用力は、並進と回転
とによって生じる次の運動方程式によって求められる。
が否かを判断する。２相系の解析を行なわない場合、ス
テップ１８で、各粒状体４２及び境界線Ｂにおける粒状
体への作用力の決定を行う。この作用力は、並進と回転
とによって生じる次の運動方程式によって求められる。
【００１６】
【数１】
【００１７】なお、図５のτ≦ｃ＋σｔａｎΦは、粒状
体の接触とは関係ない塑性平衡条件式であり、τはせん
断強さ、ｃは粒状体間に働く粘性力、σは作用面の垂直
応力、Φは内部摩擦力である。
体の接触とは関係ない塑性平衡条件式であり、τはせん
断強さ、ｃは粒状体間に働く粘性力、σは作用面の垂直
応力、Φは内部摩擦力である。
【００１８】次に、ステップ２０で各粒状体４２の加速
度を求め、ステップ２２で粒状体の位置を決定する。こ
の計算は、上記の運動方程式を前進的に数値解析するこ
とで、各粒状体４２の加速度を求め、次に、時間毎の粒
状体の移動量や回転角を求めことができる。次に、ステ
ップ２４で解析を継続するか否か判断し、解析を継続す
るなら、ステップ１４に移行する。
度を求め、ステップ２２で粒状体の位置を決定する。こ
の計算は、上記の運動方程式を前進的に数値解析するこ
とで、各粒状体４２の加速度を求め、次に、時間毎の粒
状体の移動量や回転角を求めことができる。次に、ステ
ップ２４で解析を継続するか否か判断し、解析を継続す
るなら、ステップ１４に移行する。
【００１９】また、ステップ１６で２相系の解析を行う
場合、ステップ２６で流体の諸定数を決定する。この流
体の諸定数とは、初期流量、流体の抵抗係数、流体の密
度等である。
場合、ステップ２６で流体の諸定数を決定する。この流
体の諸定数とは、初期流量、流体の抵抗係数、流体の密
度等である。
【００２０】次いで、ステップ２８で境界線Ｂで囲まれ
た領域に流入する流体の全流量Ｑを設定する（図３参
照）。ステップ３０では、各メッシュで囲まれた粒状体
部分の間隙の面積を計算し、仮想管路の管径とする。こ
の仮想管路の管径は、以下のように求められる。
た領域に流入する流体の全流量Ｑを設定する（図３参
照）。ステップ３０では、各メッシュで囲まれた粒状体
部分の間隙の面積を計算し、仮想管路の管径とする。こ
の仮想管路の管径は、以下のように求められる。
【００２１】図６（Ａ）に示すようにメッシュ（Δｘ×
Δｙ）で分割された一領域において、奥行きをΔｚと考
えると、一メッシュでの体積は、Ｖ＝Δｘ・Δｙ・Δｚ
となる。
Δｙ）で分割された一領域において、奥行きをΔｚと考
えると、一メッシュでの体積は、Ｖ＝Δｘ・Δｙ・Δｚ
となる。
【００２２】また、粒状体部分の円柱体全部の体積Ｖs
は、粒状体の半径をｒi とすれば、
は、粒状体の半径をｒi とすれば、
【００２３】
【数２】
【００２４】従って、間隙の体積Ｖe は、
【００２５】
【数３】
【００２６】一方、飽和状態において粒状体が流体に接
する全面積Ｓv は、
する全面積Ｓv は、
【００２７】
【数４】
【００２８】この時、隙間をｘ軸、ｙ軸方向へ移動する
流体の径深Ｒは、
流体の径深Ｒは、
【００２９】
【数５】
【００３０】今、図６（Ｂ）に示すような円柱体をｘ
軸、ｙ軸方向へ延びる等価な管路と考えれば、この管路
の等価な管径Ｄijは、
軸、ｙ軸方向へ延びる等価な管路と考えれば、この管路
の等価な管径Ｄijは、
【００３１】
【数６】
【００３２】次に、ステップ３２で、仮想管路網の組み
立てを以下のような方法で行なう。図７に示すように、
メッシュの中心を結ぶ線を１つの管路と考えて管路網を
設定し、管路へ流入、流出する流体の流量をｑとし、ま
た、各管路を流れる流体の流量をＱijとする。ここで、
i は回路番号、j は管路番号を示す。
立てを以下のような方法で行なう。図７に示すように、
メッシュの中心を結ぶ線を１つの管路と考えて管路網を
設定し、管路へ流入、流出する流体の流量をｑとし、ま
た、各管路を流れる流体の流量をＱijとする。ここで、
i は回路番号、j は管路番号を示す。
【００３３】隣接する２節点間の管路のエネルギー損失
をｈijとすれば、２節点間は等流と考えることができる
ので、
をｈijとすれば、２節点間は等流と考えることができる
ので、
【００３４】
【数７】
【００３５】ここで、ｒij、ｎはマニング及びヘーズン
ウイリアムスの公式による定数であり、マニングの公式
に従うものと考えると、ｒij、ｎは次の式で表される。
ウイリアムスの公式による定数であり、マニングの公式
に従うものと考えると、ｒij、ｎは次の式で表される。
【００３６】
【数８】
【００３７】だだし、Ｌijは各管路の管路長さ、Ｄijは
各管路の径であり、粗度係数ｆijは
各管路の径であり、粗度係数ｆijは
【００３８】
【数９】
【００３９】各閉回路毎に管路内の流れ方向が時計回り
を正、反時計回りを負とすると、回路一回りの損失水頭
の和は０となり、以下の式で表される。
を正、反時計回りを負とすると、回路一回りの損失水頭
の和は０となり、以下の式で表される。
【００４０】
【数１０】
【００４１】次に、ステップ３４で、ハンデイクロス法
によって管路計算を行なう。
によって管路計算を行なう。
【００４２】
【数１１】
【００４３】
【数１２】
【００４４】次に、ステップ３６で各管路の流速を決定
する。この流速Ｖijは、Ｖij＝Ｑij／Ａijで求められ
る。
する。この流速Ｖijは、Ｖij＝Ｑij／Ａijで求められ
る。
【００４５】次に、ステップ３８で流速分布を求める。
すなわち、各管路の流速から、この管路に囲まれる粒状
体位置での流速を平均化して求める。
すなわち、各管路の流速から、この管路に囲まれる粒状
体位置での流速を平均化して求める。
【００４６】次に、ステップ４０で各粒状体が流体から
受ける作用力を求める。この作用力は、流体の粘性によ
る力Ｆｖと、流体の空間的な圧力変化によるＦｐとがあ
り、次の式で求められる。
受ける作用力を求める。この作用力は、流体の粘性によ
る力Ｆｖと、流体の空間的な圧力変化によるＦｐとがあ
り、次の式で求められる。
【００４７】
【数１３】
【００４８】ここで、ＣD は抵抗係数、ρは流体の密
度、Ｖは流速である。 Ｆp ＝Ｆpx＋Ｆpy＝πｒ2 （∂Ｐ／∂ｘ＋∂Ｐ／∂ｙ） ここで、ｒは粒状体の半径、Ｐは流体の圧力である。
度、Ｖは流速である。 Ｆp ＝Ｆpx＋Ｆpy＝πｒ2 （∂Ｐ／∂ｘ＋∂Ｐ／∂ｙ） ここで、ｒは粒状体の半径、Ｐは流体の圧力である。
【００４９】この流体が粒状体に及ぼす作用力は、各時
間ステップ毎に加えられステップ１８に移行する。この
ようにして、次のステップでの粒状体の配列が決定され
る。
間ステップ毎に加えられステップ１８に移行する。この
ようにして、次のステップでの粒状体の配列が決定され
る。
【００５０】
【発明の効果】本発明は上記構成としたので、間隙中の
流体の挙動を正確に分析でき、流体から粒状体が受ける
力を正しく評価し、粒状体の挙動を正確に把握して、地
盤の動きを分析することができる。
流体の挙動を正確に分析でき、流体から粒状体が受ける
力を正しく評価し、粒状体の挙動を正確に把握して、地
盤の動きを分析することができる。
【図１】本実施例に係る流体と粒状体の２相系解析方法
の手順を示したフローチャートである。
の手順を示したフローチャートである。
【図２】モデルをメッシュ分割した状態を示した平面図
である。
である。
【図３】モデルに仮想管路網を設定した状態を示した平
面図である。
面図である。
【図４】粒状体の接触状態を示した力学モデルの概念図
である。
である。
【図５】粒状体の接触状態を示した力学モデルの概念図
である。
である。
【図６】（Ａ）はメッシュ分割の１要素の３次元的な概
念図、（Ｂ）はメッシュ分割で３次元的に管路径を示し
た概念図である。
念図、（Ｂ）はメッシュ分割で３次元的に管路径を示し
た概念図である。
【図７】仮想管路網の計算模式図である。
４２ 粒状体
フロントページの続き (58)調査した分野(Int.Cl.7，ＤＢ名) E02D 1/00 - 1/08 G06F 19/00
Claims (1)
- 【請求項１】 水の流体と土粒子の粒状体が共に存在す
る領域に外力が作用したときに、流体と粒状体が相互に
干渉しながら挙動する状態を時間毎に捉えることで地盤
の液状化、クイックサンド、又は土石流を分析する地盤
分析方法において、各々粒状体の初期配列と各々粒状体を囲む境界線 を初期
条件として設定した後、粒状体同士の接触判定を行な
い、粒状体が相互に受ける作用力と、前記境界線で粒状
体が受ける作用力を求める第１ステップと、粒状体に作用する２つの作用力によって粒状体が並進或
は回転する運動を表す 運動方程式を前進的に数値解析し
て、時間毎の粒状体の移動量や回転角を求める第２ステ
ップと、 時間毎の粒状体の移動量や回転角から各々粒状体の次の
配列及び粒状体同士の接触判定を行う第３ステップと、境界線で囲まれた領域を所定サイズのメッシュで区分け
して 、各メッシュ内での流体の初期流量、抵抗係数、密
度を設定した後、各メッシュ内に流入する流体の全流量
を求める第４ステップと、 前記メッシュ内の粒状体部分の間隙の面積を計算し、各
メッシュの中心を結ぶ線を管路として管路網を設定した
ときの、管路面積とする第５ステップと、 管路網の各節点の流量の連続性と流体の運動方程式より
各管路の流体の流速を決定する第６ステップと、 各管路の流速から管路内の粒状体位置での流速を求め、
粒状体が流体から受ける作用力を求める第７ステップ
と、 前記第１ステップから前記第７ステップにより、流体か
ら粒状体が受ける作用力により粒状体の配列状態を時間
毎に求めて粒状体の挙動を時間毎に捉えることで地盤の
液状化、クイックサンド、又は土石流を分析することを
特徴とする土の挙動を捉える地盤分析方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12021693A JP3153385B2 (ja) | 1993-05-21 | 1993-05-21 | 土の挙動を捉える地盤分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12021693A JP3153385B2 (ja) | 1993-05-21 | 1993-05-21 | 土の挙動を捉える地盤分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06330518A JPH06330518A (ja) | 1994-11-29 |
| JP3153385B2 true JP3153385B2 (ja) | 2001-04-09 |
Family
ID=14780770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP12021693A Expired - Fee Related JP3153385B2 (ja) | 1993-05-21 | 1993-05-21 | 土の挙動を捉える地盤分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3153385B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3488848B2 (ja) * | 1999-10-20 | 2004-01-19 | 株式会社 地震工学研究所 | 地震動の動的応答解析法 |
| JP4692074B2 (ja) * | 2005-05-24 | 2011-06-01 | 富士ゼロックス株式会社 | 粉体挙動解析装置及び粉体挙動解析方法、並びにコンピュータ・プログラム |
| JP4725320B2 (ja) * | 2005-12-28 | 2011-07-13 | 富士ゼロックス株式会社 | シミュレーション装置及びシミュレーション方法、画像形成装置及び画像形成方法、並びにコンピュータ・プログラム |
-
1993
- 1993-05-21 JP JP12021693A patent/JP3153385B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH06330518A (ja) | 1994-11-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nelson et al. | Flow and sediment‐transport modeling | |
| Ozdemir et al. | Evaluating scale and roughness effects in urban flood modelling using terrestrial LIDAR data | |
| Hungr et al. | Two numerical models for landslide dynamic analysis | |
| Liu et al. | Numerical simulation of debris flow with application on hazard area mapping | |
| Kang et al. | Numerical simulation of 2D granular flow entrainment using DEM | |
| Tayyebi et al. | Two-phase SPH numerical study of pore-water pressure effect on debris flows mobility: Yu Tung debris flow | |
| Pastor et al. | A depth integrated, coupled, two-phase model for debris flow propagation | |
| Manzanal et al. | Application of a new rheological model to rock avalanches: an SPH approach | |
| Shen et al. | A modified finite difference model for the modeling of flowslides | |
| Abe et al. | Numerical simulation for runout process of debris flow using depth-averaged material point method | |
| Hromadka II et al. | A two-dimensional dam-break flood plain model | |
| Jiang et al. | Experimental validation of a new semi-empirical impact force model of the dry granular flow impact against a rigid barrier | |
| Sarker | Essence of mike 21c (fdm numerical scheme): Application on the river morphology of bangladesh | |
| JP3153385B2 (ja) | 土の挙動を捉える地盤分析方法 | |
| Ceccato et al. | The effect of the front inclination on the impact forces transmitted by granular flows to rigid structures | |
| Youngs | Soil physical theory and heterogeneity | |
| Issakhov et al. | Numerical simulation of the water surface movement with macroscopic particles on movable beds | |
| KR101813903B1 (ko) | 대변형 해석기법을 이용한 사방댐 설계 시스템 및 방법 | |
| Romanoa et al. | Significance of unsaturated flow and seepage faces in the simulation of steady‐state subsurface flow | |
| Flora et al. | Uncertainty quantification of large-eddy simulation results of riverine flows: a field and numerical study | |
| Fathani et al. | A numerical model for the analysis of rapid landslide motion | |
| Khairi et al. | Smoothed particle hydrodynamics simulation for debris flow: a review | |
| Fathi et al. | Rapid strain demand estimation of pipelines deformed by lateral ground movements | |
| JP3960425B2 (ja) | 土砂災害の氾濫解析システム | |
| Gill et al. | Extending the strain path method analogy for modelling penetrometer installation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |