We formulate a theory of the aging phenomena in shape memory alloys. The order parameter of the martensitic transformation is coupled with a secondary variable which follows the evolution of the order parameter but is assumed to be an extremely slow variable. From this set of equations, the stress-strain relation due to the motion of twin boundaries is obtained by taking account of the aging effect. The results are found to be consistent, at least qualitatively, with experiments.