A numerical study is undertaken to simulate failure of solder joint caused by cyclic shear deformation. A progressive ductile damage model is incorporated into the rate-dependent elastic-viscoplastic finite element analysis, resulting in the capability of simulating damage evolution and eventual failure through crack formation. It is demonstrated that quantitative information of fatigue life, as well as the temporal and spatial evolution of fatigue cracks, can be explicitly obtained.