In this study, vibration characteristics of a piezoelectric nanobeam embedded in a viscoelastic medium are investigated based on nonlocal Euler–Bernoulli beam theory. In doing this, the governing equations of motion and boundary conditions for vibration analysis are first derived using Hamilton’s principle, where nonlocal effect, piezoelectric effect, flexoelectric effect, and viscoelastic medium are considered simultaneously. Subsequently, the transfer function method is employed to obtain the natural frequencies and corresponding mode shapes in closed form for the embedded piezoelectric nanobeam with arbitrary boundary conditions. The proposed mechanics model is validated by comparing the obtained results with those available in the literature, where good agreement is achieved. The effects of nonlocal parameter, boundary conditions, slenderness ratio, flexoelectric coefficient, and viscoelastic medium on vibration responses are also examined carefully for the embedded nanobeam. The results demonstrate the efficiency and robustness of the developed model for vibration analysis of a complicated multi-physics system comprising piezoelectric nanobeam with flexoelectric effect, viscoelastic medium, and electrical loadings.