自然界中石頭或是砂土顆粒充滿不同形狀，有角狀顆粒，也有球狀顆粒，二者常常混合在自然界或是工程應用上，例如混凝土攪拌中，預拌材料總是混合著角狀及球狀顆粒。前人曾研究利用頂點對面及稜邊對稜邊二種接觸情況發展出三維多面塊體程式及四圓弧逼近法發展三維混合型橢球顆粒堆積程式，本研究將二者所發展的多面體程式與橢球堆積程式進行結合並修改，藉由加入橢球與多面塊體的接觸判斷，使得修改之程式能增加執行混合型剛體運動數值模擬之功能。
本研究加入橢球與多面塊體的接觸判斷方法，以四圓弧逼近法所形成的橢球，可先行判斷球面與弧面何者可能與多面塊體接觸，若是可能與球面接觸，則以球面球心到接觸點距離跟球面半徑大小比較判斷橢球是否入侵多面塊體，反之若是弧面可能與多面體接觸，橢球以弧面球心到接觸點距離與弧面球心半徑大小比較，以此方法分順序與多面體的(1)頂點(2)稜線(3)面做接觸判斷，而稜線的還須判斷所接觸點是否在稜線之中，面的接觸判斷須判斷接觸點是否在面的範圍之中，以此方法結合橢球與多面塊體的接觸碰撞。
完成修改後的程式進行各種不同情況模擬及驗証，其中包含長扁橢球對(1)面、(2)稜線與(3)頂點的碰撞驗証，並進行混合多顆剛體自由落體碰撞測試，確認程式無誤後，應用程式模擬陀螺儀運動及單剪試驗。

The stone or sand particle are full of different forms in the nature. There are angle or spherical particles . The two particles often mix in the nature or apply to engineer .For example in the concrete is mixed , mix materials always mix angle and spherical particles. Forefather who use contact types are categorized by contacts of vertex-to-face and edge-to-edge only write 3-D polyhedral blocks program. And ever develop granular assemblies with mixed ellipsoidal particles program. In this research, we assemblies two programs.Add ellipsoidal and polyhedral contact detection let can simulate mixed 3D ridge body
In the study, we add ellipsoid and polyhedron contact method. we use the ellipsoid with four arcs. First ,we determine which arcs will contact polyhedron. If polyhedron may touch with the sphere ,we determine contact or not which use comparing distance of contact point to sphere center and radius of sphere . If polyhedron may touch with the arc face , we compare distance of contact point to arcs face center and radius of arcs face. We use this kinds of method in order contact detection : (1) vertices (2 ) edges (3) planes , but edges contact must check contact point in the edge line or not. In the planes contact must check contact point in the planes range.
Finish correctional program, we simulate contact include: (1) planes (2 ) edges (3)vertices contact prove . And simulate freely falling body and contact test.After confirming , we use program to simulate simple shear test and gyroscope motion.

1.Cundall, P. A., “Formulation of Three-dimensional Distinct Element Model—PartⅠ. A Scheme to Detect and Represent Contacts in a System Composed of Many Polyhedral Blocks,” Int. J. Rock Mech. Min. Sci. & Geomech. Abstr., Vol. 25, No. 3, pp. 107-116 (1988).
2.Cundall, P.A., “Formulation of a Three-dimensional Distinct Element Model—Part Ⅱ. Mechanical Calculations for Motion and Interaction of a System Composed of Many Polyhedral Blocks.” Int. J. Rock Mech. Min. Sci. & Geomech. Abstr., Vol. 25, No. 3, pp. 117-125 (1988).
3.Cundall, P.A., and O.D.L. Strack, “A Discrete Numerical Model for Granular Assemblies,” Geotechnique, Vol. 29, No. 1, pp. 47-65 (1979).
4.Das, B.M., Principle of Soil Dynamics, PWS-KENT, Boston, USA, pp. 5-51 (1993).
5.Oldenburg, M. and L. Nilsson, “The Position Code Algorithm for Contact Searching,” Int. J. Numer. Meth. Engng., Vol. 37, pp. 359-368 (1994).
6.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).
7.Vaid, Y. P. and J. C. Chern, “Cyclic and Montonic Undrained Response of Saturate Sands,” Advances in the Art of Testing Soils under Cyclic Conditions, pp. 120-147 (1985)
8.Wang, S.P. and E. Nakamachi, “The Inside-Outside Contact Search Algorithm,” Int. J. Numer. Meth. Engng., Vol. 40, pp. 3665-685(1997).
9.Zhong, Z.H. and L. Nilsson, “A Contact Searching Algorithm for General 3-D Contact-Impact Problems,” Computers and Structures, Vol. 34, No. 2, pp.327-335 (1990).
10.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)
11.Ng, T.-T. and X. Lin, “Short Communication Contact Detection Algorithms for Three Dimensional Ellipsoids in Discrete Element Modelling,” International Journal forNumerical and Analytical Methods in Geomechanics, vol. 19, pp. 653-659 (1995)
12.Ng, T.-T. and E. Petrakis, “Small-Strain Response of Spheres Using Discrete Element Method,” Journal of Engineering Mechenics, ASCE, Vol. 122, No. 3, pp. 239-244 (1996)
13.Rothenburg, L. and R. J. Bathurst, “Analytical Study of Induced Anisotropy in Idealied Granular Materials,” Geotechnique, Vol. 39, No. 4, pp. 601-614 (1989).
14.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)
15.Rothenburg, L. and R. J. Bathurst, “Micromechanical Features of Granular Assemblies With Planar Elliptical Particles,” Geotechnique, Vol. 42, No. 1, pp. 79-95 (1992)
16.Lin, J.S., and S. Chen, “Discrete Element Modeling of Direct Shear and Simple Shear Tests,” Proceedings of 2nd Int. Conference on DEM, IESL Publicatons, pp. 245-252 (1993).
17.Wang, C.Y. and V.C., Liang, “A Packing Generation Scheme for the Granular Assemblies with Planar Elliptical Particle,” International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 21, pp. 347-358 (1997).
18.Tran, T. X. and R. B. Nelson, “Analysis of Disjoint Two-Dimensional Particle Assemblies.” Journal of Engineering Mechenics, ASCE, Vol. 122, No. 12, pp. 1139-1148 (1996)
19.王泉博，「二維圓磐顆粒之堆積模擬與應用」，碩士論文，國立中央大學土木工程系，中壢(1996)
20.李宗龍，「橢圓顆粒介質體運動之數值模擬」，碩士論文，國立中央大學土木工程系，中壢(1997)。
21.李承恩，「三維球形顆粒介質體變形分析」，碩士論文，國立中央大學土木工程系，中壢(1997)。
22.陳俊全，「混合型橢球顆粒介質體運動之數值模擬」，碩士論文，國立中央大學土木工程系，中壢(1999)。
23.許秀真，「多面塊體系統之運動模擬」，碩士論文，國立中央大學土木工程系，中壢(2001)。
24.鄭凱文，「三維顆粒介質與變形體互制行為之數值模擬」，碩士論文，國立中央大學土木工程系，中壢(2003)。
25.李世昌，「三維多面塊體系統之數值模擬」，碩士論文，國立中央大學土木工程系，中壢(2004)。