NiTiHf alloys are high temperature Shape Memory Alloys (SMAs) suitable for actuation applications in a wide range of temperatures. Aging heat treatments in NiTiHf SMAs cause the formation of nano-sized, non-transforming and elastic precipitates, which modify phase transformation behavior. Experimental observations of the aged alloys showed a counterintuitive trend of increasing phase transformation strain with precipitation. In order to predict such microstructure effects on constitutive response of precipitated NiTiHf SMAs, in this work a finite element based micromechanical modeling framework is developed. Two types of representative volume elements are considered:(i) exact 3-D reconstructions of precipitates using transmission electron microscopy, and (ii) less expensive ellipsoidal representation of precipitates. The composition changes in the material are modeled, and their effects are incorporated by considering composition dependency of phase transformation properties in the matrix. The resulting model is used to predict the constitutive responses for different aging conditions in Ni 50.3 Ti 29.7 Hf 20 and Ni 50.3 Ti 34.7 Hf 15 SMAs. The model predictions are in good agreement with the experimental measurements, and experimental and modeling results are presented with confidence intervals from sample-to-sample variations. The present work constitutes the first micromechanical modeling of high temperature SMAs with consideration of real microstructures.