To address the issues of slow convergence and poor accuracy in existing dynamics methods caused by non-uniform variable cross-sections and layered material stiffness variations in variable stiffness acoustic black hole structure (VSABH). A multiscale Gaussian expansion method (MSGEM) is proposed in this paper. The structural dimensions are taken into account, and the structural parameters of multiple shape functions are adaptively selected based on stiffness and cross-sectional variation parameters. This results in the formation of shape function groups of various scales, which in turn creates the matrix of multi-scale Gaussian function groups. The feature information of the displacement field is realized to be extracted from multiple scales, and the fitting error of MSGEM to the displacement field of VSABH beam is within 2.02%. The missing solution of eigenfrequency of Gaussian expansion method is avoided, and the modal vibration patterns are basically matched, which verifies the validity of MSGEM. Meanwhile, the advantages of low-frequency multimodal vibration reduction of VSABH beams are highlighted from several aspects, and the trend of the adjustment of the vibration reduction characteristics of VSABH beams by different parameters is clarified, which provides valuable design guidance for ABH metamaterials oriented to low-frequency vibration reduction.