Using a nonlinear thermodynamic theory, we describe equilibrium polarization states and the macroscopic dielectric response of nanocrystalline ferroelectric ceramics with single-domain grains. The elastic clamping of individual crystallites by the surrounding material is explicitly taken into account via the introduction of a specific thermodynamic potential. Aggregate material properties are calculated with the aid of an iterative procedure based on the method of effective medium. The numerical calculations, performed for unpolarized BaTiO3 and Pb(Zr1−xTix)O3 ceramics, demonstrate that the equilibrium phase states of nanocrystalline ceramics may differ drastically from those of single crystals and coarse-grained materials. Remarkably, the theory predicts the coexistence of rhombohedral and tetragonal crystallites in nanocrystalline Pb(Zr1−xTix)O3 ceramics in a wide range of compositions and temperatures. For BaTiO3 ceramics, a mixture of rhombohedral and orthorhombic crystallites is found to be the energetically most favorable state at room temperature. The calculations also show that the dielectric properties of nanocrystalline ferroelectric ceramics may be very different from those of conventional materials due to the elastic clamping of single-domain crystallites.