Quantifying the dynamic near-fault stress evolution associated with fault shear deformation is crucial for revealing the earthquake nucleation mechanisms of rough faults. However, it is still challenging via traditional methods because of the difficulties in resolving continuously evolved stresses and high stress heterogeneity. In this study, we present a novel method to quantify the dynamic full-field stress evolution using three-dimensional printed models and photoelastic techniques. The dynamic full-field stresses in a circular disc determined by this method and an analytical model are compared to verify the accuracy of the proposed method. A new method to count the fringe orders near the rough fault is proposed to achieve the stress quantification. The full-field principal stress difference in a rough fault model under shear forces is derived. Our findings indicate that the proposed method can significantly improve the efficiency of stress calculation and effectively quantify the dynamic full-field stress near rough faults.