Abstract: A modified Reissner-type incremental variational functional is proposed by adopting a updated Lagrangian description. The variational functional involves displacements and second Piola-Kirchhoff stress increments as independent variables. To enhance computational efficiency, the mechanical subelement model is adopted to describe nonlinear material behaviors. Three example problems are solved to test the proposed method, and the reliability and effectiveness are illustrated by the results.