The two-dimensional (2D) delamination growth in fiber-reinforced polymer (FRP) laminates with in-plane isotropy under Mode I loading condition was numerically investigated using finite element analyses. Two sizes of plate models were developed, focusing on different fracture stages. Cohesive elements were employed to simulate the fracture behavior in the presence of large-scale bridging (LSB). The influences of the pre-crack shape/area, loading zone shape/area and fracture resistance were parametrically studied. It was found that either a flatter pre-crack shape or a flatter loading zone shape could result in higher initial structural stiffness and less uniform distribution of the strain energy release rate (SERR) along the pre-crack perimeter during crack initiation and early propagation. However, they had only a minor effect on the stiffness after full fiber bridging development in all directions. The plates finally achieved constant stiffness, which increased linearly with the fracture resistance. The final crack shape was dependent on the loading zone shape and area, but the effects were relatively weak.