  With the increasing demand for the accuracy of numerical simulation of
pavement mechanics, the variational inequality model and its induced finite
element method which can simulate the interlayer contact state becomes a
potential solution. In this paper, a layer decomposition algorithm for solving
variational inequality models of multi-layer elastic contact systems with
interlayer Tresca friction conditions is studied. Continuous and discrete
versions of the algorithm and their convergence theorems have been proposed and
proved successively. Then, the algebraic form of the executable optimization
algorithm and the numerical experimental results verify the practicability of
the variational inequality model and its algorithm in the pavement mechanics
modeling.
