Observations and semianalytical galaxy formation and evolution models (SAMs) have suggested the existence of a stellar mass-stellar metallicity relation (MZR), which is shown to be universal for different types of galaxies over a large range of stellar masses ($M_*\sim 10^3$-$10^{11}M_\odot$) and dark matter (DM) halo masses ($M_{\rm halo}\sim 10^9$-$10^{15}h^{-1}M_\odot$). In this work, we construct a chemical evolution model to investigate the origin of the MZR, including both the effects of gas inflows and outflows in galaxies. We solve the MZR from the chemical evolution model, by assuming that the cold gas mass ($M_{\rm cold}$) and the stellar feedback efficiency ($\beta$) follow some power-law scaling relationships with $M_*$ during the growth of a galaxy, i.e., $M_{\rm cold}\propto M_*^{\alpha_{\rm gs}}$ and $\beta\propto M_*^{\alpha_{\beta{\rm s}}}$. We use the SAM to obtain these power-law scaling relations, which appear to be roughly universal over a large range of stellar masses for both satellites and central galaxies within a large range of halo masses. The range of the MZRs produced by our models is in a narrow space, which provides support to the universality of the MZRs. The formation of the MZR is a result caused jointly by that the cold gas fraction decreases with increasing $M_*$ and by that the stellar feedback efficiency decreases with increasing $M_*$ in the galaxy growth, and the exponent in the MZR is around $-\alpha_{\beta{\rm s}}$ or $1-\alpha_{\rm gs}$. The MZR represents an ''average'' evolutional track for the stellar metallicity of a galaxy. The comparison of our model with some previous models for the origin of MZRs is also discussed.

Observations reveal a universal stellar mass-stellar metallicity relation (MZR) existing in Local Group dwarfs of different types, $Z_*\propto M_*^{\alpha}$ with $\alpha=0.30\pm0.02$. In this work, we investigate the "universality" of the MZRs for both satellites and central galaxies in a large number of different host dark matter halos covering a large mass range of $10^9$-$10^{15}h^{-1}M_\odot$, by using a semianalytical galaxy formation and evolution model. We obtain the following results. (1) The exponents ($\alpha$) for the MZRs of the satellites in halos with the same mass as the Milky Way halo but different individual assembly histories are mostly $\sim$0.2-0.4, i.e., having a scatter of $\sim 0.2$; and the scatter of $\alpha$ increases with decreasing halo masses. (2) The MZR relations are changed little by the variation of halo masses and the classification between central galaxies and satellites, if many halos with the same mass are stacked together. (3) A double power law exists in the MZR relations for both central galaxies and stacked satellites, with $\alpha\sim$0.2-0.4 at $10^3M_\odot< M_*<10^{8} M_\odot$ and a relatively higher $\alpha\sim0.5$ at $10^8M_\odot<M_*<10^{11}M_\odot$. (4) The high-mass satellites ($M_*>10^8M_\odot$) existing mostly in high-mass halos can lead to an apparent increase of $\alpha$ (from $\sim0.2$ to $\sim0.4$) with increasing host halo masses shown in the single power law fitting results of stacked satellites. The universality of the MZR suggests the common physical processes in stellar formation and chemical evolution of galaxies can be unified over a large range of galaxy masses and halo masses.