Extending the EGP constitutive model for polymer glasses to multiple relaxation times