This analytical paper investigates application of a novel higher-order shear deformable mathematical model for investigating vibrational analysis of a doubly curved shell in which the kinematic relations are extended using the variable-thickness transverse deflection. The Hamiltonian method is extended to derive the governing motions equations in the curvilinear coordinate system. The constitutive relations are extended using the experimental and statistical relations in the literature for a copper matrix including various amount of three-dimensional graphene nanofillers named as graphene origami using hydrogenation process. The analytical results are presented using the mathematical model in order to trace the impact of graphene origami characteristics and geometric parameters of the double-curved shell. The temperature-dependency relations for material properties as well as folding and content are employed for investigating the impact of affecting parameters. The results are presented with and without thickness stretching effect to arrive at a more confidence results.