In the modeling of elastohydrodynamic lubrication problems considering mixed friction, strongly coupled dependencies occur due to piezo-viscous effects and asperities, which can make a numerical solution exceptionally difficult. A fully implicit coupled scheme for solving mixed elastohydrodynamic lubrication problems is presented. Our scheme uses finite-volume discretization and co-allocated grids for hydrodynamic pressure and elastic deformation. To provide strong coupling between pressure and deformation even in the highly loaded zone, a correction term that adds numerical diffusion is used. The resulting linear equation system of this scheme can be efficiently solved by Krylov subspace methods. This results in an improved accuracy and computational efficiency compared to the existing methods. This approach was validated and has been shown to be accurate.