Continuum damage mechanics (CDM) [12] provides an effective method to describe the degradation of the material. Based on thermodynamics, damage-coupled elasto-plastic or visco-plastic constitutive models combined with damage evolution laws can be used to model the evolution of ductile plastic damage, fatigue, creep, and creep-fatigue interaction. The details of those methods can be found in the publications of Kachanov [13], Lemaitre and Chaboche [14], Lemaitre and Rodrigue [15], as well as Voyiadjis and Kattan [16]. Moreover, the continuum damage mechanics has been extended to describe the evolution of anisotropic damage in metal matrix composites [17,18]. Besides many theoretical models have been achieved, and CDM has also been widely used in engineering applications [19,20]. Especially, the studies performed by Basaran group have achieved abundant results on the research of fatigue and visco-plastic interaction behavior of solder joints to simulate the fatigue failure prediction [21,22] and to investigate the viscoplastic-fatigue interaction behavior of solder joints [23–25]. Two methods usually are applied to the fatigue life prediction. One is the uncoupled method, which obtains the fatigue life by integrating the damage evolution equation with ignoring the coupling effect between stress field and damage field. The other is the coupled method, which calculates the fatigue life by accumulating the damage cycle by cycle with considering the occurrence of coupling effect in each loading cycle. Generally, for the purpose of accurate applications or when the coupling between strains and damage is strong, the coupled CDM approach is very necessary. In those cases, the constitutive models need to be implemented by the numerical method, and the damage accumulation will be calculated cycle by cycle to take account of the material degradation and stress redistributions [15].