Isoparametric Element Analysis of Two-Dimensional Unsaturated Transient Seepage

Abstract

A FEM for unsaturated transient seepage is established by using a quadrilateral isoparametric element, considering the fact that the main permeability does not coincide with the axis situation. It creates a function by using the element’s node hydraulic head and shape function instead of the real head in the Richard seepage control equation. With the help of the Galerkin weighted residual method, a FEM equation is given for analyzing 2-dimensional transient seepage problem. Further, based on the Jacobi matrix and Gauss numerical integral, it determines the elements of stiffness and capacitance matrices. This FEM equation considers not only the anisotropic of soil but also the uncoincidence between permeability and the axis. It is a common form of transient seepage. In the end, two examples illustrate the node accuracy of the quadrilateral element and the correctness of this FEM equation.