Phase field modeling of deformation twinning in β-metastable titanium alloys