We give an overview of results for the quenched glueball spectrum from two-body bound state equations based on the 3PI effective action. The setup, which uses self-consistently calculated two-and three-point functions as input, is completely self-contained and does not have any free parameters except for the coupling. The results for J PC = 0 ±+ , 2 ±+ , 3 ±+ , 4 ±+ are in good agreement with recent lattice results where available. For the pseudoscalar glueball, we present first results from a two-loop complete calculation, rendering also the bound state calculation fully self-consistent.