Responses of marine vehicles and ocean structures in a seaway can be predicted by applying the probabilistic approach. When we consider a linear system, the responses in a random seaway can be evaluated through spectral analysis in the frequency domain. But when we treat nonlinear system in irregular waves, it is necessary to get time history of waves. In the previous study we introduced one-directional waves (long crested waves)as wave environment and carried out calculations and experiments in the waves. But the real sea in which marine vehicles and structures are operated has multi-directional waves (short crested waves). It is important to get a simulated random sea and analyse dynamic problems in the sea. We need entire sample function or probabillty density function to infer statistical value of random process. However if the process are ergodic process, we can get statistical values by analysis of one sample function. In this paper, we developed the simulation technique of multi-directional waves and proved that the time history given by this method keep ergodic characteristics by the statistical analysis.