Photonic bound states in the continuum (BICs) are spatially localized modes with infinitely long lifetimes, which exist within a radiation continuum at discrete energy levels. These states have been explored in various systems, including photonic and phononic crystal slabs, metasurfaces, waveguides, and integrated circuits. Robustness and availability of the BICs are important aspects for fully taming the BICs toward practical applications. Here, we propose a generic mechanism to realize BICs that exist by first principles free of fine parameter tuning based on non-Maxwellian double-net metamaterials (DNMs). An ideal warm hydrodynamic double plasma (HDP) fluid model provides a homogenized description of DNMs and explains the robustness of the BICs found herein. In the HDP model, these are standing wave formations made of electron acoustic waves (EAWs), which are pure charge oscillations with vanishing electromagnetic fields. EAW BICs have various advantages, such as (i) frequency-comb-like collection of BICs free from normal resonances; (ii) robustness to symmetry-breaking perturbations and formation of quasi-BICs with an ultrahigh Q-factor even if subject to disorder; and (iii) giving rise to subwavelength microcavity resonators hosting quasi-BIC modes with an ultrahigh Q-factor.