There are several models to estimate the dynamic response of a cantilever-type piezoelectric harvester in the literature. One of them being through Iso-Geometric Analysis (IGA) framework for a PEH [1]}, which stands out due to its higher accuracy reached in reduced computational time. This is the reason why this model sis selected for the present work. Figure REF  presents a schematic of the piezoelectric energy harvester adopted in this study.
<FIGURE>