IN THIS PAPER, THE FINITE ELEMENT METHOD (FEM) IS USED TO STUDY THE DELAMINATION GROWTH AND POST BUCKLING OF COMPOSITE LAMINATES CONTAINING EMBEDDED DELAMINATION UNDER THE COMPRESSIVE LOADING. DELAMINATION DECREASES AN ELASTIC BUCKLING LOAD OF THE COMPOSITE STRUCTURES AND LEADS UNEXPECTED STRUCTURAL FAILURE AT LOADS BELOW THE DESIGN LEVEL. THE INVESTIGATED COMPOSITE LAMINATES WHICH ARE MADE OF CARBON FIBER/EPOXY HAVE THE STACKING SEQUENCE OF [90/ (0/90) 3// (0/90) 14] CONTAINING A CIRCULAR AS WELL AS SQUARE DELAMINATION. THE SYMBOL // ILLUSTRATES THE POSITION OF DELAMINATION IN THE LAMINATE. FOR MODELING AND ANALYSIS THE LAMINATE, THE THREE-DIMENSIONAL CONTINUUM SHELL ELEMENTS HAVE BEEN USED VIA THE ABAQUS 6.12 SOFTWARE. FURTHERMORE, THE SURFACEBASED COHESIVE BEHAVIOR IS UTILIZED TO SIMULATE COHESIVE ZONE MODEL (CZM) WHICH IS APPLIED FOR DAMAGE PROPAGATION. IN ADDITION, THE SURFACE TO SURFACE CONTACT HAS BEEN APPLIED IN THE DELAMINATION ZONE TO PREVENT OVERLAPS BETWEEN ELEMENTS THROUGH SIMULATION. TWO DIFFERENT APPROACHES HAVE BEEN ADOPTED TO SOLVE BUCKLING PROBLEM: LINEAR, BASED ON THE SOLUTION OF EIGENVALUES PROBLEM, AND NONLINEAR. THE RIKS METHOD WAS USED TO SOLVE THE NONLINEAR SYSTEM. THE LOAD/DISPLACEMENT BEHAVIOR OF EMBEDDED DELAMINATION GROWTH FOR COMPOSITE LAMINATES IS PRESENTED. RESULTS REPORTED IN THIS STUDY ARE FOCUSED ON THE EFFECTS OF THE SHAPES AND GEOMETRY IMPERFECTIONS ON POST BUCKLING RESPONSE AS WELL AS. VALIDATION OF THE CIRCULAR DELAMINATION SPECIMEN HAS BEEN PERFORMED WITH EXPERIMENTAL AND PREVIOUS NUMERICAL DATA.