A lattice-Monte Carlo approach was developed to simulate ferroelectric domain behavior. The model utilizes a Hamiltonian for the total energy that includes electrostatic terms (involving dipole-dipole interactions, local polarization gradients, and applied electric field), and elastic strain energy. The contributions of these energy components to the domain structure and to the overall applied field response of the system were examined. In general, the model exhibited domain structure characteristics consistent with those observed in a tetragonally distorted ferroelectric. Good qualitative agreement between the appearance of simulated electrical hysteresis loops and those characteristic of real ferroelectric materials was found. ©2000 American Institute of Physics.