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.