A rate-dependent model for sand to predict constitutive response and instability onset