Aimed at the problem that the calculation results of existing elasto-plastic lubrication models are distorted under the working condition of high viscosity lubricating oil,a plasto-elastohydrodynamic lubrication (PEHL) model was established based on a semi-analytical method.The numerical simulation was successfully carried out through a stepping technique in increasing lubricating oil viscosity,which eliminated the previous fluctuation in the pressure profile in the contact,and overcame the previous problem of calculation distortion.Plastic strain was solved based on J2 flow theory under isotropic strengthening criterion and by radial return algorithm.The residual deformation caused by plastic strain was taken into account in the oil film thickness equation.The results show that the main peak of PEHL pressure is significantly lower than that of EHL,and the pressure spike changes little.The central film thickness for PEHL and EHL models are basically consistent,while the minimum film thickness of the PEHL model is a little higher.The increases in oil viscosity and entraining speed can move the plastic strain region closer to the outlet region and the surface.