The lead-free Ba(Zr0.2Ti0.8)O3-x(Ba0.7Ca0.3)TiO3 (BZT-xBCT) solution has been reported to show excellent piezoelectric performance comparable to that of soft PZT at room temperature (e.g., d33 = 580~620 pC/N for BZT-50BCT at room temperature). For the BZT-xBCT system, the piezoelectric properties can be further optimized by ion doping. In this work, we consider variants of this system by including additional dopants and develop parameterized phase field simulations to study domain evolution. We investigate switching with increasing field to study factors leading to large piezoelectric enhancement. Our simulation results indicate that the domain rotation behavior can be greatly influence by dopants. Small changes of the elastic constants due to ion doping can lead to about 50% enhancement of the piezoelectric property.