Evaluation of dielectric and piezoelectric behavior of unpoled and poled barium titanate polycrystals with oxygen vacancies using phase field method