Meshless local collocation method produces much better conditioned matrices than meshless global collocation methods. In this paper, the meshless local collocation method based on thin plate spline radial basis function and first-order shear deformation theory are used to calculate the natural frequencies and mode shapes of laminated composite shells. Through numerical experiments, the accuracy and efficiency of present method are demonstrated.