Abstract. In this paper, we construct extrapolated expanded mixed finite element approximations to approximate the scalar unknown, its gradient and its flux of semilinear Sobolev equations. To avoid the difficulty of solving the system of nonlinear equations, we use an extrapolated technique in our construction of the approximations. Some numerical examples are used to show the efficiency of our schemes.