Higher order DDA method for two-dimensional problems has been developed in both theory and computer codes by researchers like Koo and Chern, Ma et al. and Hsiung. Additionally, The DDA contact model which was originally based on penalty method was improved by adopting the Lagrange type approach reported by Lin et al.