Theoretical Modeling of Frequency Combs Generated in Piezoelectric Nanomechanical Resonators