Modeling of rate-dependent hysteresis in piezoelectric actuators using a family of ellipses