Peridynamics is a robust theory to capture failure initiation and failure propagation in both isotropic and orthotropic materials. This paper present a method to model Mod I delamination failure in composite materials using PD approach implemented in FEA. A bilinear PD damage law is formulated for PD interface by modifying the original failure formulation of Silling. That bilinear PD interface behavior is inspired by the bilinear damage law of CZM models. A MATLAB code is generated to generate PD interactions and corresponding surface correction factors. Proposed formulation is adapted to FEA code ABAQUS. Using generated MATLAB pre-processing code, PD model of a DCB specimen is generated. In addition, the same problem is solved using Cohesive Zone Modelling approach in ABAQUS. PD and FEA results are compared with results from literature. Obtained results indicate that PD solutions correlate well with CZM and results from literature. © 2018 The Authors.