Numerical simulation of fatigue crack propagation based on fracture mechanics and the conventional finite element method requires a huge amount of computational resources when the cracked structure shows a complicated condition such as the multiple site damage or thermal fatigue. The objective of the present study is to develop a simulation technique for fatigue crack propagation that can be applied to complex situations by employing the continuum damage mechanics (CDM). An anisotropic damage tensor is defined to model a macroscopic fatigue crack. The validity of the present theory is examined by comparing the elastic stress distributions around the crack tip with those obtained by a conventional method. Combined with a nonlinear elasto-plastic constitutive equation, numerical simulations are conducted for low cycle fatigue crack propagation in a plate with one or two cracks. The results show good agreement with the experiments. Finally, propagations of multiply distributed cracks under low cycle fatigue loading are simulated to demonstrate the potential application of the present method.