For nonlinear problems, it is shown that the accuracy and stability of superconvergent patch recovery (SPR) method can be remarkably improved through the modi!cations introduced in this paper. These modi!cations are (1) use of integration points as sampling points, (2) weighted average procedure and (3) introduction of additional nodes. They are validated through numerical investigations of several large deformation problems. Based on the modi!ed SPR procedure, a new scheme for submodeling !nite element analysis is developed. Comparison of the numerical results obtained from this new scheme and traditional method shows encouraging improvements.