A unified hybrid analytical-numerical method is proposed for vibration analysis of stiffened hulls coupled with non-axisymmetric internal structures. The coupled system is firstly decomposed into the submarine hull and internal structures. Then, the stiffened hull is further divided into ring stiffeners, bulkheads and shell segments, including conical and cylindrical shell segments. The shell strips and axisymmetric reinforcements are uniformly described by the Flügge formula and power series method. For non-axisymmetric interior structures, finite element method is utilized to establish the motion equations owing to their irregular shapes. To synthesize the shell and non-axisymmetric internal structures, the receptance functions at coupling nodes are respectively deduced by condensing the global stiff matrices of the hull and interior substructures. A precise coupling method fully considering all 6 DOFs at coupling nodes, rather than just 3 or 4 DOFs in previous literature, is proposed to synthesize the receptance functions. To verify the validity of this approach, a stiffened conical-cylindrical-spherical shell coupled with two foundations is manufactured and the forced vibration characteristics of this coupled model subjected to multiple types of excitations are tested. Results obtained by semi-analytical method are simultaneously compared with the ones measured from this experiment and calculated through finite element method (FEM). Effects of the inner structures and different external excitations on forced vibrations are discussed by experiment. Furthermore, how non-axisymmetric internal structures modify wave propagation and conversion characteristics of axisymmetric hulls is investigated through the present method.