Modeling the muscle response to functional electrical stimulation (FES) is an important step during model-based FES control system design. The Hammerstein structure is widely used in simulating this nonlinear biomechanical response. However, a fixed relationship cannot cope well with the time-varying property of muscles and muscle fatigue. In this paper, we proposed an adaptive Hammerstein model to predict ankle joint torque induced by electrical stimulation, which uses variable forgetting factor recursive least squares (VFFRLS) method to update the model parameters. To validate the proposed model, ten healthy individuals were recruited for short-duration FES experiments, ten for long-duration FES experiments, and three stroke patients for both. The isometric ankle dorsiflexion torque induced by FES was measured, and then the test performance of the fixedparameter Hammerstein model, the adaptive Hammerstein model based on fixed forgetting factor recursive least squares (FFFRLS) and the adaptive Hammerstein model based on VFFRLS were compared. The goodness of fit, root mean square error, peak error and success rate were applied to evaluate the accuracy and stability of the model. The results indicate a significant improvement in both the accuracy and stability of the proposed adaptive model compared to the fixed-parameter model and the adaptive model based on FFFRLS. The proposed adaptive model enhances the ability of the model to cope with muscle changes.