The perovskite ferroelectric materials, such as Pb(Zr_xTi_<1-x>)O_3 (PZT), consist of multiple crystal systems at morphotropic phase boundary (MPB) and they exert a huge piezoelectric response depending on domain switching and structural phase transition. In the present work, a multiscale nonlinear finite element method based on crystallographic homogenization theory was developed to evaluate nonlinear hysteresis behaviors at MPB. We employed the energy-based criteria and the incremental constitutive law taking into consideration with material property changes caused by domain switching and structural phase transition. The developed computation was applied to a PZT polycrystal and then the relation between macrostructural hysteresis behavior and microstructural crystal orientation change was discussed.