Design/methodology/approach The research adopts the Weierstrass–Mandelbrot (W-M) function for simulating bidirectional surface profiles. Through the application of elastic-plastic contact theory, it evaluates the contact area and load of a singular asperity across elasticity, elastoplasticity and plasticity phases. The contact load and stiffness of the rough surface are determined using a refined asperity density distribution function, and these findings are juxtaposed with extant models to validate their precision and rationality. Findings The study delineates the influence of fractal dimension (D), surface roughness (G), ellipse eccentricity (e) and friction coefficient (µ) on the contact area, load and stiffness. It reveals that the contact area enlarges with the fractal dimension (D) yet diminishes with increased eccentricity (e), roughness (G) and friction coefficient (µ). These elements considerably affect the contact load and stiffness, underscoring their significance in comprehending surface interactions. Originality/value This study applies fractal theory to analyze the contact mechanics of bidirectional anisotropic surfaces, considering the geometry and mechanics of ellipsoidal asperities on rough surfaces to develop a contact mechanics model. This model clarifies the deformation of an asperity in normal contact, presenting a more rational alternative to current models.