台灣地質條件多斷層、褶皺、破裂面，當颱風、地震等災害發生時，常常伴隨著邊坡的落石掉落。因此能夠事先知道落石掉落的影響範圍是值得探討的課題。相關邊坡落石數值模擬之研究已經有不少學者提出。在眾多數值模擬方法中，不連續體變形分析法能夠模擬塊體的大變位及塊體間接觸之問題，是相當適合模擬落石的數值方法之一。
本研究以Microsoft Visual C++搭配OpenGL繪圖函式庫撰寫三維顆粒狀不連續變形分析法之程式碼。本研究將落石假設為球狀，邊坡由一連串三角形組合而成，坡上的樹木假設為圓柱狀，深入探討球與球以及球與邊界元素(地形、圓柱)的碰撞行為，並透過自行假設之虛擬案例來驗證程式的正確性。

Because the geology of Taiwan is full of faults, folds, and fractures, rockfalls on slope may occur due to typhoons or earthquakes. Therefore, it is worth noting that forecasting the zone of rockfalls. There have been a lot of numerical rockfall simulation researches proposed. Compared with other numerical simulation methods, discontinuous deformation analysis (DDA) can stimulate the large deformation of blocks and the contact problems between blocks. DDA is one of numerical methods that is suitable to simulate rockfalls.

This study creates a program based on three-dimensional particle discontinuous deformation analysis using Microsoft Visual C++ and OpenGL. In this study, rockfalls are assumed to be spheres, slope consists of a series of triangles, and trees on slope are considered cylinders. The contact mechanisms of both sphere-sphere and sphere-boundaries (trees or terrain) have been discussed. The validity of the program is verified through several fictitious cases assumed by ourselves.

1. 石根華(1997)。數值流形方法與不連續體變形分析(裴覺民，譯)。清華大學出版社，北京，中國。(原著第一版出版於1997年)
2. 李錦育(2010)。山坡地社區落石潛勢危害分析。黃埔學報，58，pp. 15-24。
3. 陳春華(2008)。不連續體變形分析法模擬集集地震引致草嶺邊坡崩塌模擬之研究。成功大學土木工程學系學位論文，台南，台灣。
4. 陳修、顧承宇(1995)。不連續變形分析法於岩石工程上之研究與應用。 財團法人中興工程顧問社，台北，台灣，pp. 6-8。
5. 馮康(1965)。基於變分原理的差分格式。應用數學與計算數學，2(4)，pp. 238~262。
6. 葛德治(1997)。The Issue of Rock Deformability in the Stability Analyses of Rock Slope。第七屆大地工程學術研究討論會，金山，台灣，pp. 1259~1266。
7. 葛德治 (1996)。Development of Rockfalling Simulation Program - ROCKPATH I。1996 岩盤工程研討會論文集，新竹，台灣，pp. 355-364。
8. 蕭沛佳(2005)。不連續變形分析法應用於礫石型土石流運動之研究。中興大學水土保持學系學位論文，台中，台灣。
9. 羅佳明、鄭錦桐、冀樹勇、紀宗吉、林銘郎(2011) 。落石型山崩危害範圍劃設之應用。中興工程季刊，112，pp. 23-33。
10. 顧承宇、翁孟嘉、高憲彰、陳建忠、李怡先(2006)。三維雷射掃描技術於岩坡落石分析之應用。2006岩盤工程研討會論文集，台南，台灣，pp. 387-396。
11. 顧承宇、蔡勇賢、蘇群川(2009)。非連續體模式應用於土石流問題之分析。，第十三屆大地工程學術研究討論會，宜蘭，台灣。
12. Beyabanaki, S. Amir Reza, Mikola, Roozbeh Geraili, Hatami, Kianoosh. (2008). Three-dimensional Discontinuous Deformation Analysis (3-D DDA) using a New Contact Resolution Algorithm. Computers and Geotechnics, 35(3), pp. 346-356.
13. Beyabanaki, S. Amir Reza., Bagtzoglou, Amvrossios C. (2012). Three-Dimensional Discontinuous Deformation Analysis (3-D DDA) Method for Particulate Media Applications. Geomechanics and Geoengineering, 7(4), pp. 239-253.
14. Beyabanaki, S. Amir Reza, Bagtzoglou, Amvrossios C. (2015). Sphere-Boundary Edge and Sphere-Boundary Corner Contacts Model in DDA for Simulating Particulate Media in 3-D. Geomechanics and Geoengineering, 10(2), pp. 83-94.
15. Cundall, Peter A. (1971). A Computer Model for Simulating Progressive, Large Scale Movements in Blocky Rock Systems. Proceedings Symposium, International Society of Rock Mechanics, Nancy 2,7.
16. Cundall, Peter A., Hart, Roger D. (1989). Numerical Modelling of Discontinua, in Proceedings of the 1st U.S. Conference on Discrete Element Methods (Golden, Colorado, October 1989), pp. 1-17. G.G.W. Mustoe, M. Henriksen and H-P. Huttelmaier, Eds. Golden, Colorado: CSM Press.
17. Desai, C. S., Zaman, M. M., Lightner, J. G., Siriwardane, H. J. (1984). Thin-layer Element for Interplane and Joints, International Journal for Numerical and Analytical Methods in Geomechanics, 8, pp. 19- 43.
18. Doolin, D. M., Sitar, N. (2004). Time integration in discontinuous deformation analysis. Journal of Engineering Mechanics, 130, pp. 249–258.
19. Goodman, Richard E., Taylor, Robert, L., Brekke, Tor L. (1968). A Model for Mechanics of Jointed Rock, Journal of the Soil Mechanics and Foundations Division, A. S.C. E., 94, pp. 637- 659.
20. Guzzetti, F., Crosta, G., Detti, R., and Agliardi, F. (2002). STONE: A Computer Program for the Three-Dimensional Simulation of Rock-falls. Computer and Geosciences, 28(9), pp. 1081-1095.
21. Hinton, E, Owen, DRJ. (1979). An Introduction to Finite Element Computations. Pineridge Press Limited, Swansea, UK.
22. Itasca Consulting Group Inc. (2002). PFC3D (Particle Flow Code in Three Dimensions) Theory and Background, Version 3.0, Minneapolis, USA, pp. 1-26.
23. Jiang, Qing-Hui, Feng, Ding-Xiang. (2000). Study on an Angle-Plane Model of Three-Dimensional Discontinuous Deformation Analysis. Chinese Journal of Rock Mechanics and Engineering, 19, pp. 930-935.
24. Jiang, Qing-Hui, Yeung, MR. (2004). A Model of Point-to-Plane Contact for Three Dimensional Discontinuous Deformation analysis. Rock Mechanics and Rock Engineering, 37, pp. 95- 116.
25. Jiao, Yu-Yong, Huang, Gang-Hai1, Zhao, Zhi-Ye, Zheng, Fei1, Wang, Long. (2015). An Improved Three-Dimensional Spherical DDA Model for Simulating Rock Failure. Science China Technological Sciences, 58(9), pp. 1533-1541.
26. Jing, Lanru. (1998). Formulation of Discontinuous Deformation Analysis (DDA) – An Implicit Discrete Element Model for Block Systems. Engineering Geology, 49(3), pp. 371-381.
27. Liu, Jun, Kong, Xian-Jing, Lin, Gao. (2004). Formulations of the Three-Dimensional Discontinuous Deformation Analysis Method. Acta Mechanica Sinica, 20(3), pp. 270-282.
28. Ke, Te-Chih, Bray, Jonathan. (1995). Modeling of Particulate Media using Discontinuous Deformation Analysis, Journal of Engineering Mechanics, 121(11), pp. 1234-1243.
29. Ku, Cheng-Yu, Chern, Jin-Ching. (1998). Modification of the DDA Method for Rigid Block Problems, International Journal of Rock Mechanics and Mining Science, 35(6), pp. 683-693.
30. Liu, Jun, Kong, Xian-Jing. (2002). Three Dimensional Continuous and Discontinuous Deformation Analysis. ACTA Mech Sinica, 34, pp.941-94 8 .
31. Munjiza, A., Andrews, K. R. F. (1998). NBS Contact Detection Algorithm for Bodies of Similar size. International Journal for Numerical Methods in Engineering, 43(1), pp. 131-149.
32. Perkins, Eric, William, John R. (2001). A Fast Contact Detection Algorithm Insensitive to Object sizes. Engineering Computations, 18(1-2), pp. 48-62.
33. Piteau, Douglas R., Clayton, R. (1977). Discussion of Paper on Computerized Design of Rock Slope Using Interactive Graphics for the Input and Output of Geometrical Data by P. A. Cundall, M. D. Voegele
and C. Farihurst. Proceedings of the 16th Symposium on Rock Mechanics, Univ. of Minnesota, Minneapolis, pp. 62-63.
34. Shi, Gen‐Hua, Goodman, Richard E. (1984). Discontinuous Deformation Analysis, Proceedings of the 25th U. S. Symposium on Rock Mechanics, Evanston, Illinois, USA, pp. 269-277.
35. Shi, Gen‐Hua. (1988). Discontinuous Deformation Analysis: A New Numerical Model for the Statics and Dynamics of Block Systems. University of California, Berkeley, USA.
36. Shi, Gen‐Hua. (1989), Block System Modeling by Discontinuous Deformation Analysis, Ph.D. Thesis, University of California, Berkeley, USA.
37. Shi, Gen‐Hua. (2001). Three-Dimensional Discontinuous Deformation Analysis. In: Bicenic N, Editor. Proceedings of the Forth International Conference on Analysis of Discontinuous Deformation (ICADD-4), Glasgow, Scotland, June 6–8, pp. 1–21.
38. Shi, Gen‐Hua. (2001). Three-Dimensional Discontinuous Deformation Analysis. In DC Rocks 2001, The 38th US Symposium on Rock Mechanics (USRMS). American Rock Mechanics Association, Washington, D.C.,USA, pp. 1421-1428.
39. Belytschko, Ted. (1983). An Overview of Semidiscretization and Time Integration Procedures. in Computational Methods for Transient Analysis, North-Holand, Ch. 1, pp. 1-65.
40. Wang, Guo-Rong, Wei, Yi-Min, Qiao, San-Zheng. (2004). Equation Solving Generalized Inverses. Generalized Inverses: Theory and Computations. Beijing: Science Press. pp.6.
41. Wang, Xiao-Bo, Xu, Wen-Jie, Zhang, Bing-Yin, Sun, Qi-Cheng. (2014). Particle Crushing Simulations with Improved Discontinuous Deformation Analysis. Engineering Computations, 31(7), pp. 1321-1341.
42. Williams, John R., Hocking, G., Mustoe, Graham G.W. (1985). The Theoretical Basis of the Discrete Element Method. NUMETA, Numerical Methods of Engineering, Theory and Applications Rotterdam, A.A. Balkema, Rotterdam, The Northland.
43. Wu, Jian‐Hong, Juang, C. Hsein, Lin, Hung‐Ming. (2005). Vertex-to-Plane Contact Searching Algorithm for Three-Dimensional Frictionless Contact Problems. International Journal for Numerical Methods in Engineering, 63(6) , pp. 876-897.
44. Wu, Jian‐Hong, Ohnishi, Yuzo, Nishiyama, Satoshi. (2005). A Development of the Discontinuous Deformation Analysis for Rock Fall Analysis. International Journal for Numerical and Analytical Methods in Geomechanics, 29(10), pp. 971-988.
45. Wu, Jian‐Hong. (2010). Compatible Algorithm for Integrations on a Block Domain of Any Shape for Three-Dimensional Discontinuous Deformation Analysis. Computers and Geotechnics, 37(1), pp. 153-163.
46. Yeung, MR, Jiang, Q.-H., Sun N. (2003). Validation of Block Theory and Three Dimensional Discontinuous Deformation Analysis as Wedge Stability Analysis Methods. International Journal of Rock Mechanics and Mining Sciences, 40(2), pp. 265-275.
47. Zhang, Hong, Chen, Guang-Qi, Zheng, Lu, Han, Zheng, Zhang, Ying-Bin, Wu, Yan-Qiang, Liu, Shu-Guang. (2015). Detection of Contacts between Three-Dimensional Polyhedral Blocks for Discontinuous Deformation Analysis. International Journal of Rock Mechanics and Mining Sciences, 78, pp. 57-73.