Loss phenomena in piezoceramics presents a major challenge in modern resonant and off-resonant driven actuator systems for high-precision positioning tasks. Hysteretic material behavior, which manifests as phase delay in deflection and current, requires complex control strategies in piezo-based systems. To achieve an increase in performance, various modeling approaches can be used to describe the loss characteristics of the material. A simple and popular approach for engineering applications is the use of complex material constants, which are representative of the phase delay. This frequency-domain approach, however, has the disadvantage that it is only valid for small dissipation rates and is also unsuitable for extended multiscale analyses. For a mechanism-based investigation of losses and the use of arbitrarily material models coupled to the microscale, a consideration in the time domain is necessary. The frequency-, temperature- and amplitude-dependent loss behavior of piezoceramics can thus be more accurately represented using sophisticated rate-dependent models - however, this requires a high computational effort in dynamic applications. One way to avoid these high computational costs is to use the harmonic balance method, which allows for the prediction of the complex, frequency-dependent loss behavior of piezoelectric components with moderate computational effort.