The paper deals with numerical analysis of double cantilever beam (DCB) predefined to Mode I Interlaminar Fracture Tests of GRFP unidirectional laminates. The numerical analyses were performed in the ANSYS program based on the finite element. In geometrically nonlinear analysis, two algorithms, responsible for initiation and propagation of delamination front, were applied: Virtual Crack Closure Technique (VCCT) and Cohesive zone method (CZM). Due to the unidirectional arrangement of layers of the laminate, the problem of DCB test was solved with the use of one- and three-dimensional models with the implementation of linear interface element and contact element. The present study highlights the limitations of existing formulae used to reliably reflect the behavior of DCB. The use of three-dimensional models allowed confirming the curved shape of the delamination front observed in experimental studies. The application of the VCCT in the three-dimensional model led to an underestimation of the global response (force-opening displacement curve) recorded during numerical DCB test.