Stability, Evolution and Switching of Ferroelectric Domain Structures in Lead-free BaZr$_{0.2}$Ti$_{0.8}$O$_3$-Ba$_{0.7}$Ca$_{0.3}$TiO$_3$ System: Thermodynamic Analysis and Phase-field Simulations