Nonlinearities in piezoelectric systems can arise from internal factors such as nonlinear constitutive laws or external factors like realizations of boundary conditions. It can be difficult or even impossible to derive detailed models from the first principles of all the sources of nonlinearity in a system. This paper introduces adaptive estimator techniques to approximate the nonlinearities that can arise in certain classes of piezoelectric systems. Here an underlying structural assumption is that the nonlinearities can be modeled as continuous functions in a reproducing kernel Hilbert space (RKHS). This approach can be viewed as a data-driven method to approximate the unknown nonlinear system. This paper introduces the theory behind the adaptive estimator, discusses precise conditions that guarantee convergence of the function estimates, and studies the effectiveness of this approach numerically for a class of nonlinear piezoelectric composite beams.