The very large vessels like VLCC and container ship have been built recently and those vessels have smaller structural strength in comparison with the other convectional skips. As a result the fatigue destruction of upper deck occurs a frequently due to the springing phenomenon at the encountering frequencies. In this study, the hydrodynamic loads are calculated by three-dimensional source distribution method with the translating and pulsating Green function. A ship is longitudinally divided into 23 sections and the added mass, damping and hydrodynamic force of each section is calculated. focusing only on the vertical motion. Stiffness matrix is calculated by the Euler beam theory. The calculation is carried out for Esso Osaka.