A FEA model of two-dimensional (2D) sliding between two interfering fractal rough solids was build, and the results was presented. The trends in shear stress between the rough solids were provided by the FEA. Combined with the delaminating theory of wear, the Maximum shear stresses and its distance are presented when loading and sliding. The result showed that the Maximum shear stresses was located in the subsurface which range about 2 to 8 um from the surface. In addition, with different plastic characteristics, the location of Maximum shear stresses were different. With modeling the contact and sliding between rough surfaces, the friction, sub-surface crack initiation and propagation and wear mechanisms can be understood deeply.