This paper concerns a square lattice, Ising-type model with interactions between the four spins at the corners of each face. These may include nearest and next-nearest-neighbor interactions, and interactions with a magnetic field. Provided the Hamiltonian is symmetric with respect to both row reversal and column reversal, a rapidly convergent sequence of variational approximations is obtained, giving the free energy and other thermodynamic properties. For the usual Ising model, the lowest such approximations are those of Bethe and of Kramers and Wannier. The method provides a new definition of corner transfer matrices.