A nonlocal layerwise 3D bending solution for multilayered functionally graded shallow nanoshells is presented. Eringen’s nonlocal elasticity theory is used for approaching the small length scale effect. Curvilinear orthogonal coordinate system is employed for writing the constitutive and equilibrium equations. The middle-plane of displacements are assumed as a Navier solution which are only valid for simply supported spherical, cylindrical panels and rectangular plates. The thickness domain is discretized by the well-known Chebyshev–Gauss–Lobatto, and its derivatives are calculated numerically by the Differential Quadrature method (DQM). Lagrange interpolation polynomials are utilized as the basis functions for DQM. The traction conditions at the top and the bottom of the shell panel and continuity conditions of the interlaminar adjacent layers for transverse stresses are considered, so this method presents layerwise capabilities. The results are compared with higher order theories proposed in the literature, and due to the 3D nature of the presented formulation, benchmark solutions are provided.