This paper introduces a numerical model to estimate fatigue life under step‐stress conditions, using the Weibull and lognormal distributions. The maximum likelihood method was used to estimate the free parameters of the distributions. The model was fitted to an experimental data on fatigue life in the specimens of steel SAE 8620, by using evolutionary computation to optimize the likelihood function. Results are reported on the values of the parameters and their confidence interval. Also, a validation of the model is discussed using analysis of residuals.