Shape memory alloys (SMAs) are a subset of broad class of smart materials which exhibit unique properties of shape memory effect and pseudoelasticity/superelasticity. Shape memory effect is the capability to recover apparently permanent strains upon heating, while superelasticity is the ability to recover very large strains (up to 8%) during a mechanical loading/unloading cycle. The unique properties of shape memory and superelasticity are due to a solid state martensitic transformation from a high symmetry, high temperature phase known as austenite to a low temperature, low symmetry martensite phase. These recently-revealed materials also benefit from other properties such as high force-to-weight ratio, biocompatibility, silent response, smooth and lifelike operations which cause them to gain an immense attention in many real-life engineering applications from biomedical to aerospace fields. Although they can be manufactured in different shapes, SMAs are mainly used in the form of wire or ribbon. Therefore, 1-D simulation of SMAs is usually required in the investigation of real-life problems. Many constitutive models have been presented for SMAs by researchers since the last three decades, and different numerical approaches for the implementation of these models have been so far proposed. Most of the existing solutions are only applicable for SMA wires lonely or for particular simple types of structures consisting SMA wires. Since smart structures consisting SMA wires are various and complicated, these special case studies are no longer effective in such cases. Beside these case studies, there are also a few general works showing several drawbacks. In the present study, in order to provide a general solution, a common 1-D constitutive model for SMAs capturing both pseudoelasticity and shape memory effect was implemented into ABAQUS commercial finite element package via a user material subroutine (UMAT). One of the main benefits of this 1-D thermomechanical UMAT is its capability to investigate any complex mechanical as well as thermal loading path. Hence, it can be utilized in finite element simulation of any complicated smart structure consisting of SMA wires. Moreover, this proposed UMAT can be employed in studying the static, coupled thermomechanical (dynamic), and coupled electrothermomechanical behavior of SMAs. In order to validate the proposed numerical approach, its predictions were shown to be in a good agreement with those obtained from analytical, experimental and other approved numerical solutions under different thermomechanical loading and boundary conditions. The effect of heating rate, conductive and convective coefficients on the actuation response of SMA actuators in static condition was investigated. Also it was shown that when mechanical loadings are applied relatively fast, the isothermal assumption would be no longer valid and as loading rate increases, thermal effects can lead to very different transformation characteristics. In this case, the temperature variation is dependent on factors such as strain rate, convective coefficient, and wire diameter. Moreover, the influence of the latent heat of transformation on the performance of SMA actuators was numerically explored and it was demonstrated that the latent heat will decrease the resultant actuation due to reduced martensitic transformation. Finally the electrothermomechanical behavior and the effective parameters in this case were studied. It was shown that if the electrical resistivity assumed to be constant and equal to its martensitic value, the resulting temperature rise will be slower. Keywords SMA wire, electrothermomechacial behavior, Finite Element simulation, ABAQUS, subroutine