The polyvinyl-alcohol-engineered cementitious composite (PVA-ECC) is a superior cementitious material when used for tension and flexural loading. The utilization of PVA-ECC in the tension zone can prevent the development of wide cracks and increase the flexural resistance of reinforced PVA-ECC members. In this paper, a nonlinear finite element model is established to simulate the behavior of PVA-ECC beams in bending. In the model, the constitutive models for PVA-ECC in compression and tension are employed by simplifying them as piece-wise linear models, and the bond between the reinforcing bar and PVA-ECC is also considered. The load–deflection curve and failure mode of beams can be obtained from the finite element model. Comparisons between numerical and experimental results show that the developed numerical model can estimate the ultimate load and failure mode of beams with reasonably good accuracy. After evaluating the accuracy of the finite element model, parameter analysis is conducted to investigate the effects of the reinforcement ratio, steel strength grade, and mechanical properties of PVA-ECC on the flexural behavior of reinforced PVA-ECC beams. The numerical results conclude that the effects of reinforcement ratio on the peak load, stiffness, and deflection are obvious while the influence of steel grade is mainly on the peak load. The tensile localization strain of PVA-ECC mainly affects the ductility of the beam. Furthermore, a design method is proposed based on the plane-section assumption to calculate the ultimate load of reinforced PVA-ECC beams, in which the contribution of PVA-ECC to the moment resistance of beam sections is considered. Comparisons between existing design methods and the proposed method indicate that the ultimate load of beams can be predicted more accurately by considering the tensile strength of PVA-ECC in the tension zone.