Abstract: In this work, Brinson's one-dimensional constitutive relation and the refined kinetic relations for phase transformation of shape memory alloy (SMA) are used to write UMAT subroutine, so that SMA becomes one member of the material library of ABAQUS. Firstly, the stress-strain relation at various temperatures, the temperature-stress relation under constrained condition and the stress-strain relation under cyclic loading of a SMA bar are calculated to verify the UMAT code. Then the thermoviscoelastic property of the beam matrix is introduced into UTRS subroutine. The bending of the SMA reinforced epoxy resin beam under lateral load is simulated by FEM code of ABAQUS and the results coincide well with experiment.