Simulations of domain evolution in morphotropic ferroelectric ceramics under electromechanical loading using an optimization-based model