In this study, a Monte Carlo simulation (MCs)-based isogeometric stochastic Finite Element Method (FEM) is proposed for uncertainty quantification in the vibration analysis of piezoelectric materials. In this method, deterministic solutions (natural frequencies) of the coupled eigenvalue problem are obtained via isogeometric analysis (IGA). Moreover, MCs is employed to solve various uncertainty parameters, including separate elastic and piezoelectric constants and their combined cases.