Desire to accurately predict the deformation behaviour throughout industrial forming processes, such as thermoforming and stretch blow moulding, has led to the development of mathematical models of material behaviour, with the ultimate aim of embedding into forming simulations enabling process and product optimization. Through the use of modern material characterisation techniques, biaxial data obtained at conditions comparable to the thermoforming process was used to calibrate the Buckley material model to the observed non-linear viscoelastic stress/strain behaviour. The material model was modified to account for the inherent anisotropy observed between the principal directions through the inclusion of a Holazapfel–Gasser–Ogden hyperelastic element. Variations in the post-yield drop in stress values associated with deformation rate and specimen temperature below the glass transition were observable, and facilitated in the modified model through time-temperature superposition creating a linear relationship capable of accurately modelling this change in yield stress behaviour. The modelling of the region of observed flow stress noted when above the glass transition temperature was also facilitated through adoption of the same principal. Comparison of the material model prediction was in excellent agreement with experiments at strain rates and temperatures of 1–16 s−1 and 130–155 °C respectively, for equal-biaxial mode of deformation. Temperature dependency of the material model was well replicated with across the broad temperature range in principal directions, at the reference strain rate of 1 s-1. When concerning larger rates of deformation, minimum and maximum average error levels of 6.20% and 10.77% were noted. The formulation, and appropriate characterization, of the modified Buckley material model allows for a stable basis in which future implementation into representative forming simulations of poly-aryl-ether-ketones, poly(ether-ether-ketone) (PEEK) and many other post-yield anisotropic polymers.