A general fatigue life prediction methodology, based on a unified creep plasticity damage (UCPD) model, was developed for predicting fatigue cracks in 95.5Sn–3.9Ag–0.6Cu (wt %) solder interconnects. The methodology was developed from isothermal fatigue tests using a double-lap-shear specimen. Finite element analysis model geometries, mesh densities, and assumptions were detailed for both a full model (an octant-symmetry slice of the entire ball grid array (BGA) assembly) and a submodel (the solder joint deemed most likely to fail and the surrounding package layers) to facilitate fatigue prediction. Model validation was based on the thermal mechanical fatigue of plastic BGA solder joints (250–4000 thermal cycles, , and ). Metallographic cross sections were used to quantitatively measure crack development. The methodology generally underpredicted the crack lengths but, nonetheless, captured the measured crack lengths within a error band. Possible shortcomings in the methodology, including inaccurate materials properties and part geometries, as well as computational techniques, are discussed in terms of improving both the UCPD constitutive model and the fatigue life prediction methodology fidelity and decreasing the solution time.