Summary: In this article, we are concerned with the interactions of delta shock waves with contact discontinuities for the relativistic Euler equations for Chaplygin gas by using split delta functions method. The solutions are obtained constructively and globally when the initial data consists of three piecewise constant states. The global structure and large time-asymptotic behaviors of the solutions are analyzed case by case. During the process of the interaction, the strengths of delta shock waves are computed completely. Moreover, it can be found that the Riemann solutions are stable for such small perturbations with special initial data by letting perturbed parameter \(\epsilon\) tends to zero.