A method is developed to pred ict numerically the damage of composite laminates with more plies under low velo city impact loading. The Puck criterion for 3D stress states is adopted to mode l the intralaminar damage including matrix cracking and fiber breakage, and to o btain the orientation of the fracture plane due to matrix failure. According to interlaminar delamination mechanism, a new delamination criterion is proposed. The influence of transverse and through-thickness normal stress, interlaminar sh ear stress and damage conditions of adjacent plies on delamination is considere d. In order to quickly and efficiently predict the impact-induced damage of composite laminates with mor e plies, an approach, which can predict the specific dam age of several plies in a single solid element, is proposed by interpolation on the strains of element integration points. Moreover, the proposed model can pre dict specific failure modes. A good agreement between the predicted delamination shapes and sizes and the experimental results shows correctness of the develope d numerical method for predicting low velocity impact damage on composite laminates.