This paper presents an improved mixed finite element method of analyzing elastic contact problems, in which a new numerical algorithm is developed to calculate the flexibility matrix of the contact region. With this method the contact problem of two teeth-engaging is initially solved as a three-dimensional elastic finite element analgsis, In comparison with mixed fimte element method, the improved mixed finite element method is much more efficiency in solving contact problems.