In variable amplitude fatigue of high strength low alloy (HSLA) steel components, overloads can severely retard subsequent crack propagation for a number of cycles. In order to be able to predict fatigue crack propagation with a reduced degree of conservatism, retardation has to be taken into account. Of all numerical models that have been developed over time, crack tip plasticity models are selected based on the need for a detailed and fast cycle-by-cycle simulation of high cycle. After introducing the load interaction zone concept, common to all crack tip plasticity models, the Wheeler and Willenborg models are discussed, implemented and compared to experimental data. It is concluded that the Modified Wheeler model provides the most promising results, whereas the main limitation of Willenborg models is the need for extensive experimental data.