Coupled phase field simulations of ferroelectric and ferromagnetic layers in multiferroic heterostructures