Finite-Element Modeling of the Hysteresis Behavior of Tetragonal and Rhombohedral Polydomain Ferroelectroelastic Structures.

The influence of the domain structure's initial topology and its evolution on the hysteresis curves of tetragonal and rhombohedral polydomain structures of ferroelectroelastic materials is studied. Based on the analysis of electrical and mechanical compatibility conditions, all possible variants of representative volume elements of tetragonal and rhombohedral second-rank-domain laminate structures were obtained and used in simulations. Considerable local inhomogeneity of stress and electric fields within the representative volume, as well as domain interaction, necessitates the use of numerical methods. Hysteresis curves for laminated domain patterns of the second rank were obtained using finite-element homogenization. The vector-potential finite-element formulation as the most effective method was used for solving nonlinear coupled boundary value problems of ferroelectroelasticity. A significant anisotropy of the hysteresis properties of domain structures was established both within individual phases and when comparing the tetragonal and rhombohedral phases. The proposed approach describes the effects of domain hardening and unloading nonlinearity.