The model formulation is based on the Kirchhoff-Love plate theory and Hamilton's generalized principle for electro-mechanical bodies. The model is solved numerically using the IGA method, where B-Splines \(N_I\)  describe the device’s domain, and the deflection field \(\textbf {w}\)  is approximated using the same basis functions \(N_I\) . The framework was presented in [1]}, where the relative deflection solution \(\textbf {w}\)  can be approximated by a truncated expansion of the first \(K\)  mode shape vectors \(\textbf {w}_o = \Phi _o \eta \) . Here, \(\Phi _o\in \mathbf {R}^{N \times K}\)  is the matrix which contains the \(K\)  first mode shape vectors \(\phi _i\)  and \(\eta \in \mathbf {R}^{K\times 1}\)  denotes the modal coordinates. Therefore, the procedure leads to a coupled system of differential equations of the form,
\(\ddot{\eta }+\textbf {c}_o\dot{\eta }+ \textbf {k}_o\eta  - \theta _ov = \textbf {f}_o a_b(t)\) 
\(C_p\dot{v}(t) + \frac{v(t)}{R_l} + \Theta ^T\Phi _o \dot{\eta }=0\) 
