Abstract. In this paper, we construct extrapolated expanded mixed nite element approximations to approximate the scalar unknown, itsgradient and its ux of semilinear Sobolev equations. To avoid the di-culty of solving the system of nonlinear equations, we use an extrapolatedtechnique in our construction of the approximations. Some numerical ex-amples are used to show the eciency of our schemes. 1. IntroductionIn this paper, we consider the following semilinear Sobolev equation(1.1) u t (u x + u tx ) x = f(x;t;u); in (0;L) (0;T];(1.2) u x (0;t) + u tx (0;t) = u x (L;t) + u tx (L;t) = 0; in (0;T];(1.3) u(x;0) = u 0 (x); in (0;L);where 0 <T<1and u 0 (x) and f(x;t;u) are suciently smooth functions sothat (1.1)-(1.3) has a suciently smooth unique solution. The problem (1.1)-(1.3) arises in the various areas, for examples, in the ow of uids through ssured materials [1] and thermodynamics [4]. For details about the physicalsigni cance and the existence and uniqueness of the solutions of the Sobolevequations, we refer to [1, 2, 3, 4, 6, 7, 12].The advantages of mixed nite element formulations can be summarizedas follows: one can simultaneously approximate both the displacement and theux(or the stress or the pressure) and one can approximate the ux to the sameorder of convergence as the unknown scalar u(x;t). Due to these advantages,the mixed nite element method has been applied to some types of Sobolevequations by the authors in [8, 10, 11].