For computation of large amplitude motions of ships fastened to a dock, a fast evaluation scheme is implemented for computation of the time-domain Green function for finite water depth. Based on accurate evaluation of the Green function directly, a fast approximation method for the Green function is developed by use of Chebyshev polynomials. Examinations are carried out of the accuracy of the Green function and its derivatives from the scheme. It is shown that when an appropriate number of polynomial terms are used, very accurate approximation can be obtained.