A phantom node method with mixed mode cohesive law is proposed for the simulation of splitting in laminates. With this method, a discontinuity in the displacement field can be modeled at arbitrary locations. The micromechanical phenomenon that splitting cracks grow parallel to the fiber, is incorporated on the mesolevel, i.e., in the homogenized ply, by setting the direction of the crack propagation equal to the fiber direction. A new mixed mode cohesive law is introduced for increased robustness of the incremental-iterative solution procedure. The model is validated with mixed mode bending tests, and its utility is illustrated with examples for a single ply and for a laminate.