The paper is concerned with the analysis and the numerical solution of the multimaterial topology optimization problems for bodies in unilateral contact. The contact phenomenon with Tresca friction is governed by the elliptic variational inequality. The structural optimization problem consists in finding such topology of the domain occupied by the body that the normal contact stress along the boundary of the body is minimized. The original cost functional is regularized using the multiphase volume constrained Ginzburg-Landau energy functional. The first order necessary optimality condition is formulated. The optimal topology is obtained as the steady state of the phase transition governed by the generalized Allen-Cahn equation. The optimization problem is solved numerically using the operator splitting approach combined with the projection gradient method. Numerical examples are provided and discussed.