|本論文主要運用軟球模式的三維離散元素法（DEM）來模擬顆粒體在振動床和剪力設備中的運動，探討主題包含速度場分佈、迴流運動、擴散運動及各種與顆粒混合有關的流動性質。 本文首先系統性地提出六項簡單碰撞系統的測試，以動量守恆定律和動能守恆定律來驗證DEM程式的準確性。同時，再模擬振動床中的顆粒運動，其中明顯的對稱迴流現象與相關文獻的實驗結果極為吻合。其實，因為邊壁摩擦作用所產生的顆粒迴流運動，是振動顆粒床中非常重要的運動機制，所以本文詳細探討振動條件（包含振動加速度及振動頻率）對於對稱迴流形成的影響。文中為了有效描述顆粒對流運動和擴散運動的特性，分別定義無因次化對流率Jconv和垂直方向自我擴散係數Dyy，並且以培利數Pe來作為這兩種效應的比值。培利數在對稱迴流的形成過程中，扮演非常重要而有趣的角色，因此文中有詳盡的討論。另外，為了探討顆粒的混合行為，兩群等量而顏色不同的的玻璃珠以上下分置的初始擺放方式，放入振動床中。同時，為了衡量不同顆粒的混合狀態，本文使用Lacey index來作為混合度M以量化混合狀態，而且還可以利用最小平方法的線性擬合，從M隨時間的變化歷程中獲知混合速率。模擬的結果顯示，混合速率隨Dyy的加大而變快，且呈現指數關係的增快趨勢。另外一方面，為了衡量靜電力對顆粒流動的效應，本文定義靜電力和顆粒重量的比值為靜電數Es。模擬的結果發現，粒子溫度隨靜電數的加大而線性上升，而且混合速率常數是以冪次法則隨靜電數的加大而增快。在顆粒剪力流的模擬中，顏色不同的玻璃珠隨機堆積且上下分置，下邊壁等速移動後就會造成顆粒的混合。對於混合層厚度隨時間的變化歷程，DEM的模擬結果與利用擴散方程式的計算結果頗為吻合，這亦表示擴散運動是剪力流內顆粒混合的重要機制。 This thesis examines the mixing behaviors of granular materials subjected to external vertical vibration and sheared force. Three-dimensional discrete element computer simulation is used to study the velocity distribution, convective flow, diffusive motion and granular mixing in both vibrated granular bed and sheared granular flow. A series of systematic validation tests for the DEM program, including six tests of simple collision system and a test of macroscopic phenomenon. The conservation laws of momentum and kinetic energy are used to verify tests of collisions between particles or between particles and boundaries. Also, the simulation result of symmetric convection rolls are compared with the experimental result. With frictional sidewalls, the convection flow is a very important phenomenon in the vibrated granular bed. The influence of vibrating conditions, including vibration acceleration and frequency, on the formation of symmetric convection flow is investigated in this work. In order to characterize the convective flow and the diffusive motion of granular materials, the dimensionless convection flow rate, Jconv, and the vertical self-diffusion coefficient, Dyy, are defined, respectively. The Péclet number, Pe, is employed to characterize the ratio of the convective flow to the diffusive motion in vertical direction. The role of Pe in the formation of symmetric convection flow is discussed in detail. Moreover, the top-bottom initial loading pattern of two groups of glass beads with different colors is employed to investigate mixing behavior of granular materials. The well-known Lacey index is employed as the mixing degree, M, to quantify the mixing quality. The mixing rate is calculated from a least-square fit using the time evolution of M. The simulation results demonstrate that the mixing rate increases with increasing Dyy in exponential relation. In order to characterize the effect of electrostatic force on the granular flow, the Electrostatic Number Es is defined as the ratio of the electrostatic force to the particle weight. The simulation results demonstrate that the granular temperatures increase linearly with the increasing Es number. Meanwhile, the mixing rate constants increase with the increasing Es number in power law relations. In the simulation of sheared granular flow, the initial loading of identical glass beads with different colors is also arranged in a top-bottom loading pattern. The transverse mixing of particles is caused by the moving bottom wall with a constant velocity along the x direction. The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.