粒子欠損を考慮した個別要素法
【課題】個別要素法において、簡易な粒子欠損モデルにより、土粒子の破砕の影響を再現する。
【解決手段】連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、前記個別要素を土粒子として見立て、これら個別要素iが受ける接触力の最大値Fiが設定された限界荷重Fuを超えたときに、個別要素iの半径Riを所定量だけ縮小させる(縮小後半径:R*i)粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重Fuを段階的に大きく設定し、個別要素の半径縮小比率[R*i/Ri]を段階的に小さく設定する。
【解決手段】連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、前記個別要素を土粒子として見立て、これら個別要素iが受ける接触力の最大値Fiが設定された限界荷重Fuを超えたときに、個別要素iの半径Riを所定量だけ縮小させる(縮小後半径:R*i)粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重Fuを段階的に大きく設定し、個別要素の半径縮小比率[R*i/Ri]を段階的に小さく設定する。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、連続体を個別要素に分割して離散モデル化し、その挙動を解析する個別要素法において、簡易な粒子欠損モデルにより、土粒子の破砕の影響を再現し得る個別要素法に関する。
【背景技術】
【0002】
従来より、静的荷重や動的荷重を受ける連続体の挙動を解析する手法として、連続体を有限の単純な要素の集合体にモデル化して数値解析を行う有限要素法が用いられている。この有限要素法による数値解析は、連続体の微小変形や微小振動などに対しては精度よく再現できるが、大変形や破壊現象をシミュレーションすることは困難であった。
【0003】
そのため、近年では、かかる連続体の大変形や破壊などのシミュレーションについては、連続体を離散可能な個別要素に分割して解析を行う個別要素法が用いられるようになっている。前記個別要素法は、連続体を複数の個別要素に分割し、各要素間に作用するバネ作用などを考慮した離散モデルにモデル化して数値解析を行う手法である。前記離散モデルについて、図12に示されるように、任意の個別要素50、51を取り出して具体的に詳述すると、これら個別要素50、51間の接続は、個別要素50、51同士の接触によるバネ作用や粘性作用が、バネ52及びダッシュポット53でモデル化されている。また、個別要素50、51間の接続は、図13に示されるように、個別要素50,51に引張荷重Tが作用し、所定の引張荷重Ty以上となって、両要素間距離σが2つの要素の半径の和に相当する距離σuより大きくなったとき、両要素が分離して相互に無関係となる状態が、ディバイダー54でモデル化されている。なお、前記引張荷重Tyまでの前記要素間距離σは、図示例のように線形となっている。
【0004】
前記個別要素法では、実際の土に比べてせん断初期圧縮しにくく、その後のダイラタンシーは大きいという問題があった。その結果、後述の実施に示されるように、地震によって沈下が生じた地盤を個別要素法によって解析すると、初期のモデル状態で生じている先行圧縮応力が開放されたためと思われるが、地盤沈下とは逆に地盤が隆起する結果となることがあった。この原因は、現状の個別要素法では粒子の破砕による影響が考慮されていないためであると思われる。
【0005】
一方、下記非特許文献1では、本来の地盤材料を個別要素法によって解析するには、粒子破砕のような粒子の塑性変形を取り入れる必要があるとし、図14に示されるように、土粒子を球要素のかたまり60として表現するとともに、幾つかの球要素間に結合強度(引張力及びせん断力に対して)を与える一方で、残りの球要素間にはすべり強度を与え、ある限界値を超えた場合には結合が解かれ(破砕)、土粒子の形状に変化が起こるとした個別要素法が提案されている。この文献では、図15に示されるように、前述した土粒子モデルを389個用いて集合化したモデル61によって、等方圧縮試験や三軸圧縮試験などの力学試験をシミュレーションしている。
【先行技術文献】
【非特許文献】
【0006】
【非特許文献1】中田幸男、M.D.Bolton、”破砕性土の個別要素法シミュレーション”、第37回地盤工学会研究発表会、2002.7
【発明の概要】
【発明が解決しようとする課題】
【0007】
上記非特許文献1に係る解析手法のように、土粒子を球要素の固まり60として表現すれば破砕の影響をシミュレーションできるものと考えられるが、実際の地盤沈下等をこの個別要素法で解析する場合には、要素数が膨大となりモデルの作成が煩雑であるとともに、計算に多くの時間を要するなどの問題があった。
【0008】
そこで本発明の主たる課題は、連続体を個別要素に分割して離散モデル化し、その挙動を解析する個別要素法において、簡易な粒子欠損モデルにより、土粒子の破砕の影響を再現し得るようにすることにある。
【課題を解決するための手段】
【0009】
前記課題を解決するために請求項1に係る本発明として、連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定することを特徴とする粒子欠損を考慮した個別要素法が提供される。
【0010】
上記請求項1記載の本発明では、各個別要素を土粒子として見立て、これに他の土粒子や境界から作用する圧縮荷重(法線方向接触力)が限界荷重を超えると、個別要素の半径を所定量だけ縮小させる粒子欠損操作を行うものである。つまり、設定された限界荷重を超えた際に、強制的に個別要素の半径を所定量だけ縮小させれば、要素の体積減少によって、土粒子の破砕の影響が簡易的に再現されるとの考え方を導入したものである。
【0011】
また、後述する実施例1に記載されるように、砕石試験体を圧縮荷重によって破砕した実験によれば、初期の小さな荷重で小さく欠損し、次の荷重で大きく破砕するものが多かった点に鑑み、前記粒子欠損操作は、複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定するするようにした。すなわち、最初は小さい荷重で小さく破壊し、順次大きな荷重で大きく破壊が生じる現象に近づけることで精度の向上を図ったものである。
【0012】
請求項2に係る本発明として、前記粒子欠損操作を2段階とする請求項1記載の粒子欠損を考慮した個別要素法が提供される。
【0013】
上記請求項2記載の発明は、解析を省力化するとともに、計算時間の短縮を図るために、前記粒子欠損操作を2段階とするものである。
【0014】
請求項3に係る本発明として、前記粒子欠損操作は、下式(12)により行う請求項1,2いずれかに記載の粒子欠損を考慮した個別要素法が提供される。
【数12】
【発明の効果】
【0015】
以上詳説のとおり本発明によれば、連続体を個別要素に分割して離散モデル化し、その挙動を解析する個別要素法において、簡易な粒子欠損モデルにより、土粒子の破砕の影響を再現し得るようになる。
【図面の簡単な説明】
【0016】
【図1】本発明に係る個別要素法のフローチャート図である。
【図2】本発明に係る個別要素法の解析原理説明図(その1)である。
【図3】本発明に係る個別要素法の解析原理説明図(その2)である。
【図4】本発明に係る個別要素法の解析原理説明図(その3)である。
【図5】砕石の破砕実験結果を示す図である。
【図6】三軸圧縮試験をシミュレーションする個別要素法の解析モデル図である。
【図7】実施例2における比較例の解析結果を示す図である。
【図8】実施例2における本発明解析結果を示す図である。
【図9】実施例3における地盤沈下解析対象の図である。
【図10】実施例3における解析の入力地震波(加速度、速度、変位)である。
【図11】実施例3における解析結果を示す図である。
【図12】個別要素法における要素間のモデル図[(A)は法線方向、(B)は接線方向]である。
【図13】個別要素50,51間の引張荷重Tと要素間距離σの関係を示すグラフである。
【0017】
【図14】非特許文献1における土粒子モデル図である。
【図15】非特許文献1における土粒子の集合化モデル図である。
【発明を実施するための形態】
【0018】
以下、本発明の実施の形態について図面を参照しながら詳述する。
【0019】
本発明は、連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きく設定し、個別要素の半径縮小比率を段階的に小さく設定することにより、簡易に土粒子の破砕の影響を再現し得るようにしたものである。
【0020】
具体的な手順は、図1のフローチャート図の内、破線で囲んだ部分が本発明の粒子欠損操作部分であり、その他の部分は一般的な個別要素法のアルゴリズムである。
(1)要素の運動方程式
質量mi、慣性モーメントIiのある要素iについての運動方程式は次式で表される。
【数1】
【0021】
(2)接触間の接触判定
要素間の相互作用が働く条件として、要素間の接触判定を行う。時刻tにおける要素iと要素jの中心の座標をそれぞれ([X]ti、[Y]ti)、([X]tj、[Y]tj)とし、回転変位量を([θ]ti、[θ]tj)、半径をそれぞれri,rjとした時,下式(2)を満たすと要素iと要素jは接触したものとする。
【数2】
【0022】
(3)要素間接触力(法線方向)の計算
図12に示されるように、接触した要素間には相互作用力が働く相互作用力の粒子間に弾性バネ52と粘性ダッシュポッド53を仮定している。
【0023】
弾性バネ52による相互作用力は下式(3)で表される。
【数3】
【0024】
また、粘性ダッシュポッド53による相互作用力は下式(4)で表される。
【数4】
【0025】
従って、法線方向に働く力は下式(5)で示される。
【数5】
【0026】
本発明において、接触力の最大値が限界荷重を超えるか否かの判定を行う接触力は、上記法線方向に作用する荷重[Fn]tである。
【0027】
(4)要素間接触力(接線方向)の計算
図12に示されるように、弾性ばね52によって作用する接線方向の力は,接触した時点での値を0とし下式(6)に従い接線方向力を要素が接触している間,毎ステップ増減させることによって計算を行う。
【数6】
【0028】
また、滑りが生じた場合には,下式(7)で示される接線力が働く。
【数7】
【0029】
粘性ダッシュポッド53によって接線方向に働く力は,下式(8)で表される。
【数8】
【0030】
結局,接線方向に働く力は下式(9)で示される。
【数9】
【0031】
(5)粒子欠損操作について
本発明では、個別要素を土粒子として見立て、これら個別要素iが受ける接触力の最大値Fiが限界荷重Fuを超えたときに、個別要素iの半径Riを所定量だけ縮小させるものである(縮小後半径:R*i)。この際、土粒子の破砕形態は、1回の欠損で破砕が完了するわけではなく、欠損が何度か繰り返されるとともに、この破砕は初期の小さな荷重では小さく欠損し、次の大きな荷重で大きく破砕する傾向を示すとの知見の下、前記粒子欠損操作を、接触力Fiが漸次大きくなる過程で複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重Fuを段階的に大きく設定し、個別要素iの半径縮小比率[R*i/Ri]を段階的に小さく設定するものである。
【0032】
以下、更に具体的に説明を行う。
【0033】
図に示されるように、時間tにおいて、個別要素iが周辺の要素j,k,l,m等から接触荷重(法線方向接触力)Fji(t),Fki(t),Fli(t),Fmi(t)を受けているとする。この力は,互いの要素の位置関係からバネの縮みを求めて計算される。
【0034】
また,周辺の要素j,k,l,m等に対しては反作用があり、これらの力の不釣り合い力により,要素iは加速度を得てこれまでの速度を増加させ、次のステップで,別の場所へ移動する。周辺の要素j,k,l,m等も同様に別の場所へ移動する。すると、次のステップ時間t+Δtでは,要素iは,前ステップとは異なる力Fji(t+Δt), Fki(t+Δt),Fli(t+Δt) ,Fmi(t+Δt)を受けることになる。このように,要素iは,周囲のいくつかの要素から,時々刻々変化する力を受ける。
【0035】
個別要素iが、時刻(t)に周囲の要素から受ける接触力Fji(t), Fki(t),Fli(t) ,Fmi(t)の内、最大値をFi(t)とし、この接触力Fi(t)が限界荷重Fu(j)を超えたとき、下式(10)に示されるように、個別要素iの半径Riを所定量だけ縮小させる粒子欠損操作を行うようにする。なお、限界荷重Fuの添字(j)は粒子欠損操作回数時を表す。
【数10】
【0036】
この際、粒子欠損操作を終えた要素iについては、次回の欠損時の限界荷重Fu(j+1)を下式(11)に設定する。
【数11】
【0037】
以上詳述した粒子欠損操作を纏めて一般式で表すと下式(12)となる。
【数12】
【0038】
粒子欠損操作は、複数回段階的に繰り返し行われる。前記粒子欠損操作は、回数の限度なく行うことも可能であるが、計算が煩雑になるため、実用上は2段階で一定の精度を確保することが可能である。また、限界荷重係数au、半径を縮小する際の減少分の基本比率α は、土の物性に合わせて試行的に決定される。
【実施例】
【0039】
〔実施例1〕
実施例1では、砕石の欠損を個別要素法に考慮するための基礎実験として、砕石の破砕実験を行った。その結果を図5に示す。試験は、アムスラーで試験体を押し潰す方法で行った。砕石は初期の小さな荷重で小さく欠損し、次の大きな荷重で大きく欠損するものが多かった。図中、●は局所欠損を示し、○が大きな破壊を示す。
【0040】
上記実験は砕石を対象としているものであるが、土粒子レベルでも同様の現象が再現されるとの仮定の下、「最初は小さい荷重で小さく破壊し、次に大きな荷重で大きな破壊が生じる」現象を個別要素法に取り入れるべく、この粒子欠損の特徴を表現する式の1つの例として、上式(12)を導いた。
【0041】
〔実施例2〕
実施例2では、原地盤材料に対する三軸圧縮試験を対象に本発明のシミュレーション検討を行った。
【0042】
三軸圧縮試験の拘束圧は、50kPa、100kPa、300kPaの3種類とした。また、個別要素法の解析モデルを図6に示す。また、解析パラメータを下表1に示す。粒子欠損のパラメータは、下表2のとおりとした。
【0043】
なお、比較例として粒子欠損操作を行わない通常の個別要素法解析を行った。
【0044】
【表1】
【0045】
【表2】
【0046】
比較例の個別要素法解析結果を図7に示し、本発明の個別要素法解析結果を図8に示す。また、各グラフ中に三軸圧縮試験結果をプロットした。
【0047】
比較例ではせん断初期において原地盤よりも圧縮性がかなり小さくなるという傾向が見られたのに対して、本発明の場合はこの傾向がかなり改善されていることが分かる。
【0048】
〔実施例3〕
実際の地震によって地盤沈下が生じた地盤を対象として、本発明の個別要素法を適用してシミュレーションを行った。対象は図9に示されるように、建屋に挟まれた35m区間の地盤である。また、入力地震波を図10に示した。更に、個別要素法における地盤構成要素諸元を下表3に示し、解析パラメータを下表4に示す。なお、初期最大接触力の平均値Foは、地盤を10mの深度毎に分割し、それぞれの領域で求めた。もちろん、全深度の平均値としてもよい。
【0049】
【表3】
【0050】
【表4】
【0051】
図11に個別要素法によるシミュレーション結果(本発明法)を実際に生じた地盤沈下及び粒子欠損操作を行わない通常の個別要素法解析の結果(比較例)と共に示す。
【0052】
図11において、粒子欠損操作を行わない比較例の場合は、沈下せずに逆に膨張する結果となっている。これは、初期のモデル状態で生じている先行圧縮荷重が逆に開放されたためであると思われる。
【0053】
これに対して、本発明法の場合は、建屋近傍の地盤沈下は少なめにシミュレートされているが、建屋から離れるに従ってかなり精度良くシミュレートされている。従来の個別要素法の場合は逆に膨張していることを考えれば、建屋近傍で誤差が生じているものの十分に実用に供するものと思われる。
【符号の説明】
【0054】
1…個別要素
【技術分野】
【0001】
本発明は、連続体を個別要素に分割して離散モデル化し、その挙動を解析する個別要素法において、簡易な粒子欠損モデルにより、土粒子の破砕の影響を再現し得る個別要素法に関する。
【背景技術】
【0002】
従来より、静的荷重や動的荷重を受ける連続体の挙動を解析する手法として、連続体を有限の単純な要素の集合体にモデル化して数値解析を行う有限要素法が用いられている。この有限要素法による数値解析は、連続体の微小変形や微小振動などに対しては精度よく再現できるが、大変形や破壊現象をシミュレーションすることは困難であった。
【0003】
そのため、近年では、かかる連続体の大変形や破壊などのシミュレーションについては、連続体を離散可能な個別要素に分割して解析を行う個別要素法が用いられるようになっている。前記個別要素法は、連続体を複数の個別要素に分割し、各要素間に作用するバネ作用などを考慮した離散モデルにモデル化して数値解析を行う手法である。前記離散モデルについて、図12に示されるように、任意の個別要素50、51を取り出して具体的に詳述すると、これら個別要素50、51間の接続は、個別要素50、51同士の接触によるバネ作用や粘性作用が、バネ52及びダッシュポット53でモデル化されている。また、個別要素50、51間の接続は、図13に示されるように、個別要素50,51に引張荷重Tが作用し、所定の引張荷重Ty以上となって、両要素間距離σが2つの要素の半径の和に相当する距離σuより大きくなったとき、両要素が分離して相互に無関係となる状態が、ディバイダー54でモデル化されている。なお、前記引張荷重Tyまでの前記要素間距離σは、図示例のように線形となっている。
【0004】
前記個別要素法では、実際の土に比べてせん断初期圧縮しにくく、その後のダイラタンシーは大きいという問題があった。その結果、後述の実施に示されるように、地震によって沈下が生じた地盤を個別要素法によって解析すると、初期のモデル状態で生じている先行圧縮応力が開放されたためと思われるが、地盤沈下とは逆に地盤が隆起する結果となることがあった。この原因は、現状の個別要素法では粒子の破砕による影響が考慮されていないためであると思われる。
【0005】
一方、下記非特許文献1では、本来の地盤材料を個別要素法によって解析するには、粒子破砕のような粒子の塑性変形を取り入れる必要があるとし、図14に示されるように、土粒子を球要素のかたまり60として表現するとともに、幾つかの球要素間に結合強度(引張力及びせん断力に対して)を与える一方で、残りの球要素間にはすべり強度を与え、ある限界値を超えた場合には結合が解かれ(破砕)、土粒子の形状に変化が起こるとした個別要素法が提案されている。この文献では、図15に示されるように、前述した土粒子モデルを389個用いて集合化したモデル61によって、等方圧縮試験や三軸圧縮試験などの力学試験をシミュレーションしている。
【先行技術文献】
【非特許文献】
【0006】
【非特許文献1】中田幸男、M.D.Bolton、”破砕性土の個別要素法シミュレーション”、第37回地盤工学会研究発表会、2002.7
【発明の概要】
【発明が解決しようとする課題】
【0007】
上記非特許文献1に係る解析手法のように、土粒子を球要素の固まり60として表現すれば破砕の影響をシミュレーションできるものと考えられるが、実際の地盤沈下等をこの個別要素法で解析する場合には、要素数が膨大となりモデルの作成が煩雑であるとともに、計算に多くの時間を要するなどの問題があった。
【0008】
そこで本発明の主たる課題は、連続体を個別要素に分割して離散モデル化し、その挙動を解析する個別要素法において、簡易な粒子欠損モデルにより、土粒子の破砕の影響を再現し得るようにすることにある。
【課題を解決するための手段】
【0009】
前記課題を解決するために請求項1に係る本発明として、連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定することを特徴とする粒子欠損を考慮した個別要素法が提供される。
【0010】
上記請求項1記載の本発明では、各個別要素を土粒子として見立て、これに他の土粒子や境界から作用する圧縮荷重(法線方向接触力)が限界荷重を超えると、個別要素の半径を所定量だけ縮小させる粒子欠損操作を行うものである。つまり、設定された限界荷重を超えた際に、強制的に個別要素の半径を所定量だけ縮小させれば、要素の体積減少によって、土粒子の破砕の影響が簡易的に再現されるとの考え方を導入したものである。
【0011】
また、後述する実施例1に記載されるように、砕石試験体を圧縮荷重によって破砕した実験によれば、初期の小さな荷重で小さく欠損し、次の荷重で大きく破砕するものが多かった点に鑑み、前記粒子欠損操作は、複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定するするようにした。すなわち、最初は小さい荷重で小さく破壊し、順次大きな荷重で大きく破壊が生じる現象に近づけることで精度の向上を図ったものである。
【0012】
請求項2に係る本発明として、前記粒子欠損操作を2段階とする請求項1記載の粒子欠損を考慮した個別要素法が提供される。
【0013】
上記請求項2記載の発明は、解析を省力化するとともに、計算時間の短縮を図るために、前記粒子欠損操作を2段階とするものである。
【0014】
請求項3に係る本発明として、前記粒子欠損操作は、下式(12)により行う請求項1,2いずれかに記載の粒子欠損を考慮した個別要素法が提供される。
【数12】
【発明の効果】
【0015】
以上詳説のとおり本発明によれば、連続体を個別要素に分割して離散モデル化し、その挙動を解析する個別要素法において、簡易な粒子欠損モデルにより、土粒子の破砕の影響を再現し得るようになる。
【図面の簡単な説明】
【0016】
【図1】本発明に係る個別要素法のフローチャート図である。
【図2】本発明に係る個別要素法の解析原理説明図(その1)である。
【図3】本発明に係る個別要素法の解析原理説明図(その2)である。
【図4】本発明に係る個別要素法の解析原理説明図(その3)である。
【図5】砕石の破砕実験結果を示す図である。
【図6】三軸圧縮試験をシミュレーションする個別要素法の解析モデル図である。
【図7】実施例2における比較例の解析結果を示す図である。
【図8】実施例2における本発明解析結果を示す図である。
【図9】実施例3における地盤沈下解析対象の図である。
【図10】実施例3における解析の入力地震波(加速度、速度、変位)である。
【図11】実施例3における解析結果を示す図である。
【図12】個別要素法における要素間のモデル図[(A)は法線方向、(B)は接線方向]である。
【図13】個別要素50,51間の引張荷重Tと要素間距離σの関係を示すグラフである。
【0017】
【図14】非特許文献1における土粒子モデル図である。
【図15】非特許文献1における土粒子の集合化モデル図である。
【発明を実施するための形態】
【0018】
以下、本発明の実施の形態について図面を参照しながら詳述する。
【0019】
本発明は、連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きく設定し、個別要素の半径縮小比率を段階的に小さく設定することにより、簡易に土粒子の破砕の影響を再現し得るようにしたものである。
【0020】
具体的な手順は、図1のフローチャート図の内、破線で囲んだ部分が本発明の粒子欠損操作部分であり、その他の部分は一般的な個別要素法のアルゴリズムである。
(1)要素の運動方程式
質量mi、慣性モーメントIiのある要素iについての運動方程式は次式で表される。
【数1】
【0021】
(2)接触間の接触判定
要素間の相互作用が働く条件として、要素間の接触判定を行う。時刻tにおける要素iと要素jの中心の座標をそれぞれ([X]ti、[Y]ti)、([X]tj、[Y]tj)とし、回転変位量を([θ]ti、[θ]tj)、半径をそれぞれri,rjとした時,下式(2)を満たすと要素iと要素jは接触したものとする。
【数2】
【0022】
(3)要素間接触力(法線方向)の計算
図12に示されるように、接触した要素間には相互作用力が働く相互作用力の粒子間に弾性バネ52と粘性ダッシュポッド53を仮定している。
【0023】
弾性バネ52による相互作用力は下式(3)で表される。
【数3】
【0024】
また、粘性ダッシュポッド53による相互作用力は下式(4)で表される。
【数4】
【0025】
従って、法線方向に働く力は下式(5)で示される。
【数5】
【0026】
本発明において、接触力の最大値が限界荷重を超えるか否かの判定を行う接触力は、上記法線方向に作用する荷重[Fn]tである。
【0027】
(4)要素間接触力(接線方向)の計算
図12に示されるように、弾性ばね52によって作用する接線方向の力は,接触した時点での値を0とし下式(6)に従い接線方向力を要素が接触している間,毎ステップ増減させることによって計算を行う。
【数6】
【0028】
また、滑りが生じた場合には,下式(7)で示される接線力が働く。
【数7】
【0029】
粘性ダッシュポッド53によって接線方向に働く力は,下式(8)で表される。
【数8】
【0030】
結局,接線方向に働く力は下式(9)で示される。
【数9】
【0031】
(5)粒子欠損操作について
本発明では、個別要素を土粒子として見立て、これら個別要素iが受ける接触力の最大値Fiが限界荷重Fuを超えたときに、個別要素iの半径Riを所定量だけ縮小させるものである(縮小後半径:R*i)。この際、土粒子の破砕形態は、1回の欠損で破砕が完了するわけではなく、欠損が何度か繰り返されるとともに、この破砕は初期の小さな荷重では小さく欠損し、次の大きな荷重で大きく破砕する傾向を示すとの知見の下、前記粒子欠損操作を、接触力Fiが漸次大きくなる過程で複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重Fuを段階的に大きく設定し、個別要素iの半径縮小比率[R*i/Ri]を段階的に小さく設定するものである。
【0032】
以下、更に具体的に説明を行う。
【0033】
図に示されるように、時間tにおいて、個別要素iが周辺の要素j,k,l,m等から接触荷重(法線方向接触力)Fji(t),Fki(t),Fli(t),Fmi(t)を受けているとする。この力は,互いの要素の位置関係からバネの縮みを求めて計算される。
【0034】
また,周辺の要素j,k,l,m等に対しては反作用があり、これらの力の不釣り合い力により,要素iは加速度を得てこれまでの速度を増加させ、次のステップで,別の場所へ移動する。周辺の要素j,k,l,m等も同様に別の場所へ移動する。すると、次のステップ時間t+Δtでは,要素iは,前ステップとは異なる力Fji(t+Δt), Fki(t+Δt),Fli(t+Δt) ,Fmi(t+Δt)を受けることになる。このように,要素iは,周囲のいくつかの要素から,時々刻々変化する力を受ける。
【0035】
個別要素iが、時刻(t)に周囲の要素から受ける接触力Fji(t), Fki(t),Fli(t) ,Fmi(t)の内、最大値をFi(t)とし、この接触力Fi(t)が限界荷重Fu(j)を超えたとき、下式(10)に示されるように、個別要素iの半径Riを所定量だけ縮小させる粒子欠損操作を行うようにする。なお、限界荷重Fuの添字(j)は粒子欠損操作回数時を表す。
【数10】
【0036】
この際、粒子欠損操作を終えた要素iについては、次回の欠損時の限界荷重Fu(j+1)を下式(11)に設定する。
【数11】
【0037】
以上詳述した粒子欠損操作を纏めて一般式で表すと下式(12)となる。
【数12】
【0038】
粒子欠損操作は、複数回段階的に繰り返し行われる。前記粒子欠損操作は、回数の限度なく行うことも可能であるが、計算が煩雑になるため、実用上は2段階で一定の精度を確保することが可能である。また、限界荷重係数au、半径を縮小する際の減少分の基本比率α は、土の物性に合わせて試行的に決定される。
【実施例】
【0039】
〔実施例1〕
実施例1では、砕石の欠損を個別要素法に考慮するための基礎実験として、砕石の破砕実験を行った。その結果を図5に示す。試験は、アムスラーで試験体を押し潰す方法で行った。砕石は初期の小さな荷重で小さく欠損し、次の大きな荷重で大きく欠損するものが多かった。図中、●は局所欠損を示し、○が大きな破壊を示す。
【0040】
上記実験は砕石を対象としているものであるが、土粒子レベルでも同様の現象が再現されるとの仮定の下、「最初は小さい荷重で小さく破壊し、次に大きな荷重で大きな破壊が生じる」現象を個別要素法に取り入れるべく、この粒子欠損の特徴を表現する式の1つの例として、上式(12)を導いた。
【0041】
〔実施例2〕
実施例2では、原地盤材料に対する三軸圧縮試験を対象に本発明のシミュレーション検討を行った。
【0042】
三軸圧縮試験の拘束圧は、50kPa、100kPa、300kPaの3種類とした。また、個別要素法の解析モデルを図6に示す。また、解析パラメータを下表1に示す。粒子欠損のパラメータは、下表2のとおりとした。
【0043】
なお、比較例として粒子欠損操作を行わない通常の個別要素法解析を行った。
【0044】
【表1】
【0045】
【表2】
【0046】
比較例の個別要素法解析結果を図7に示し、本発明の個別要素法解析結果を図8に示す。また、各グラフ中に三軸圧縮試験結果をプロットした。
【0047】
比較例ではせん断初期において原地盤よりも圧縮性がかなり小さくなるという傾向が見られたのに対して、本発明の場合はこの傾向がかなり改善されていることが分かる。
【0048】
〔実施例3〕
実際の地震によって地盤沈下が生じた地盤を対象として、本発明の個別要素法を適用してシミュレーションを行った。対象は図9に示されるように、建屋に挟まれた35m区間の地盤である。また、入力地震波を図10に示した。更に、個別要素法における地盤構成要素諸元を下表3に示し、解析パラメータを下表4に示す。なお、初期最大接触力の平均値Foは、地盤を10mの深度毎に分割し、それぞれの領域で求めた。もちろん、全深度の平均値としてもよい。
【0049】
【表3】
【0050】
【表4】
【0051】
図11に個別要素法によるシミュレーション結果(本発明法)を実際に生じた地盤沈下及び粒子欠損操作を行わない通常の個別要素法解析の結果(比較例)と共に示す。
【0052】
図11において、粒子欠損操作を行わない比較例の場合は、沈下せずに逆に膨張する結果となっている。これは、初期のモデル状態で生じている先行圧縮荷重が逆に開放されたためであると思われる。
【0053】
これに対して、本発明法の場合は、建屋近傍の地盤沈下は少なめにシミュレートされているが、建屋から離れるに従ってかなり精度良くシミュレートされている。従来の個別要素法の場合は逆に膨張していることを考えれば、建屋近傍で誤差が生じているものの十分に実用に供するものと思われる。
【符号の説明】
【0054】
1…個別要素
【特許請求の範囲】
【請求項1】
連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定することを特徴とする粒子欠損を考慮した個別要素法。
【請求項2】
前記粒子欠損操作を2段階とする請求項1記載の粒子欠損を考慮した個別要素法。
【請求項3】
前記粒子欠損操作は、下式(12)により行う請求項1,2いずれかに記載の粒子欠損を考慮した個別要素法。
【数12】
【請求項1】
連続体である被対象物を個別要素に分割し、前記個別要素間のバネ作用を考慮した離散モデルにモデル化し解析を行う個別要素法において、
前記個別要素を土粒子として見立て、これら個別要素が受ける接触力の最大値が設定された限界荷重を超えたときに、個別要素の半径を所定量だけ縮小させる粒子欠損操作を複数回段階的に行うとともに、粒子欠損操作毎に前記限界荷重を段階的に大きくし、個別要素の半径縮小比率を段階的に小さく設定することを特徴とする粒子欠損を考慮した個別要素法。
【請求項2】
前記粒子欠損操作を2段階とする請求項1記載の粒子欠損を考慮した個別要素法。
【請求項3】
前記粒子欠損操作は、下式(12)により行う請求項1,2いずれかに記載の粒子欠損を考慮した個別要素法。
【数12】
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【公開番号】特開2011−248826(P2011−248826A)
【公開日】平成23年12月8日(2011.12.8)
【国際特許分類】
【出願番号】特願2010−124150(P2010−124150)
【出願日】平成22年5月31日(2010.5.31)
【出願人】(000003687)東京電力株式会社 (2,580)
【出願人】(000221546)東電設計株式会社 (44)
【公開日】平成23年12月8日(2011.12.8)
【国際特許分類】
【出願日】平成22年5月31日(2010.5.31)
【出願人】(000003687)東京電力株式会社 (2,580)
【出願人】(000221546)東電設計株式会社 (44)
[ Back to top ]