プログラムの大枠

Pythonで個別要素法プログラムの作成



プログラムの大枠

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
ページトップへ