By numerical simulations and experiments of fully chaotic billiard lasers, we show that single-mode lasing states are stable, whereas multi-mode lasing states are unstable when the size of the billiard is much larger than the wavelength and the external pumping power is sufficiently large. On the other hand, for integrable billiard lasers, it is shown that multi-mode lasing states are stable, whereas single-mode lasing states are unstable. These phenomena arise from the combination of two different nonlinear effects of mode-interaction due to the active lasing medium and deformation of the billiard shape. Investigations of billiard lasers with various shapes revealed that single-mode lasing is a universal phenomenon for fully chaotic billiard lasers.