This study aims to evaluate the applicability of failure criteria and damage evolution methods in the finite element analysis of composite laminates under lowvelocity impact. Implemented by the user-defined VUMAT subroutine in ABAQUS, various progressive damage models are used to predict damage initiation and accumulation in a T700GC/M21 composite laminate. Cohesive elements are inserted between adjacent plies to capture interface delamination. The applicability of damage models is investigated by comparing the global mechanical response and distribution of various damage modes. A new variable, equivalent damage volume, is introduced to quantitatively describe the predicted damage when using different models. The numerical results establish that Hashin and Puck failure criteria generate matrix compression damage in more layers of the composite. Maximum stress and Tsai-Wu criteria are not preferred due to their improper predictions in terms of damage area and permanent deformation of the laminate. As for damage evolution laws, the equivalent strain method provides faster stiffness degradation of the laminate and a smaller area 2 of matrix damage compared with the predictions of the equivalent displacement method.