This work is a collaboration between EDF R&D and the Laboratory of Mechanics and Acoustics. The objective is to develop theoretical and numerical tools to compute nonlinear normal modes (NNMs) of structures with localized nonlinearities. We use an approach combining the harmonic balance and the asymptotic numerical methods, known for its robustness principally for smooth systems. Regularization techniques are used to apply this approach for the study of nonsmooth problems. Moreover, several aspects of the method are improved to allow the computation of NNMs for systems with a high number of degrees of freedom (DOF). Finally, the method is implemented in Code_Aster, an open-source finite element solver developed by EDF R&D. The nonlinear normal modes of a two degrees-of-freedom system are studied and some original characteristics are observed. These observations are then used to develop a methodology for the study of systems with a high number of DOFs. The developed method is finally used to compute the NNMs for a model U-tube of a nuclear plant steam generator. The analysis of the NNMs reveals the presence of an interaction between an out-of-plane (low frequency) and an in-plane (high frequency) modes, a result also confirmed by the experiment. This modal interaction is not possible using linear modal analysis and confirms the interest of NNMs as a diagnostic tool in structural dynamics.