The prediction of failure and optimisation of the initial blank shape of a component formed under the complex thermal and mechanical conditions could not simply rely on its thinning or strain. In contrast, the prediction of forming limit of material as a comprehensive result of temperature, strain rate and loading path is an essential indicator in a non-isothermal hot stamping process. However, it is challenging to implement an advanced forming limit model into conventional finite element (FE) software to simulate the failure of the formed component and thus optimise its initial blank shape. In this work, a unified viscoplastic-Hosford-MK model was developed to post-process the simulated evolutionary thermomechanical characteristics and then numerically predict the forming limit criterion of AA7075 in a hot stamping process, of which results were demonstrated on a dedicated visualiser that independent of FE software. Subsequently, the initial blank shape was optimised by analysing the evolutionary temperature, strain rate and loading path of failure regions. Furthermore, the developed forming limit model and optimised initial blank shape were experimentally verified with an error of less than 10%.