To design very large ships, such as very large drillships, we have to estimate the hydroelastic responses of the very large ships in waves. A numerical procedure is described for estimating the hydroelastic responses of very large ships advancing with slow speed in waves. The developed numerical approach is based on a combination of the three-dimensional source distribution method and the finite element method, including fluid-structure interaction by regarding a very large ship as many hull elements connected with elastic beam elements. Numerical results are compared with experimental and numerical ones obtained in the literature. The results of comparison confirmed the validity of the proposed approach.