A numerical procedure is described for estimating the effects of the multi-directional irregular waves on the motion responses and tension variations of the ISSC-TLP. The numerical approach is based on a three-dimensional source distribution method and a spectral analysis technique of directional waves. The spectral description for the linear system of ISSC-TLP in the frequency domain is sufficient to completely define the motion responses and tension variations. This is because both the wave inputs and responses are stationary Gaussian random processes, of which the statistical properties in the amplitude domain are well known. The numerical results for the linear motion responses and tension variations in regular waves are compared with the experimental and numerical ones, which are obtained in the literature. The results of comparison confirmed the validity of the proposed approach.