Current methodologies for variation simulation of compliant sheet metal assemblies and parts are simplified by assuming linear relationships. From the observed physical experiments, it is evident that plastic strains are a source of error that is not captured in the conventional variational simulation methods. This paper presents an adaptation toward an elastoplastic material model with isotropic hardening in the method of influence coefficients (MIC) methodology for variation simulations. The results are presented in two case studies using a benchmark case involving a two-dimensional (2D) quarter symmetric plate with a centered hole, subjected to both uniaxial and biaxial displacement. The adaptation shows a great reduction in central processing unit time with limited effect on the accuracy of the results compared to direct Monte Carlo simulations.