In multibody dynamic analysis including contact and impact, there are two major analysis methods, i.e., piecewise analysis and continuous analysis. Modeling of contact phenomena is mainly classified with a Kelvin-Voigt model or a model of Hertz contact model. In this paper, a contact module fur AutoDyn7 program was developed and implemented. Both the Kelvin-Voigt model and a model of Hertz contact law were developed. The process of this module is composed of contact distinction and the contact force calculation. Two examples were verified and compared to the commercial program DADS.