The simulation of shipping of green water and prediction of slamming loads are fundamental issues in the field of naval architecture. However, severe green water and slamming are rarely researched in freaking waves due to the complexity associated with strong nonlinear interactions between waves and ships. In this study, the numerical method was applied for the simulation of green water and slamming on ships advancing in freaking waves. Firstly, the numerical simulation of freaking waves, green water and ship motions were compared with theoretical and experimental data respectively to verify the reliability and accuracy of the numerical method in this study. Secondly, a three-dimensional numerical tank was built to investigate the phenomenon of multiple green water and slamming caused by freaking waves in the time domain, which indicated the second green water had a longer duration and more violent intensity compared to the first. Finally, the FFT method was adopted to analyze the features of motions and slamming loads in spectral. The results indicated that the frequencies of peak pressure in the bow and deck area were affected significantly by high frequency components of the freaking waves, while peak pressure in the superstructure area was mainly affected by low frequency.