“…Motivations of the different responses in the initial elastic branch between the numerical solutions and the analytical one are well known in literature (see References 16 and 56) and are due to the two different delamination conditions, based on the linear elastic fracture mechanics theory for the analytical solution, with an ideally brittle traction‐separation law, whereas in the CZMs, the delamination condition produces a bilinear response with a smooth transition from the elastic behavior to the fully debonded one. The maps of stress , and , obtained by the HEE with embedded interface at four different loading conditions of imposed displacement u = 5, 10, 15, 20mm, are plotted in the deformed configuration in Figures 10 to 12.…”