Fast method for the calculation of surface bending on circular multilayered piezoelectric structures