Actuators based on piezoelectric materials have ideal characteristics for applications such as acoustic transmission and micromanipulation. However, the inherent nonlinearities of those actuators, such as hysteresis and creep, greatly increase the challenge to control such devices. Furthermore, the increasing need for more precise and faster actuators, allied with frequent changes in the environmental and operational conditions, further worsens the problem. Analytical models are application-specific, meaning that they are not easily and efficiently scalable to all systems. Also, with increased complexity, the understating of underlying phenomena is not fully documented, making it difficult to develop such models. This work investigates those challenges from the perspective of the system identification methodology and data-driven models for piezoelectric actuators. The black-box approach is tested with experimental data acquired in a laboratory setting for micromanipulator and acoustic transmission case studies. In some datasets, general-purpose signals were employed as the excitation input of the system to accelerate the data acquisition of the whole system dynamic and estimation process. Additionally, some models were validated on a separate dataset. In both cases, preprocessing was employed to optimize the amount of data. The tested models include the AutoRegressive Moving Average with eXogenous inputs (ARMAX), Nonlinear AutoRegressive with eXogenous inputs (NARX) with an artificial neural network structure, and Nonlinear AutoRegressive Moving Average with eXogenous inputs (NARMAX). The results show a good ability to predict the nonlinearities of the micromanipulator and, therefore, the hysteresis at different input frequencies. The acoustic transmission system was successfully modeled. Although the results show that there is still room for improvements, it provides insights into possible optimizations for the setup as the models here devised are useful for short prediction windows.