We report a new model, based on harmonic balancing, to predict the resonance peak shape and extract pre-strain and gauge factors of carbon nanotube-resonators. In addition, the model is able to separate transconductive and piezoresistive current contributions. The model predictions agree well with measurements performed on tunable carbon nanotube resonators throughout a wide range of frequencies [70, 85] MHz and static gate biases [-10, +10] V. In particular, the model correctly describes the changes in the peak shape when the sign of the static gate voltage is flipped, which is not the case with previous models in the literature that use a Lorentzian fit with an adhoc phase.