A model for second-harmonic generation from materials with inversion symmetry is developed. The reflecting medium is treated as an ordered array of localized, polarizable entities and the self-consistent relations between induced multipole moments and local fields are derived. The resulting equations are evaluated for a simple model of an entity's response, yielding microscopic estimates of both surface and bulk susceptibility components. The difficulties of applying the model to describe real materials are discussed.