Context. Massive and intermediate-mass stars reside in binary systems much more frequently than low-mass stars. At the same time, binaries containing massive main-sequence (MS) component(s) are often characterised by eccentric orbits, and can therefore be observed as eccentric ellipsoidal variables (EEVs). The orbital phase-dependent tidal potential acting on the components of EEV can induce tidally excited oscillations (TEOs), which can affect the evolution of the binary system. Aims. We investigate how the history of resonances between the eigenmode spectra of the EEV components and the tidal forcing frequencies depends on the initial parameters of the system, limiting our study to MS. Each resonance is a potential source of TEO. We are particularly interested in the total number of resonances, their average rate of occurrence and their distribution in time. Methods. We synthesised 20,000 evolutionary models of the EEVs across the MS using Modules for Experiments in Stellar Astrophysics (MESA) software for stellar structure and evolution. We considered a range of masses of the primary component from 5 to 30 M . Later, using the GYRE stellar non-adiabatic oscillations code, we calculated the eigenfrequencies for each model recorded by MESA. We focused only on the l = 2, m = 0, +2 modes, which are suspected of being dominant TEOs. Knowing the temporal changes in the orbital parameters of simulated EEVs and the changes of the eigenfrequency spectra for both components, we were able to determine so-called 'resonance curves', which describe the overall chance of a resonance occurring and therefore of a TEO occurring. We analysed the resonance curves by constructing basic statistics for them and analysing their morphology using machine learning methods, including the Uniform Manifold Approximation and Projection (UMAP) tool.Results. The EEV resonance curves from our sample are characterised by striking diversity, including the occurrence of exceptionally long resonances or the absence of resonances for long evolutionary times. We found that the total number of resonances encountered by components in the MS phase ranges from ∼10 2 to ∼10 3 , mostly depending on the initial eccentricity. We also noticed that the average rate of resonances is about an order of magnitude higher (∼10 2 Myr −1 ) for the most massive components in the assumed range than for EEVs with intermediate-mass stars (∼10 1 Myr −1 ). The distribution of resonances over time is strongly inhomogeneous and its shape depends mainly on whether the system is able to circularise its orbit before the primary component reaches the terminalage MS (TAMS). Both components may be subject to increased resonance rates as they approach the TAMS. Thanks to the lowdimensional UMAP embeddings performed for the resonance curves, we argue that their morphology changes smoothly across the resulting manifold for different initial EEV conditions. The structure of the embeddings allowed us to explore the whole space of resonance curves in terms of their morphology and to iso...