In this paper an efficient procedure for obtaining a cohesive law for Mode I timber fracture (crack opening), based on the Double Cantilever Beam (DCB) tests is given. DCB tests were performed on ten European spruce specimens in order to determine the energy release rate vs crack length (R curves). Two crucial parameters - crack length during the experiment and the crack tip opening displacement were obtained using 2D Digital Image Correlation (DIC) technique. In order to determine accurate fracture resistance (R curve), procedure which includes calculating cumulative released energy was employed. The cohesive law for Mode I fracture of wood was obtained by differentiation of the strain energy release rate as a function of the crack tip opening displacement. This cohesive law is further implemented in the successful numerical modelling of failure modes in large-scale end-notched glulam beams which were experimentally tested in four-point bending configuration.