Abstract: The interfacial debonding shear strength and ultimate strength of short-cut shape memory alloy (SMA) reinforced epoxy composites were determined based on SMA single fiber pull-out test and uniaxial tensile test, respectively. An algorithm for the automatic generation of unidirectional random distribution short-cut SMA fiber/epoxy composites was developed, using Monte-Carlo method and boundary condition control equation via algorithmic processor description language (APDL). The interfacial debonding process between SMA fibers and epoxy matrix was simulated by means of exponential cohesive zone model (CZM). The results indicate that the macroscopic elastic modulus gradually reduce with the decrease of fiber slendevness ratio and the overall elastic modulus of composite are enhanced owing to the driving of SMA fiber, under the same fiber volume.