In this paper, we study a finite element scheme of some semi-linear elliptic boundary value problems in high-dimensional space. With uniform mesh, we find that, the numerical scheme derived from finite element method can keep a preserved multisymplectic structure.