The study of the present paper is inspired by Gangl, Kaneko and Zagier's result of the connection with double zeta values and modular forms.
We introduce double Eisenstein series $E_{r,s}$ in positive characteristic with double zeta values $\zeta_A(r,s)$ as their constant term and compute the t-expansions of the double Eisenstein series. Moreover, we derive the shuffle relations of double Eisenstein series which match the shuffle relations of double zeta values in \cite{Huei2015}.