Based on the 3D continuous shell theory and the 1st-order shear deformation theory,a moving least squares meshless method was proposed for solving the linear bending problem of functionally graded sandwich shells with the layered method.With the mapping technology,the 2D meshless node information on the convec-ted coordinate system was mapped on the 3D shell,and the shape function of moving least squares(MLS)ap-proximation was formed on the convected coordinate system.Due to the inability of shell numerical solutions based on the 3D continuous shell theory to provide an explicit expression for the thickness direction of a specif-ic shell,the portion of material parameter changes in the functionally graded sandwich material shell structure was divided into several layers,and the material parameters of each layer were obtained as constants.The gov-erning meshless equation for linear bending of functionally graded sandwich shells was derived under the princi-ple of minimum potential energy.Through introduction of a linear transformation in the thickness direction,the Gaussian integral in the thickness direction of each layer was bounded within the range of-1 to 1,without vio-lation of the 1st-order shear deformation theory(FSDT).The essential boundary conditions were employed with the complete transformation method.Finally,the effects of different gradient coefficients,diameter to thickness ratios,and curvature radii on numerical results were discussed through examples of functionally gra-ded sandwich plates,cylindrical shells,and hyperbolic shallow shells with classical geometric shapes.The cal-culated results were compared with the literature solutions.The results show that,the proposed method has the characteristics of good convergence and high computation accuracy in solving linear bending problems of func-tional graded sandwich shells with different shapes.