A hypoplastic constitutive model is implemented in a smoothed particle hydrodynamics code for the first time. An improved wall boundary treatment is presented for better performance. The proposed approach is first validated by comparing the numerical results with the analytical solutions for oedometer test. Two more problems involving large deformation, i.e., the collapse of sand between two parallel plates and the failure of a homogeneous slope, are analyzed to demonstrate the performance of the proposed method
