個別要素法のソースを表示

'''個別要素法'''（こべつようそほう、{{lang-en-short|Distinct Element Method}}、DEM）または'''離散要素法'''（{{lang|en|Discrete Element Method}}、DEM）は、解析の対象を自由に運動できる多角形や円形・球の要素の集合体としてモデル化し、要素間の接触・滑動を考慮して、各時刻におけるそれぞれの要素の運動を逐次追跡して解析する手法である。もとは岩盤工学に適用するためにPeter A. Cundall (1971)およびCundall and Strack (1979)<ref>P. A. Cundall, O. D. L. Strack, "A discrete numerical model for granular assembles," Geotechnique 29 (1979) 47-65.</ref>により発表された論文に端を発しており、現在は液状化や土石流など地盤の挙動解析や[[コンクリート]]構造物、[[粉体]]（[[化学工学]]、[[リチウムイオン電池]]、[[薬学]]、[[農学]]など）、磁気相互作用力を有する電子写真システムの[[トナー]]の挙動解析などに用いられている。

== 概略 ==
以下に、円形要素を用いた際の運動方程式を示す。

質量<math>m_i</math>、慣性モーメント<math>I_i</math>のある円形要素<math>i</math>について、次の運動方程式が成り立つ。
:<math>m_i \ddot{\textbf{u}}+C_i \dot{\textbf{u}}+F_i =0 </math>
:<math>I_i \ddot{ \phi }+D_i \dot{ \phi }+M_i =0 </math>

ここに<math>F_i</math>：要素に働く合力、<math>M_i</math>：要素に働く合モーメント、<math>C_i</math>、<math>D_i</math>：減衰定数、<math>\textbf{u}</math>：要素の変位ベクトル、<math>\phi</math>：要素の回転変位である。

要素同士が接触しているときは<math>F_i = K \textbf{u}</math>（<math>K</math>は弾性定数）及び<math>M_i = K r^2</math>（<math>r</math>は円の半径）、離れているときは、<math>F_i = M_i = 0</math>で表される。ただし、重力を考える場合は合力の項で考慮する。

上式を数値積分することで、逐次変位ベクトルと回転変位を得ることができる。

== 商用コード ==

* Granuleworks
* EDEM
* Rocky DEM

== 参考文献 ==
{{reflist}}
* {{cite|和書 |author=伯野元彦 |title=破壊のシミュレーション －拡張個別要素法で破壊を追う－ |publisher=森北出版 |year=1997 |isbn=}}
* {{cite|和書 |author=酒井幹夫 |title=粉体の数値シミュレーション |publisher=丸善出版 |year=2012 |isbn=978-4-621-08582-0 |pages=1-11}}
* Catherine O'Sullivan, 鈴木 輝一 (訳):「粒子個別要素法」、森北出版、ISBN:978-4627915817、（2014年5月27日発売）


{{DEFAULTSORT:こへつようそほう}}
[[Category:数値微分方程式]]
[[Category:数学に関する記事]]