In this research the coupling problem between the elastic structure and the fluid, specially the hydroelastic harmonic vibration problem, is studied. In order to couple the domains, i.e., the structural domain and the fluid domain, the boundary integral method(direct boundary integral formulation) is used in the fluid domain in combination with the finite element method for the structure. The boundary integral method has been widely developed to apply it to the hydroelastic vibration problem. The hybrid boundary integral method using eigenfunctions on the radiation boundaries and the boundary integral method using the series form image-functions to replace the even bottom and free surface boundaries in case of high frequencies have been developed and tested. According to the boundary conditions and the frequency ranges the different boundary integral methods with the different idealizations of the fluid boundaries have been studied. Using the same interpolation functions for the pressure distribution and the displacement the two domains have been coupled and using Hamilton principle the solution of the hydroelastic have been obtained through the direct minimizing process. It has become evident that the finite-boundary element method combining with the eigenfunction or the image-function method give good results in comparison with the experimental ones and the other numerical results by the finite element method.