A fast boundary element method (FBEM) based on multilevel fast multipole algorithm (MLFMA) is introduced for simulating acoustic emission from large vibrating structures. The vibration of a structure is first simulated using a finite element method (FEM). The surface velocity distribution, extracted from the FEM result, is then used a source term in FBEM. The fast acoustic BEM solver utilizes a special broadband version of MLFMA that allows efficient and accurate modeling of large scale acoustics problems on a broad range of frequencies.