The paper presents numerical and experimental research on glulam delamination in a double lap connection with predominant shear stresses. Laboratory tests and wide literature survey enabled to determine timber and glue joint parameters. Cohesive zone theory, generally used for epoxy matrix and fiber reinforced composites, was adopted to modelling glue layer delamination in glulam elements. Numerical models were validated with laboratory tests.