プログラムの大枠

DEMの計算の流れ

①初期設定
　↓
②各要素の合力を計算
　↓
③合力から加速度を計算
　↓
④加速度から速度を計算
　↓
⑤座標更新 →②へ
（速度から微小時間後の座標を計算）

プログラムの流れ

＜初期設定＞
・ステップ時間の設定
・範囲の設定
・Line要素の追加
・Particle要素の追加
・Interfaceの設定
・セットアップ

＜ステップ毎の処理＞
・合力の計算
・速度、加速度、座標の計算

例

dem = DEM()
dem.dt = 0.01
dem.area = [0,300,0,200]
for line in lines:
    dem.lines.append(line)
for p in particles:
    dem.particles.append(p)
dem.setInterface(inf1,1,1)
dem.setInterface(inf2,1,2)
dem.setup()

while dem.step<1000:
   dem.calc_force()
   dem.update_coord()
   dem.step += 1