  The elastostatic Green function tensor approach, which was recently used to
treat electrostriction in numerical simulation of domain structure formation in
cubic ferroelectrics, is reviewed and extended to the crystals of hexagonal and
orthorhombic symmetry. The tensorial kernels appearing in the expressions for
effective nonlocal interaction of electrostrictive origin are derived
explicitly and their physical meaning is illustrated on simple examples. It is
argued that the bilinear coupling between the polarization gradients and
elastic strain should be systematically included in the Ginzburg-Landau free
energy expansion of electrostrictive materials.
