To date, the complex behaviour of small punch creep test (SPCT) specimens has not been completely understood, making the test hard to numerically model and the data difficult to interpret. This paper presents a novel numerical model able to generate results that match the experimental findings. For the first time, pre-strained uniaxial creep test data of a P91 steel at 600 • C have been implemented in a conveniently modified Liu and Murakami creep damage model in order to simulate the effects of the initial localised plasticity on the subsequent creep response of a small punch creep test specimen. Finite element (FE) results, in terms of creep displacement rate and time to failure, obtained by the modified Liu and Murakami model are in good agreement with experimental small punch creep test data. The rupture times obtained by the FE calculations which make use of the non-modified creep damage model are one order of magnitude shorter than those obtained by using the modified constitutive model. Although further investigation is needed, this novel approach has confirmed that the effects of initial localised plasticity, taking place in the early stages of small punch creep test, cannot be neglected. The new results, obtained by using the modified constitutive model, show a significant improvement with respect to those obtained by a 'state of the art' creep damage constitutive model (the Liu and Murakami constitutive model) both in terms of minimum load-line displacement rate and time to rupture. The new modelling method will