The aim of this research is to determine the relaxation and creep modulus of 3D printed materials, and the numerical research is based on the finite volume method. The basic material for determining these characteristics is ABS (acrylonitrile butadiene styrene) plastic as one of the most widely used polymeric materials in 3D printing. The experimental method for determining the relaxation functions involved the use of a creep test, in which a constant increase of the stress of the material was performed over time to a certain predetermined value. In addition to this test, DMA (dynamic mechanical analysis) analysis was used. Determination of unknown parameters of relaxation functions in analytical form was performed on the basis of the expression for the storage modulus in the frequency domain. The influence of temperature on the values of the relaxation modulus is considered through the determination of the shift factor. Shift factor is determined on the basis of a series of tests of the relaxation function at different constant temperatures. The shift factor is presented in the form of the WLF (Williams-Landel-Ferry) equation. After obtaining such experimentally determined viscoelastic characteristics with analytical expressions for relaxation modulus and shift factors, numerical analysis can be performed. For this numerical analysis, a mathematical model with an incremental approach was used, as developed in earlier works although with a certain modification. In the experimental analysis, the analytical expression for relaxation modulus in the form of the Prony series is used, and since it is the sum of exponential functions, this enables the derivation of a recursive algorithm for stress calculation. Numerical analysis was performed on several test cases and the results were compared with the results of the experiment and available analytical solutions. A good agreement was obtained between the results of the numerical simulation and the results of the experiment and analytical solutions.