SUMMARYThis paper concerns the buckling load analysis of laminates with a pre-existing delamination, using the finite element method based on the Mindlin plate theory. To deal with the contact problem in the buckling mode, an effective algorithm is presented. In this method, an iterative updating process based on the first-order sensitivity analysis and the quadratic programming technique is proposed to compute the fictitious forces in contacting areas at first. These fictitious forces are then transferred into the stiffness parameters of some artificial springs. The original stiffness matrix of system can be modified, using these artificial springs. Finally, the penetration between two delaminated layers in the buckling mode can be prevented effectively. Numerical examples show that this method is very efficient to solve the contact problem in eigenvalue analysis from the viewpoint of its accuracy, stability and convergence speed. The effects of contact and delamination size on the buckling load analysis are also investigated.