A new crack opening/closing simulation method is developed utilizing the crack tip stress field parameters evaluated by FEM. By combining this method and the CP-System, the growth behavior of fatigue cracks can be computed for large-scale structures under variable amplitude loading. In order to demonstrate the usefulness of the developed program, fatigue crack propagation in a ship structure is computed under realistic wave load sequence.