In this work, modifications to existing energy-based nonlinear normal mode (NNM) methodologies are developed in order to investigate internal resonances. A formulation for computing resonant NNMs is developed for discrete, or discretized for continuous systems, sets of weakly nonlinear equations with uncoupled linear terms (i.e systems in modal, or canonical, form). By considering a canonical framework, internal resonance conditions are immediately recognized by identifying commensurable linearized natural frequencies. Additionally, the canonical formulation allows for a single (linearized modal) coordinate to parameterize all other (modal) coordinates during a resonant modal response. Energy-based NNM methodologies are then applied to the canonical equations and asymptotic solutions are sought. In order to account for the resonant modal interactions, it will be shown that high-order terms in the 0(1) solutions must be considered. Two applications (‘3:1’ resonances in a two-degree-of-freedom system and ‘3:1’ resonance in a hinged-clamped beam) are then considered by which to demonstrate the application of the resonant NNM methodology. Resonant normal mode solutions are obtained and the stability characteristics of these computed modes are considered. It is shown that for some responses, nonlinear modal relations do not exist in the context of physical coordinates and thus the transformation to canonical coordinates is necessary in order to define appropriate NNM relations.