Discrete damage modeling techniques are superior in predicting failure of composite laminates involving massive delamination. In order to reproduce the influence of transverse crack opening on local stress distribution, identically orientated ply groups are often meshed with multiple layers of elements when applying these techniques to model transverse crack evolution processes of laminates. To save computation costs, researchers also tried to use only a single layer of elements in cracking plies whose accuracy, however, still needs to be improved. In the current study, an improved modeling strategy is proposed where adjacent identically-oriented plies are meshed with a single layer of elements. Herein, a phantom node solid element is adopted for matrix cracking modeling and a phantom node cohesive element is proposed to match the intraply element. Soft interfaces and a Weibull model of transverse strength are introduced in order to reproduce experimental phenomena and a parameter correction method for interface properties is proposed based on a shear lag model to improve predictions. The proposed strategy is verified by modeling crack evolution processes of cross-ply laminates of two different material systems. Besides, influences of model parameters and the applicability of identified parameters to cases of other layups are investigated.