CFD計算完成后，先輸出為case and data，然后點擊Rocky Export下的Export one-way data＝＞Export current data to Rocky，這是穩態流場與DEM耦合會用到的。此時會在工作目錄中輸出.f2r文件，作為單向耦合時的流場條件輸入給離散元。同時還有各個邊界的幾何的.stl文件。現在可以關閉Fluent，打開Rocky DEM。

         新建項目后右鍵點擊Data瀏覽器中的Geometries，點擊Import Geometry，導bot入上步保存的.cas文件。可以看到各部件以fluent中的邊界命名。選擇其中的upper-inlet、bot-inlet1、channel-inlet、outlet、對稱面sym-sta和對稱面面右鍵點擊刪除，因為這些面在實際物理中是不存在的，可能會發生物質的流進流出。

        圖形界面如下。

        接下來需要創建顆粒的注入入口的幾何與位置。右鍵Geometries，點擊Create Inlet。創建bot_inlet1和upper_inlet兩個顆粒入口，每個顆粒入口屬性如下。

        在Data欄中點擊Physics，設置重力加速度為y-方向9.81m/s2。

        由于模型中有個運動部件——料滾，以52rpm的速度順時針旋轉，首先需定義旋轉運動系，右鍵點擊Motion Frames，選擇Create Motion Frame，然后按下方屬性設置旋轉系屬性，這里的Relative Position就是Fluent里Frame Motion設置的旋轉的origin。Rotation Vector根據右手定則因為料滾在面向我們的這一側是順時針旋轉，所以轉軸指向z-方向，那么vector就是（0,0,-1）。在仿真過程中料滾一直在旋轉，所以Stop Time設置長一些，只要超過你設置的仿真時長就行。Type選擇Rotation，初始速度為-52rpm，而且這里我們認為料滾是勻速轉動，故旋轉加速度為0。

        要想讓料滾運動，還得將這個旋轉運動系賦予到部件上。點擊Geometries下的blade_wall，在Data Editor中點擊Geometry選項卡，在Motion Frame中選擇剛才設置好的旋轉運動系的名稱。

        想預覽部件運動是否正確，回到Motion Frames項目，點擊Preview。此時會彈出一個新的preview的圖形窗口，并且在料滾的軸線處有一個局部坐標系，此時點擊上方的播放按鈕，可以預覽部件播放動作，如果順時針旋轉，則說明剛才的設置無誤。

        接下來點擊右側Windows中的3D Views<01>，退出部件運動預覽模式并回到常規的3D圖形界面。

        創建水稻顆粒、清雜（長短不一的莖稈）和設備壁面材料屬性。右擊Materials，選擇Create Material，分別創建稻種顆粒（rice）、清雜（stem）和設備壁面（equip）材料屬性如下。

        再到Geometries這里，點擊blade_wall，在Material這里選擇equip，給料滾賦予equip材料。同理所有其他設備壁面都這樣進行操作。

        各種材料之間需要設置接觸參數，點擊Material Interactions，坐下角兩個下拉欄可以選擇已設置的任意兩種材料之間的配對，比如這里選擇rice和equip，設置稻種和壁面之間的靜摩擦系數、動摩擦系數和碰撞恢復系數等。同理設置其他材料對的接觸參數如下。

        之后需要設置顆粒種類和形狀。以稻種顆粒為例，右擊Particles，點擊Create Particle，在Daata Editors的Particle中命名為rice，Shape選擇Sphero-Polyhedron（多角球體）。下方的Size中的Size Type選擇篩分尺寸（Sieve Size），點擊+添加尺寸分類，單位設置為m，短軸尺寸設置為0.007（因硬件條件有限顆粒巨多，所以在原尺寸基礎上進行了一定倍數的粗?；?。在Shape標簽中的Verticle Aspect Ratio中填入2。點擊View按鈕在圖形界面查看顆粒形狀。同理設置短、中、長莖稈（分別命名為sstem、mstem、lstem），Shape為直纖維（Straight Fiber），Sieve Size均為0.002，Verticle Aspect Ratio則分別為8、15、25，它們也進行了同樣倍數的粗粒化。

        點擊右側WIndows中的3D View<01>，返回模型界面?，F在要設置每個顆粒注入口的顆粒注入種類、注入流率、注入持續時間等。右擊Data欄中的Inputs，點擊Create Continuous Injection。首先為右上方入料口創建粒子注入參數，命名為upper_inlet，Entity Point設置為upper_inlet，點擊下面particle選項卡的+號添加注入的顆粒種類并填入相對應的流量如下。

        切換到Time選項卡，輸入開始時間為0s，結束時間為1000s。

        同理設置一個名為bot_inject 的Continuous Injection，Entry Point為bot_inlet1，各種顆粒的流量如下，開始和結束時間與upper_inlet相同。

        另外，我們在這里加一個體積填充條件。右擊Inputs，選擇Create Volume Fill，命名為stockgrain，設置Particles和Regoin如下。Region當中的Seed坐標指的是填充顆粒最初發生處的坐標，這個保證在右側腔體內即可。然后幾何中選中所有wall作為限制種子填充的硬邊界。同時在box bounds中設定一個長方體區域將這個顆粒體積填充區域限制其中，需要輸入長方體的中心坐標和各維的長度。相應的區域會在圖形界面中顯示。

        CFD耦合設置：點擊CFD Coupling，下方的Coupling Mode下拉列表中選擇1-Way：Constant（Fluid-Particle），這便是穩態流場單向耦合離散元的方式。此時自動彈出一個文件選擇的對話框，選擇之前由Fluent導出的.f2r文件即可。我們會發現CFD Coupling下多出一個1-Way Fluent，點擊它左下角會出現一個設置粒子受力模型的地方。這里對于非球形顆粒，所有的Drag Law均設置為Schiller&Nauman，Lift Law選擇Saffman。

        隨后點擊Domain Settings，勾選Use Boundary Limits，軟件根據幾何范圍自動確定計算域的范圍。

        點擊Solver，設置Simulation Duration為5.8s，輸出頻率為0.05s輸出一個結果，在General選項卡中可設置CPU或GPU。貌似無需設置仿真的時間步長，這個可能是軟件根據模型情況自動決定的。最后點擊Start按鈕進行計算。

        計算過程中可點擊上方的顆粒相關變量按顏色顯示的下拉列表，查看顆粒當前的各種變量云圖分布。

        根據我的使用狀況，并結合以前使用EDEM的經歷，總結下來：

        1）Rocky在非球形顆粒的形狀設置上比EDEM要方便一些。用過EDEM的都知道，如果是形狀稍微復雜的顆粒可能需要用一個個元顆粒填充，還要用到bond模型，工作量不小。

        2）Rocky允許在計算前預覽部件運動，這一點方便判斷Frame Motion是否設置正確。

        3）畢竟現在同屬于Ansys的模塊，相比EDEM，Rocky和Fluent的耦合更加絲滑。而EDEM與Fluent需要通過UDF耦合，并且有版本限制，配置較為繁瑣。關鍵是這個耦合UDF的更新并不是Ansys的人員執行的，而是Altair的工程師，而且更新得并不勤。當然EDEM和HyperWorksCFD算一對CP，但是我個人非常不愛用HyperWorksCFD，覺得它的內包面功能比較差。

        4）Rocky在計算時可以查看有關顆粒的變量云圖分布，而EDEM基本只能等計算結束。

        5）因為EDEM的市場占有率較高，現在已經開發出了較多的和其他軟件的耦合方法，如AnsysWorkbench、Recurdyn、SimSolid等等；Rocky目前只能和Ansys自己的模塊耦合。

        6）求解效率上感覺EDEM高一點，而且EDEM允許同時使用GPU和多個CPU；Rocky不能兩者同時使用。

        7）EDEM提供的材料庫對于做顆粒物料的對標試驗有極大的便利，它可以根據輸入的休止角提供若干種接觸系數的組合；而Rocky似乎沒看到有這樣的功能。

        8）也許是我Rocky使用還不熟，在后處理上Rocky似乎不太好測量每種顆粒在某處的質量流量，這一點EDEM更方便，EDEM后處理能夠設置質量流量傳感器，還能通過設置Grid來顯示物料不同分層的速度等。