於程式中可在有限空間生成不同形狀之顆粒，霣降、堆積後進行安息角試驗模擬，以 MATLAB 7.1 軟體作為後處理，可展現出顆粒在空間的運動狀況、分佈情況與剖面圖，以利分析。
孔隙比、顆粒間之接觸點數、顆粒形狀、顆粒基本摩擦角等堆積體有關參數，為影響顆粒狀材料力學行為之重要因素。本論文針對不同顆粒形狀與安息角的關係做研究，安息角試驗方式分為移除單一邊牆與單一邊牆斜傾90度兩種，依試驗分析結果得知，顆粒形狀確實會影響顆粒堆積體安息角之大小，且扁橢球的安息角最大，卵球次之，長橢球最小。

The domain to be filled with particles can be a polyhedron of any shape. Simulation results can be processed by MATLAB 7.1 to display the distribution of particles in space and the fabric structure at any instant. Parameter studies of the influence of the void ration, the contact numbers among the gravels, particle shape and the internal friction angle between particles on the macroscopic behaviors of particles assemblies are conducted in this thesis. In this study, it emphasizes the relationship between the angle of repose and the particle shape. The reactions of angle of repose are taken account by the relative particle shape through remove one side of wall and set a wall spinning slowly into 90 degree. It is found that the particle shape will influence the angle of repose, ellipsoidal particles are the most, and then egg-shaped particles and the long ellipsoidal particles are the least.

1. 王泉博，「二維圓盤顆粒之堆積模擬與應用」，碩士論文，國立中央大學土木工程研究所，中壢(1996)。
2. 李宗龍，「橢圓顆粒介質體運動之數值模擬」，碩士論文，國立中央大學土木工程研究所，中壢(1997)。
3. 李承恩，「三維球形顆介質體變形分析」，碩士論文，國立中央大學土木工程研究所，中壢(1997)。
4. 劉建宏，「橢球顆粒介質體運動之數值模擬」，碩士論文，國立中央大學土木工程研究所，中壢(1998)。
5. 陳俊全，「混合型橢球顆粒介質體運動之數據模擬」，碩士論文，國立中央大學土木工程研究所，中壢(1999)。
6. 林宏冠，褚炳麟，「複合圓形顆粒元素之應用與建立」，碩士論文，中興大學土木工程研究所，台中(1999)。
7. 陳泰安、王仲宇、盛若磐，「運動顆粒接觸判斷快速搜尋法」，電子計算機於土木水利工程應用研討會88 年論文集，台中，第 503-512 頁(2000)。
8. 許秀真，「剛性多面體系統之數值模擬」，碩士論文，國立中央大學土木工程研究所，中壢(2000)。
9. 郭殷孝，褚炳麟，「複合圓形顆粒落門與安息角試驗之數值模擬」，碩士論文，中興大學土木工程研究所，台中(2000)。
10. 林世洪，「以二維類多邊形顆粒堆積體探討粒狀材料行為」，博士論文，台灣科技大學營建工程研究所，台北，第1～69頁(2002)。
11. 李世昌，「三維多面塊體系統之數值模擬」，碩士論文，中央大學土木研究所，中壢(2003)。
12. 張義隆，「以複合圓形顆粒模式探討粒狀材料之力學行為」，博士論文，中興大學土木工程學研究所，台中， (2003)。
13. 林佳宏，「混合型剛體顆粒系統之數值模擬」，碩士論文，國立中央大學土木工程研究所，中壢(2006)。
14. Brewer, R., Fabric and Mineral Analysis of Soils, John Wiley and Sons, Inc., pp. 129-158 (1964).
15. Cundall, P. A., “A Computer Model of Simulating Progressive Large
Scale Movement in Blocky Rock System,” Symposium of International Society of Rock Mechanics, Nancy, France, pp. 11-19(1971).
16. Cundall, P. A. and O. D. L. Strack, “A Discrete Numerical Model for Granular Assemblies,” Geotechnique, Vol. 29, No. 1, pp. 47-65 (1979).
17. Cundall, P. A. and O. D. L. Strack, “Modeling of Microscopic
Mechanisms in Granular Material.” Mechanics of Granular Materials: Mew Models and Constitutive Relations, Amsterdam, Netherlands, pp. 137-149 (1983).
18. Hakuno, M. and K. Meguro, “Simulation of the Collapse of Concrete Frames and Volcanic Eruption,” Advances in Micromechanics of Granular Materials, Elsevier Science Publishers B.V., pp.321-330(1992).
19. Lin, J.-S. and S. Chen, “Discrete Element Modeling Of Direct Shear and Simple Shear Tests,” Proceedings of the 2nd Int. Conference on DEM, IESL Publications, pp. 245-252 (1993).
20. Lin, X. and T.-T. Ng, “A three-dimensional discrete element model using arrays of ellipsoids,” Geotechnique, Vol. 47, No. 2, pp. 319-329 (1997).
21. Ng, T.-T. and X. Lin, “ Numerical Simulations of Naturally Deposited Granular Soil with Ellipsoidal Elements,” Proceedings of the 2nd Int. Conference on DEM, IESL Publications, pp. 557-567 (1993).
22. Ng, T.-T. and X. Lin, “Short Communication Contact Detection Algorithms for Three Dimensional Ellipsoids in Discrete element Modelling, ” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 19, pp. 653-659 (1995).
23. Ng, T.-T. and E. Petrakis, “Small-Strain Response of Random Arrays
of Spheres Using Discrete Element Method,” Journal of Engineering Mechanics, ASCE, Vol. 122, No. 3, pp. 239-244 (1996).
24. Oda, M., “Initial Fabrics and Their Relations to Mechanical Properties of Granular Material,” Soils and Foundations, Vol. 12, No. 1, pp.17-36 (1972).
25. Rothenburg, L. and R. J. Bathurst, “Analytical study of induced anisotropy in idealized granular materials,” Geotechnique, Vol. 39,No. 4, pp. 601-614 (1989).
26. Rothenburg, L. and R. J. Bathurst, “Numerical Simulation of Idealized Granular Assemblies With Plane Elliptical Particles,” Computers and Geotechnics, Vol. 11, pp. 315-329 (1991).
27. Shi, G.-H., “Discontinuous Deformation Analysis, A new Numerical Model for the Static and Dynamics of Block systems,” Ph. D. Thesis, Department of Civil Engineering, University of California at Berkeley(1988).
28. Tran, T. X. and R. B. Nelson, “Analysis of Disjoint Two-Dimensional
Particle Assemblies.” Journal of Engineering Mechanics, ASCE, Vol. 122, No. 12, pp. 1139-1148 (1996).
29. Wang, C.-Y. and V.-C. Liang, “A Packing Generation Scheme for the Granular Assemblies With Planar Elliptical Particles,” International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 29, pp. 347-358 (1997).
30. Vaid, Y. P. and J. C. Chern,“Cyclic and Monotonic Undrained Response of Saturate Sands,”Advances in the Art of Testing Soils under Cyclic Conditions, pp. 120-147 (1985).