In this study, we developed a new method for Eulerian hydrocode to represent a solid-solid stick interface in a voxel where multiple solids were included. In this method, individual velocity gradient tensor was assumed for each solid associated with normal vector of the contact area between solids, and determined by solving the equation of force balance on the contact area. We also investigated the accuracy of displacement around the contact area using a layered objects under shearing deformation, and compared the results with those of theoretical solutions and conventional Mixture theory. As the results, in the case that the contact area was on the voxel border, the displacement matched with the theoretical solution as well as the Mixture theory. However, in the case that the contact area was inside voxel, the displacement had still large error, which didn't derive from the individual velocity gradient tensor for each solid.