The work hardening of alloys hardened by precipitate heat treatments depends on the distribution of the precipitate sizes and the solute level left in the metal matrix. A mean field theory for precipitation is first applied for the ageing and subsequently it is coupled to a work hardening model to study the stress-strain responses of age hardened conditions of AA6xxx alloys. The predictions are compared to mechanical experiments and to TEM characterisations.