A micro-electromechanical constitutive model for ferroelectric polycrystals is introduced and applied to determine possible uniaxial remanent strain and remanent polarization states in the material. Simulations based on the model are carried out to determine the polarization saturation states under different levels of prevailing remanent strain. The simulation results clearly demonstrate how the polarization develops to the saturation and how the saturation polarization depends upon the remanent strain. The obtained information of saturation states can be applied in phenomenological constitutive models.