Simple Finite Element models for soil dynamics and earthquake engineering problems in the frequency domain are a fast and valuable tool providing a first approximation before a full non-linear analysis in the time domain is performed.Quite often the problem concerns saturated soils with very small permeability and pore fluid of neglectable compressibility. In the limit, the permeability is assumed to be zero and the pore fluid incompressible. Here, engineers use standard finite element codes formulated in terms of displacements but incompressibility may result in volumetric locking of the mesh with a severe loss of accuracy.The purpose of this paper is to present a simple mixed finite element formulation in the frequency domain based on displacements and pore pressures as main variables. A suitable stabilization technique allowing for equal order interpolation of displacements and pressures has been introduced for incompressible and zero permeability limits.Of course, the range of application is limited to those problems in which the behaviour of the material can be approximated by linear models, and therefore modelling of phenomena such as liquefaction, cyclic mobility or cavitation occur is excluded.The paper shows as well an extremely simple way of coupling solid and water domains as it occurs for instance in quay walls under dynamic loading.