Abstract: In order to overcome the disadvantage of traditional cohesive zone damage model (CZM) that it cannot present the effect of intralaminar cracks on delamination interface, a modified CZM has been developed for simulation of damage process of composite laminates subject to low-velocity impact. By means of the modification of the damage initiation criteria of cohesion constitutive model used in interface element, the influences of damage status and stress distribution of matrix and fiber in neighboring interface layers on interlaminar strength and delamination propagation have been taken into account. The proposed model, combined with intralaminar failure criterion, was implemented to ABAQUS as a user defined subroutine VUMAT to simulate the progressive damage process of composite laminate under low-velocity impact. In order to valid the accuracy and rationality of this method, the damage status of composite laminates with different ply orientations and material properties are compared with test results.