For an ideal piezoelectric bimorph, assuming a third order deflection equation for the cantilever beam, the equilibrium equation can be written as shown in attached. The first matrix is the stiffness matrix which is multiplied by the generalized coordinates vector.