---

This paper presents a finite element solution for the static analysis of a multi-layers beam with and without piezoelectric layers. The beam is under large deformation. The virtual work principle and the Lagrangian update method (LUM) have been employed to study the static behavior of piezoelectric beams. Four-nodes element with two displacement degrees of freedom and one electrical degree of freedom has been used in this analysis. Finally, in order to prove the validity of the presented formulation and the solving process, the results are compared with the other available data.