We introduce methods for modelling and simulation of anisotropic poroelastic materials in the frequency domain. Starting from the equations formulated by Biot, in their anisotropic form, we derive two different symmetric weak forms together with the boundary conditions that has to be satisfied. We employ a mixed displacement-pressure formulation and solutions are obtained by applying the finite element method to the proposed weak forms. In order to illustrate the use of the finite element method, we highlight some particular aspects related to simulations where poroelastic materials are involved. These include convergence of the discretised solution and boundary conditions at interfaces between poroelastic materials and solids/fluids. Results are given for some selected application examples of foam and plate combinations as well as a poroelastic foam with embedded inclusions.