The surface roughness induced by geometric irregularities (asperities) has substantial influence on the contact stiffness, which further affects the working performance and service life of mechanical equipment. In this study, an elastic–plastic contact law between a sinusoidal asperity and a rigid smooth flat is first studied, which is then applied on a statistical model to simulate the contact behavior of a pair of 18CrMo4 steel surfaces to investigate the influences of morphology parameters on the contact stiffness. The analysis shows that smaller shape ratios [Formula: see text] and larger wavelengths [Formula: see text] induce larger normal contact stiffness [Formula: see text] for surfaces with identical roughness, wherein the roughness is defined by the mean value of asperity heights [Formula: see text] and the standard deviation of asperity heights [Formula: see text]. The normal contact stiffness increases as [Formula: see text] decreases under the same loading conditions, while the normal contact stiffness increases as [Formula: see text] decreases for surfaces with a fixed [Formula: see text]. Besides, the normal pressure and normal contact stiffness derived from the proposed contact model are validated. The results demonstrate the potential of the proposed model in contact design due to its ability of establishing the relations between the normal contact stiffness and surface morphology parameters.