Linear modal analysis is a mature tool enjoying various applications ranging from bridges to satellites. Nevertheless, modal analysis fails in the presence of nonlinear dynamical phenomena and the development of a practical nonlinear analog of modal analysis is a current research topic. Recently, numerical techniques (e.g., harmonic balance, continuation of periodic solutions) were developed for the computation of nonlinear normal modes (NNMs). Because these methods are limited to conservative systems, the present study targets the computation of NNMs for nonconservative systems. Their definition as invariant manifolds in phase space is considered. Specifically, a new finite element technique is proposed to solve the set of partial differential equations governing the manifold geometry.