This work introduces a numerical scheme for free vibration analysis of elastically supported piezoelectric nanobeam. Based on Hamilton principle, governing equations of the problem are derived. The problem is formulated for linear and nonlinear Winkler–Pasternak foundation type. Three differential quadrature techniques are employed to reduce the problem to an Eigen-value problem. The reduced system is solved iteratively. The natural frequencies of the beam are obtained. Numerical analysis is implemented to investigate computational characteristics affecting convergence, accuracy and efficiency of the proposed scheme. The obtained results agreed with the previous analytical and numerical ones. Furthermore, a parametric study is introduced to show influence of supporting conditions, two different electrical boundary conditions, material characteristics, foundation parameters, temperature change, external electric voltage, nonlocal parameter and beam length-to thickness ratio on the values of natural frequencies and mode shapes of the problem.