A chemical model is proposed for a dusty plasma consisting of electrons, ions, neutrals, and positively charged dust particles all being at thermodynamic equilibrium. An expression is derived for the Helmholtz free energy, which comprises the ideal part, taking into account the charge of dust particles, and the excess part, evaluated in the framework of the self-consistent chemical model [Phys. Rev. E 83, 016405 (2011)PLEEE81539-375510.1103/PhysRevE.83.016405]. The ionization potential depression for a dust-free partially ionized hydrogen is analytically evaluated for weakly and strongly ionized states to consistently account for the presence of charged and neutral components. An ad hoc interpolation of the ionization potential depression, valid across the whole ionization region, is put forward and subsequent solution of the generalized Saha equation is found to be in a perfect agreement with exact calculations. Minimization of the Helmholtz free energy of dusty plasmas provides the number densities of free electrons, free ions, neutrals, and the dust electric charge as well. Based on consideration of weakly and strongly ionized states, a straightforward comparison is made of the ionization equilibrium in a partially ionized plasma with and without dust particles to demonstrate that at thermal equilibrium positively charged dusts are responsible for an increase in the number density of free electrons and a decrease in the number density of free ions. It is analytically proved that nonideality effects result in a growth of the number densities of free electrons and ions by obtaining the so-called electron and proton ionization potential depressions. Electric charge of dust particles is systematically studied as a full plasma component rather than considering a detailed balance of the electron and ion fluxes on the surface of a solitary dust grain.