The collision of cylindrical asperities is simulated by means of a commercial multiphysics software in this paper. The model combines a computational structure mechanics approach with a contact algorithm. Numerical results are first compared to classical Hertzian theory. Effects of the bulk material under the mating surfaces have been included in the simulation. Materials of top and bottom surfaces have been assigned different properties. Development of plastic deformation and friction are presented along the process of collision. Strain hardening is included in the elastic-plastic deformation. The effect of interference change and friction coefficient variation on deformation and friction has been studied. Deformed shape and subsurface residual stress in asperity are output parameters from this study.