Taiwan has a special geographical position, and is often threatened by typhoon and earthquake, etc. Hence, the purpose of this research was provided a method for simulated discrete element method (DEM) and vector form intrinsic finite element (VFIFE) and also provided contact numerical analysis procedure with particles and frame element of VFIFE, use DEM as foundation for develop 3D contact theorem was provided by Cundall and include forward and tangential contact force model, and consider five contact types: particle to particle, particle to boundary, particle to frame element, frame element to frame element, and frame element to boundary and verified different contact cases.