While it is important to design stellarators with high magnetohydrodynamic stability β-limit, it is also crucial to ensure that good magnetic surfaces exist in a large range of β values. As β increases, pressure-driven currents perturb the vacuum magnetic field and often lead to the emergence of magnetic field line chaos, which can worsen the confinement and is the cause of another kind of β-limit, the so-called equilibrium β-limit. In this paper, we explore numerically the dependence of the equilibrium β-limit on the bootstrap current strength in a classical stellarator geometry using the stepped pressure equilibrium code. We develop a diagnostic to determine whether or not magnetic islands are expected to participate significantly to radial transport, and we build an analytical model to predict the expected equilibrium β-limit, which recovers the main features of the numerical results. This research opens the possibility to include additional targets in stellarator optimization functions, provides additional understanding on the existence of magnetic surfaces at large β, and is a step forward in the understanding of the equilibrium β-limit.