The numerical simulation is an effective method for researching the relationship between macroscopic properties of ZnO varistors and their microstructural characteristics. However, the computing efficiency is a big problem for further researches. In this paper, a quick piecewise linear method for solving the complicated nonlinear resistor networks of ZnO varistors is proposed, which can significantly improve the computing efficiency by a hundred times. The simulating models for the microstructural characteristics of ZnO varistors are improved based on actual conduction mechanism of grain boundaries. The simulations of more bulk characteristics of ZnO varistors become possible, such as alternating current (AC) and surge properties.