In this paper, growths of different types of failures including adhesion, cohesion and delamination for a single lap joint (SLJ) of composite laminates were investigated using three-dimensional geometrically non-linear finite element analysis and by adopting a suitable modeling technique. A unique damage modeling method called sub-laminate modeling was employed for the modeling of damages of different failure modes so as to avoid the oscillatory stress and displacement fields around the damage front. The strain energy release rate (SERR) parameter was used for studying the damage growth and the individual and total components of the SERR along the various damage fronts are evaluated using the virtual crack closure technique (VCCT) based on the linear elastic fracture mechanics (LEFM) approach. This study reveals: that the opening mode is the dominant mode of the propagation for the adhesion and delamination damages, while the sliding mode is dominant for the cohesion failure; that the cohesion failure grows at a faster rate than the adhesion failure; and that the delamination front entrapped within the overlap region in the top adherend of the SLJ grows faster when the delamination damages are present simultaneously in both the adherends. This is particularly true when the delamination centers are exactly aligned with the overlap ends of the joint. (C) Koninklijke Brill NV, Leiden, 2011