The operating characteristics of machine parts and units are defined in many respects by the physical and mechanical properties of their surface layers. Unfortunately, it is still not quite clear what parameters and mechanisms are responsible for one or another modification of surface layer properties. In this context, computer modeling techniques can be a useful tool in studying the variation of surface properties in contact interaction and run-in. Of fundamental importance is the possibility to consider the processes occurring on nanoscales and scales of individual atoms. In the work, loading conditions in plastic surface deformation was reproduced on the macroscale (traditional approach), atomic scale, and mesoscale by computer modeling with the finite element method, movable cellular automata method, and molecular dynamics method. The modeling results are in good qualitative agreement with data of experimental measurements.