We present a computational framework to explore the effect of microstructure and constituent properties upon fracture toughness fibre-reinforced polymer composites. To capture microscopic matrix cracking fibre-matrix debonding, couples phase field method cohesive zone model in context finite element method. Virtual single-notched three point bending tests on fibre reinforced composites are cond...