Superelasticity in shape memory alloys is an important feature for actuators and medical devices. However, the function of the devices is typically limited by mechanical bandwidth and fatigue, which are dominated by the microstructures. Thus, in order to correlate the mechanical response and the microstructures, the microstructural evolution in NiTi single crystals under the compression, tensile, and shearing tests is simulated by molecular dynamics (MD) in the current study. Then, the martensite variant identification method, which identifies the crystal variants/phases of each lattice based on the transformation matrix, is used to post-process the MD results. The results with the detailed information of variants and phases reveal many features that have good agreement with those reported in the literature, such as X-interfaces and the transitional orthorhombic phase between the austenite and monoclinic phases. A new twin structure consisting of diamond and wedge-shaped patterns is also discovered. The macroscopic behavior, such as stress-strain curves and the total energy profile, is linked with the distribution of dislocation and twin patterns. The results suggest that the loading cases of shear and compression allow a low critical strain for the onset of martensitic transformation and a better superelasticity behavior. Therefore, the two loading cases are suitable to apply to the NiTi actuators. The current work is expected to provide insight into the mechanical responses and design guideline for NiTi shape memory alloy actuators.