A hybrid model is proposed in this paper to describe the static nonlinear and dynamic characteristics of rate-dependent hysteresis in piezoelectric actuators. In this model, a neural network based submodel is implemented to approximate the static nonlinear characteristic of the hysteresis while a submodel with the first-order differential operators is constructed to describe the dynamic behavior of the hysteresis. In this paper, a special hysteretic operator is proposed to extract the hysteretic feature of the hysteresis. Then, an expanded input space with such special hysteretic operator is constructed. Based on the constructed expanded input space, the neural network can be implemented to approximate the hysteresis phenomenon. The submodel to describe the dynamics is a sum of the weighted first-order differential operators. Finally, the experimental results of applying the proposed method to the modeling of hysteresis in a piezoelectric actuator are also presented.