The damage and fracture behavior of Fiber Reinforced Polymers (FRPs) is quite complex and is different than the failure behavior of the traditionally employed metals. There are various types of failure mechanisms that can develop during the service life of composite structures. Each of these mechanisms can initiate and propagate independently. However, in practice, they act synergistically and appear simultaneously. The difficulties that engineers face to understand and predict how these different failure mechanisms result in a structural failure enforce them to use high design safety factors and also increases the number of certification tests needed. Considering that the experimental investigations of composites can be limited, very expensive, and time-consuming, in this contribution the newly developed multi Phase-Field (PF) fracture model [1] is employed to numerically study the failure in different Unidirectional Fiber Reinforced Polymers (UFRPs) laminates, namely, fracture in single-edge notched laminated specimens, matrix cracking in cross-ply laminates, and delamination migration in multi-layered UFRPs. The formulation of the PF model incorporates two independent PF variables and length scales to differentiate between fiber and inter-fiber (matrix-dominated) failure mechanisms. The physically motivated failure criterion of Puck is integrated into the model to control the activation and evolution of the PF parameters. The corresponding governing equations in terms of variational formulation is implemented into the Finite Element (FE) code ABAQUS utilizing the user-defined subroutines UMAT and UEL.