In the present work, the doubly curved spherical and cylindrical laminated composite shells are modelled and analysed in the framework of non-polynomial axiomatic approach. The inverse hyperbolic shear deformation theory originally developed for the laminated composite plates is extended to model the deformation characteristics of laminated composite shells. The mathematical formulation is developed under the assumption of linear structural kinematics and linear-elastic-orthotropic material behaviour. The governing equations of the model are obtained using the principle of virtual work and solved in exact manner for simply supported boundary conditions following the Navier solution methodology. The bending response of thick and thin spherical and cylindrical shells subjected to different types of transverse loads such as point load, uniform load and sinusoidal load is analysed in the framework of developed methodology. The obtained results due to inverse hyperbolic shear deformation theory are compared with other shell theories and on the basis of this comparison, the validity and applicability of the inverse hyperbolic shear deformation theory for doubly curved spherical and cylindrical shells is ensured.