This paper is devoted to a detailed analysis of the appearance of frequency combs in the dynamics of a micro-electro-mechanical systems (MEMS) resonator featuring 1:2 internal resonance. To that purpose, both experiments and numerical predictions are reported and analysed to predict and follow the appearance of the phononic frequency comb arising as a quasi-periodic regime between two Neimark-Sacker bifurcations. Numerical predictions are based on a reduced-order model built thanks to an implicit condensation method, where both mechanical nonlinearities and electrostatic forces are taken into account. The reduced order model is able to predict a priori, i.e. without the need of experimental calibration of parameters, and in real time, i.e. by solving one or two degrees-of-freedom system of equations, the nonlinear behaviour of the MEMS resonator. Numerical predictions show a good agreement with experiments under different operating conditions, thus proving the great potentiality of the proposed simulation tool. In particular, the bifurcation points and frequency content of the frequency comb are carefully predicted by the model, and the main features of the periodic and quasi-periodic regimes are given with accuracy, underlining that the complex dynamics of such MEMS device is effectively driven by the characteristics of the 1:2 internal resonance.

The aim of this contribution is to present numerical comparisons of model-order reduction methods for geometrically nonlinear structures in the general framework of finite element (FE) procedures. Three different methods are compared: the implicit condensation and expansion (ICE), the quadratic manifold computed from modal derivatives (MD), and the direct normal form (DNF) procedure, the latter expressing the reduced dynamics in an invariant-based span of the phase space. The methods are first presented in order to underline their common points and differences, highlighting in particular that ICE and MD use reduction subspaces that are not invariant. A simple analytical example is then used in order to analyze how the different treatments of quadratic nonlinearities by the three methods can affect the predictions. Finally, three beam examples are used to emphasize the ability of the methods to handle curvature (on a curved beam), 1:1 internal resonance (on a clamped-clamped beam with two polarizations), and inertia nonlinearity (on a cantilever beam).