This study presents numerical analyses for edge chipping by impact loading. As a numerical analysis method, we extend Particle Discretization Scheme Finite Element Method (PDS-FEM) developed by the authors to be able to simulate fracture due to impact loading. We performed simulations targeting edge chipping of soda-lime glass by impact of rigid iron sphere and examined the crack morphology while varying the diameter of the impactor, the impact velocity, and the impact distance. The proposed method successfully simulates the 3D complex crack pattern on edge chipping such as Hertzian cone crack and conchoidal chip scar. The method also reproduces the change of crack morphologies depending on the impact force and the impact distance. Also, a series of numerical analyses are presented to reveal the effect of the impactor geometry on the chip dimensions. The height of chip is independent of the impactor geometry while the width of chip depends on it. According to the agreement with experimental results, it is confirmed that the proposed method is capable of realizing edge chipping due to impact loading.