本篇續《Abaqus無網格法之DEM分析案例二則，沙漏與高爾頓板，附仿真源文件》一文。

DEM分析常用于顆?；旌?、篩選等物理過程，目前Abaqus GUI還不支持DEM顆粒建模，不過可以通過編輯關鍵字*particle generator或者運行Python腳本來實現，下面詳細地介紹通過Python腳本生成DEM顆粒單元的方法。

顆粒混合：

顆粒篩選：

首先將DEM部件按照實體建模，并劃分為C3D8R六面體單元，然后生成名為Galton_Board的inp文件。

DEM分析前處理過程：

打開Abaqus Command窗口，按照下面的格式運行solidtodem.py文件，藍框為工作路徑，紅框為剛才生成inp文件名稱，確保solidtodem.py文件與inp文件都在工作路徑內。

運行腳本：

運行完畢后，會生成一個dem_Galton_Board.inp文件，用其中的離散粒子單元替換Galton_Board.inp中的C3D8R單元，并保存。

替換單元：

這樣就完成了從C3D8R到離散粒子單元PD3D的轉化，再基于最新的Galton_Board.inp文件進行修改，定義一下顆粒密度、大小、阻尼與接觸等即可進行高爾頓板的DEM分析，詳細關鍵詞見上篇文章中提供的inp文件。

DEM分析也常與多體分析、流體分析等過程進行耦合，以便計算大量離散粒子對機構、流場的影響。

挖勺力計算：

DEM與流體耦合：

solidtodem.py腳本鏈接：

https://pan.baidu.com/s/1t4ys7bSBvfRat-mOqtJfMg     

提取碼：c9tm