A three-dimensional constitutive model for shape memory alloys (SMAs) is developed along the lines of the Stebner–Brinson (SB) implementation of the Panico–Brinson model. Plastic kinematic hardening behavior is simulated in addition to elastic deformation and phase transformation. A series of finite element simulations is carried out using this model to investigate the localization effects of the stress and strain field on NiTi plates with structured arrays of pores. The application of this model on porous architectures provides insight into how geometric features influence the mechanics of the structure. The incorporation of plastic deformation shows a marked decrease in the maximum stress levels; these results are more consistent with experimental data as compared to the original SB model. Furthermore, the new results demonstrate that clustered pores lead to more distributed stresses and transformation compared to a dispersed configuration of pores, indicating the importance of pore geometry in determining the stress and strain distribution. The improved model provides a practical tool toward design and optimization of porous SMA structures.