個別要素法
ここに F i {\displaystyle F_{i}} :要素に働く合力、 M i {\displaystyle M_{i}} :要素に働く合モーメント、 C i {\displaystyle C_{i}} 、 D i {\displaystyle D_{i}} :減衰定数、 u {\displaystyle {\textbf {u}}} :要素の変位ベクトル、 φ {\displaystyle \phi } :要素の回転変位である。
要素同士が接触しているときは F i = K u {\displaystyle F_{i}=K{\textbf {u}}} ( K {\displaystyle K} は弾性定数)及び M i = K r 2 {\displaystyle M_{i}=Kr^{2}} ( r {\displaystyle r} は円の半径)、離れているときは、 F i = M i = 0 {\displaystyle F_{i}=M_{i}=0} で表される。ただし、重力を考える場合は合力の項で考慮する。
上式を数値積分することで、逐次変位ベクトルと回転変位を得ることができる。
Granuleworks
EDEM
Rocky DEM