The finite element modeling is used to study the buckling and postbuckling behavior of composite laminates with an embedded delamination. Degenerated shell element and rigid beam element are utilized for the finite element modeling. In the nonlinear finite element formulation, the updated Lagrangian description method based on the second Piola-Kirchhoff stress tensor and the Green strain tensor is used. The buckling and postbuckling behavior of composite laminates with a delamination are investigated for various delamination sizes, stacking sequences, and boundary conditions. It is shown that the buckling load and postbuckling behavior of composite laminates depend on the buckling model which is determined by the delamination size, stacking sequence and boundary condition. Also, results show that introduction of couplings can reduce greatly the buckling load.