Abstract: Based on the zig-zag deformation model, the first order shear strain effect laminated plate theory and the von Karman non-linearity assumption, a finite element method for fractural analysis of the debonded sandwich composites is developed. The total energy release rate can be estimated by a virtual crack closure technique (VCCT). The delamination fracture behavior of the sandwich plates containing an different shape interfacial crack under compression is studied. From some typical numerical examples and discussions, it can be seen that the magnitude and distribution of the energy release rate along the crack front are dependent on the crack area, shape and loading orientation.