Linear and nonlinear viscoelastic behaviors of a Zr55Al10Ni5Cu30 bulk metallic glass are investigated through experiments and described by the fictive stress model. Systematic deviations between the predicted stress-strain curves by fictive stress model and by the experimental results were found. In order to describe the flow stress curves of the Zr55Al10Ni5Cu30 BMG at different temperatures and strain rates in the supercooled liquid region more precisely, the fictive stress model was modified. The parameters of the model were optimized by the genetic algorithm, and a time relaxation factor Z' was introduced. The comparisons of the predicted compressive stress-strain curves and extrusion load-punch stroke curve by the modified fictive stress model with the experimental data show good agreements.