A new method of vibration analysis of continuous bodies is presented. In the method, the domain to be analyzed is divided into several subdomains. Analytical solutions are used for subdomains which have regular shapes and the FEM is used for irregularly shaped subdomains. The analytical solutions are transformed into relations between generalized forces and generalized displacements at nodal points located on the boundaries of the regular subdomains. This transformation is based on a kind of the energy method and the concentration of stresses. The relations are formally the same as those of the FEM. Therefore, the solutions for the regular subdomains can be joined directly to the FEM solutions for the irregular subdomains. As simple examples, some two-dimensional acoustic problems are analyzed by this method. It is confirmed that the method brings more accurate results and shorter computational time in comparison with the FEM only.