Abstract The prediction of the modal properties of structural systems with contact clearances and prestress presents a computational challenge, as the nonlinearity induced by piecewise-linear stiffness eliminates the use of efficient linear modal analysis techniques. The most common approach to obtaining the nonlinear normal modes (NNMs) of these structural systems is a numerical framework that integrates numerical integration, the shooting method, and the pseudo-arc-length continuation scheme. This numerical continuation framework (NCF) computes NNMs through iterative numerical calculations; thus, the computational cost of the nonlinear modal analysis of complex nonlinear systems, particularly piecewise-linear systems, becomes prohibitively expensive as the model size increases. In this study, a hybrid continuation framework (HCF) combining analytic and numerical methods is proposed to enable efficient computations of NNMs for systems with contact boundaries. This new hybrid framework utilizes a semi-analytic method to conduct the iterative shooting procedure; thus, the computational burden of the numerical continuation can be significantly reduced. The proposed method is demonstrated on spring-mass oscillators with contact elements, and the NNMs obtained using the proposed method are validated by those computed using the traditional numerical continuation framework. The modal properties of the systems can be computed using the proposed framework with significant speed-up. Furthermore, the modal properties, including internal resonance and sharp turning in NNM curves, of the piecewise-linear systems are identified and discussed.