The high demand of low cost wind energy needs to design large scale turbine blades with reduced weight which poses great challenges to their structural integrity while prone to extreme wind gusts. The loading can cause large-deflection bending and damage leading to significant drop in the load-bearing ability of long composite wind turbine (WT) blades. In this study, a comprehensive finite element (FE) modelling procedure is developed to simulate structural integrity and damage in composite blade using ANSYS software. The three-dimensional blade model is analyzed by carrying out geometrically nonlinear FE analysis to investigate the blade deformation and highly stressed regions leading to possible failure modes. The results show that the blade suction side is subjected to high compressive stress causing local skin buckling, which is further investigated using linear buckling analysis. Such local buckling drives interfacial debonding between skin and spar joined with a weak adhesive. Subsequently, the interfacial damage in the identified critical region is modelled by developing a damage submodel employing cohesive zone model (CZM) approach at the skin-spar interface. The analysis results indicate that buckling driven skin-spar debonding at adhesive interface is initial damage mode which can lead to progressive failure of the blade structure. Consequently, the ultimate load bearing capacity of WT blade is governed by a coupled buckling and debonding phenomenon even at load level below the ultimate design load. The simulation methodology adopted in this study can be employed to develop reliable and cost-effective computational tools for analyzing structural integrity and assessing damage in blade structure than expensive experimental testing.