The application possibility for the model of ferroelectroelastic material based on the testing of effective properties of polycrystal using the methods of finite element homogenization to describe the cyclic loading of polycrystalline piezoceramics has been investigated. The results of finite element solution for nonlinear coupled electro-mechanical problems for a representative volume element of the polycrystal were presented. They were compared with experimental results under various loading conditions.