SUMMARYA computational strategy for the evaluation wf delamination stresses in multilayered fibre-reinforced composites is presented. The approach uses special inter-laminar shear elements and an iterative procedure in the calculation of out-of-plane stresses. Both regular and singular elements for the free edge are presented. The solutions are carried on by switching between direct matrix solution methods of small systems of equations and an iterative conjugate gradient method. The proposed discretization avoids ill-conditioning problems, modelling difficulties and large computational cost associated with threedimensional analysis of delamination of holes, cut-outs, joints and areas of complex stress states.