This article discusses the rheological tests and analyses based on the Schapery non-linear viscoelasticity model that were performed to study asphalt mastic behaviour under high shear stresses. Seven mineral filler types were applied in this study, including a mixed filler with hydrated lime and fillers derived from dust extraction systems. Determination of basic properties of the fillers was followed by creep and recovery tests (DSR) at different levels of shear stress conducted in accordance with a modified MSCR procedure. The first stage in the analysis was the identification of linear viscoelastic region and the non-linear viscoelasticity model parameters such as the length of the loading period, the temperature and the stress level using TTSSP (Time-Temperature- Stress Superposition Principle). Subsequent numerical simulations of strain variation with respect to stress confirmed a high degree of agreement between the non-linear viscoelasticity model and mastic sample behaviour. A strong correlation was found between the non-linear viscoelasticity parameters and mastic properties. The proposed methodology is able to quickly identify and eliminate the fillers that may contribute to HMA deformations.