筒倉試驗常用於探討顆粒體的流動特性，包含速度、壓力、流量等，不同之筒倉幾何參數及顆粒材料參數會呈現不同的流動型態。本研究係用PFC 3D軟體並以實驗室內筒倉之幾何參數及顆粒材料參數進行數值模擬，並分析顆粒於筒倉內流動之型態、速度、流量、流動路徑軌跡，以及顆粒流出筒倉後落入下方矩形渠槽之堆積分佈。本模擬將筒倉內顆粒依垂直高度分為五種不同之顏色，依流動型態分析本模擬之結果類似於漏斗狀流動(funnel flow)，顆粒流動時流量峰值與低值會有週期性之變化；經單顆粒追蹤結果發現越靠近筒倉上方之顆粒，在流動時會有往中間集中之現象；顆粒從筒倉落入下方矩形渠槽後堆積之分佈與原本筒倉內顆粒之分佈差異很大。

Granular flows in Silos are frequently examined by both experimental work and numerical simulations. The key physical parameters in silos include the velocity, the pressure, the heap shape and the flow rate. The PFC 3D model is employed in this study and to explore the granular flows in an experimental silo. The simulated flow patterns, velocity, flow rate, and particle trajectories of the particles flow, are compared with experimental results. The flow patters depict he funnel-type flows. The flow rate from the silo shows a fluctuation due to the arching effect near the orifice. The flow paths of top layer are more concentrated them the bottom layer in the silo.

[1] 吳侑謙(2015)，「顆粒特性及水流條件對顆粒體運動及淤積型態之實驗研究」，國立中央大學土木工程學系碩士論文。

[2] 徐松圻、蔣志宏(2009)，「分離元素法應用於土石流行為之研究」，第十三屆大地工程學術研究討論會。

[3] 唐昭榮、胡植慶、羅佳明、林銘朗(2009)，「遽變式山崩之PFC3D模 擬初探－以草嶺與小林村為例」，地工技術第122 期，第143-152頁。

[4] 陳泰安(2001)，「球形顆粒系統之碰撞分析及演算法」，國立中央大學土木工程學系博士論文。

[5] 張光宗、江志文(2012)，「應用分離元素法及連續體力學法分析流槽中砂土的流動與堆積」，水土保持學報第44卷-第三期： 第251-266頁。

[6] 莊明衛(2011)，「球型與非球型顆粒體在平面儲槽內流動行為之研 究」，國立中央大學機械工程學系碩士論文。

[7] 張欽舜(2014)，「顆粒流對垂直平板撞擊力之數值模擬」，國立中央大學土木工程學系碩士論文。

[8] 黃郅軒(2011)，「儲槽內顆粒流動與發聲特性之研究」，國立中央大學土木工程學系碩士論文。

[9] 曾偉鑫(2011)，「二維堆積顆粒崩塌過程之研究」，國立中央大學土木工程學系碩士論文。

[10] 蔡輝平(2013)，「不同含水條件下堆積顆粒塊體崩塌過程之實驗與數值模擬」，國立中央大學土木工程學系碩士論文。

[11] 羅佳明、鄭錦桐、驥樹勇、紀宗吉、林銘朗(2011)，「落石型山崩危害範圍劃設之應用」，中興工程季刊第112期，第23-33頁。

[12] 美國工程工具箱網站 http://www.engineeringtoolbox.com/

[13] Ahmed, I. B., Nsom, B., Allal, A. M., & Ghomari, F. (2010) “Stress distribution during a silo filling or a discharging process,, Nature & Technology, (3), pp.29-36.

[14] Belytschko,T.(1983) “An Overiew of semidiscretization and Time Integration Procedures,, in Computational Methods for Transient Analysis, Ch.1, pp.1-65.

[15] Chau, K. T., Wong, R. H. C., & Wu, J. J. (2002) “Coefficient of restitution and rotational motions of rockfall impacts,, International Journal of Rock Mechanics and Mining Sciences, 39(1), pp.69-77.

[16] González-Montellano, C., Ayuga, F., & Ooi, J. Y. (2011) “ Discrete element modelling of grain flow in a planar silo: influence of simulation parameters,,Granular Matter, 13(2), pp.149-158.

[17] Itasca Consulting Group Inc.(2010) “ PFC3D-Particle Flow Code in three Dimensions,, Minneapolis, USA.

[18] Mao, K., Wang, M. Y., Xu, Z., & Chen, T. (2004) “DEM simulation of particle damping,, Powder Technology, 142(2), pp.154-165.

[19] Masson, S., & Martinez, J. (2000) “Effect of particle mechanical properties on silo flow and stresses from distinct element simulations,, Powder Technology,109(1), pp.164-178.

[20] Ostendorf, M., & Schwedes, J. (2005) “Application of particle image velocimetry for velocity measurements during silo discharge,, Powder Technology, 158(1), pp.69-75.

[21] Perng, A. T. H., Capart, H., & Chou, H. T. (2006)“Granularconfigurations, motions, and correlations in slow uniform flows driven by an inclined conveyor belt,, Granular Matter, 8(1), pp.5-17.

[22] Slominski, C., Niedostatkiewicz, M., & Tejchman, J. (2007) “ Application of particle image velocimetry (PIV) for deformation measurement during granular silo flow,, Powder Technology, 173(1), pp.1-18.

[23] Uñac, R. O., Vidales, A. M., Benegas, O. A., & Ippolito, I. (2012) “ Experimental study of discharge rate fluctuations in a silo with different hopper geometries,, Powder Technology, 225, pp.214-220.

[24] Zhou, G. G., & Sun, Q. C. (2013) “Three-dimensional numerical study on flow regimes of dry granular flows by DEM,, Powder Technology, 239, pp.115-127.