Based on the fluid saturated biphasic porous media model deduced from mixture theory,a finite element formulation with u S u F p variables for dynamic response analysis is given out.This method overcomes the difficulty of choosing suitable penalty parameter in penalty finite element method,and the accuracy of pressure distribution obtained with the mixed finite element method is higher than the penalty finile element method.An iterative solution method is suggested to solve the system of equations whose coefficient matrices are non positive definite.It is concluded from patch test that the order of interpolation function for pressure variables must be higher than that of displacement variables of solid phase.