The origin of ultrahigh-energy cosmic-ray (UHECR) has been a mystery for many years. In order to identify candidate sources, it is crucial to figure out the energetics of cosmic-ray sources. By fitting the experiment data of cosmic-ray energy spectrum from Pierre Auger Observatory, we calculate the energy generation rate densities for different nuclei species via scanning over various parameters that are related to the source distribution models, injection properties as well as propagation of UHECR. Assuming UHECRs are protons and heavy nuclei produced by extragalactic sources, the energy generation rate density is estimated to be $\varepsilon^2Q$ $\sim 0.66^{+0.8}_{-0.3} \times$10$^{44}$ erg Mpc$^{-3}$ yr$^{-1}$ averaging different nuclei species at power-law spectral index $s = 2.0$ and it varies on different fitting range. For a power-law scaling of the CR production rate with redshift, maximum energy and pure proton injection, $dN/dE \propto E^{-s}\text{exp}(-E/E_\text{max})(1+z)^m$, best-fit parameters are $s=2.1$, $m=5.0$, $E_\text{max}=10^{19.9}$eV. Hard spectra with strong source redshift evolution and low maximal source energy are slightly favored. We did not find a strong proof of the linear relationship between $\varepsilon^2Q$ and $s$ as shown in Katz et al 2009. We discuss cases of mixed composition models by analyzing the combined fit of spectrum and mass composition $X_\text{max}$ data.