With the effects of electronic structure and atomic size being introduced, a revised model to calculate the viscosity of the bulk metallic glass alloys was proposed and the viscosity of ternary Zr-Al-Cu, Zr-Ni-Al and quaternary Zr-Al-Ni-Cu systems are calculated in this paper, and the computed results agree well with the empirical one. The sequence of viscosity of different systems is: VZr-Al-Cu <VZr-Al-Ni<.VZr-Al-Ni-Cu. To Zr-Al-Cu and Zr-Ni-Al, the highest viscosity locates in the composition range of XZr=0.37-0.86, XCu=0-0.40 and XZr = 0.45-0.79, XAl = 0.12-0.50, respectively. And to the Zr-Ni-Al-Cu system with 66.67% Zr, the highest viscosity is obtained in the region of XAl= 0.63-0.80, XNi = 0.14-0.24.