Polyethylene terephthalate (PET) has been widely used in the stretch blow moulding (SBM) process for packaging applications. Finite element analysis has become extensively useful for assessing container designs and enabling the designers to perform analyses earlier in the design cycle to determine the best material and the best structure. However, there are several challenging issues due to various processing parameters and complex material behaviour, which is both temperature and strain-rate dependent. In this paper, we generalize the G'Sell-Jonas law in the three-dimensional (3D) case to model and simulate the elasto-visco-plastic (EVP) behaviour of PET, taking into account strain-hardening and strain-softening. In addition, it is observed that the internal pressure (inside the preform) is significantly different from the nominal pressure (imposed in the blowing device upstream) since the internal pressure and the enclosed volume of the preform are fully coupled. In order to accurately simulate this phenomenon, a thermodynamic model was used to characterize the pressure-volume relationship (PVR). The predicted pressure evolution is therefore more realistic when imposing only the machine power of the blowing device (air compressor or vacuum pump). Mechanical and temperature equilibrium equations are fully nonlinear and solved separately with implicit schemes on the current deformed configuration, which is updated at each time step. Biaxial characterization tests were used to determine the model parameters in order to simulate the SBM process using the PVR. Three industrial case studies, comparing simulated thickness predictions to experimental measurements, will be presented in order to illustrate the applicability of the proposed model.