We present a generalized variational procedure oriented to the algebraic solution of many-body Hamiltonians expressed in bosonic and fermionic variables. The method specializes in the nonperturbative regime of the solutions. As an example, we focus on the application of the method to nonperturbative QCD.