Purpose The purpose of this paper is to study the mechanical property of three-dimensional (3D) Printed photopolymer (Vero Yellow and Tango Black) with different constant strain rate. According to the experimental results, three constitutive models are used to describe the stress-strain and stress-time relation in the tension and stress relaxation process. Design/methodology/approach The Stratasys Objet 260 was used to prepare the four groups of samples with different photopolymers (Vero Yellow and Tango Black). The stress-strain and stress-time relations are obtained by the uniaxial tensile tests and relaxation tests performed at room temperature with different constant strain rates. The generalized Kelvin model (GKM), standard linearized model (SLM) and fractional order model (FOM) are used to describe experimental data by means of the curve-fitting approach. Findings Experimental results show that the tension stress increases faster at a higher strain-rate for tensile tests. Relaxation stress is influenced by the preload strain-rate for relaxation tests. For the theoretical fitting, the error comparison between three constitutive models and experimental data are calculated to demonstrate the high accuracy in describing the stress-strain relationship for tension. For stress relaxation, the error comparison confirms higher accuracy of FOM with the largest error within 3%, while the error of GKM and SLM up to 10%. Originality/value The paper confirms the viscous-elastic mechanical property of 3D printed photopolymer composites (Vero Yellow and Tango Black) for Stratasys PolyJet. As FOM shows high accuracy both in describing stress-strain and stress-time relation for tension and stress relaxation process, it can be directly used as a constitutive model to predict mechanical properties for engineering application.