The well-known tension/compression asymmetry exhibited by nearly equiatomic NiTi has been previously modeled using a hardening potential for compression and a partially softening one for tension to represent pseudoelastic phase transformations (Jiang et al., 2016b). The present study first extends this constitutive model to include anisotropy revealed in the combined axial force-internal pressure experiments on NiTi tubes of Bechle and Kyriakides (2016a). The model is then calibrated anew, implemented in a finite element analysis of tubes and used to simulate the entire range of biaxial experiments performed. Overall, the simulations reproduce well the stress-average strain hystereses and the transformation stress loci, while for hoop dominant stress paths the extents of the transformation strains are somewhat over-predicted. The evolution of localization in the form of high or low strain helical bands, the variation of helix angles with respect to the stress ratio, and the dissipated energy compare favorably. The hardening response and essentially homogeneous deformation exhibited in the neighborhood of the equibiaxial stress state is reproduced, but with reduced hardening and mild inhomogeneity. Despite some minor differences, the results demonstrate the overall success of the analysis in reproducing the phenomena-rich behavior exhibited by tubular NiTi structures under biaxial loadings.