Finite element reduced order models for nonlinear vibrations of piezoelectric layered beams with applications to NEMS

This article presents a finite element reduced order model for the nonlinear vibrations of piezoelectric layered beams with application to NEMS. In this model, the geometrical nonlinearities are taken into account through a von Kármán nonlinear strain–displacement relationship. The originality of the finite element electromechanical formulation is that the system electrical state is fully described by only a couple of variables per piezoelectric patches, namely the electric charge contained in the electrodes and the voltage between the electrodes. Due to the geometrical nonlinearity, the piezoelectric actuation introduces an original parametric excitation term in the equilibrium equation. The reduced-order formulation of the discretized problem is obtained by expanding the mechanical displacement unknown vector onto the short-circuit eigenmode basis. A particular attention is paid to the computation of the unknown nonlinear stiffness coefficients of the reduced-order model. Due to the particular form of the von Kármán nonlinearities, these coefficients are computed exactly, once for a given geometry, by prescribing relevant nodal displacements in nonlinear static solutions settings. Finally, the low-order model is computed with an original purely harmonic-based continuation method. Our numerical tool is then validated by computing the nonlinear vibrations of a mechanically excited homogeneous beam supported at both ends referenced in the literature. The more difficult case of the nonlinear oscillations of a layered nanobridge piezoelectrically actuated is also studied. Interesting vibratory phenomena such as parametric amplification or patch length dependence of the frequency output response are highlighted in order to help in the design of these nanodevices.

► Finite element reduced order model for nonlinear vibrations of nanoresonators.
► A priori computation of the nonlinear stiffness coefficients of the low-order model.
► The model is computed with an original purely harmonic-based continuation method.
► Applications to layered beams with external and parametric piezoelectric actuation.
► Parametric study of the resonator geometry influence on the nonlinear resonances.