Emergent universe (EU) cosmological models with viscosity in a modified gravity which contains a general function f(R,T)f(R,T), where RR and TT denote the curvature scalar and the trace of the energy–momentum tensor, respectively, are studied in a flat Friedmann–Robertson–Walker metric. Cosmological solutions are obtained in f(R,T)f(R,T) theory of gravity, which represented as f(R,T)=f(R)+f(T)f(R,T)=f(R)+f(T) with bulk viscosity that described by Eckart theory, truncated Israel Stewart theory and full Israel Stewart Theory. The physical and geometrical features of the EU models in f(R,T)=α1R+λ1Tf(R,T)=α1R+λ1T gravity, where α1α1 and λ1λ1 are coupling parameters, with bulk viscosity are studied in details. Constraints of the EU models parameters in f(R,T)f(R,T) gravity with bulk viscosity are estimated from observational data set. The stability analysis of the equilibrium points admitting cosmological solutions of the dynamical system associated with the evolution in the modified gravity is studied in Eckart theory, truncated Israel Stewart theory and full Israel Stewart theory.