thermo‐hydro‐mechanical process in porous media involving pollutant transport, Arch. Comput. Methods Eng., 2, 1–54.
51 Schrefler, B. A., (2002). Mechanics and thermodynamics of saturated‐unsaturated porous materials and quantitative solutions, Appl. Mech. Rev., 55(4), 351–388.
52 Schrefler, B. A. and Gawin D. (1996). The effective stress principle: incremental or finite form?, Int. J. Numer. Anal. Methods Geomech., 20, 785–814.
53 Schrefler, B. A. and Scotta, R. (2001). A fully coupled dynamic model for two phase fluid flow in deformable porous media, Comput. Methods Appl. Mech. Eng., 190, 3223–3246.
54 Schrefler, B. A. and Zhan, X. (1993). A fully coupled model for waterflow and airflow in deformable porous media, Water Resour. Res., 29, 1, 155–167.
55 Serafim, J. L. (1954). A subpressëo nos barreyens – Publ. 55, Laboratorio Nacional de Engenheria Civil, Lisbon.
56 Sheng, D., Sloan, S. W., Gens, A. (2004). A constitutive model for unsaturated soils: thermomechanical and computational aspects. Comput. Mech., 33 (6), 453–465.
57 Simoni, L. and Schrefler, B. A. (1991). A staggered finite element solution for water and gas flow in deforming porous media, Commun. Appl. Num. Meth., 7, 213–223.
58 Skempton, A. W. (1954). The pore pressure coefficients A and B, Géotechnique, 4, 143–147.
59 Terzaghi, K. von (1925). Erdbaumechanik auf bodenphysikalischer Grundlage, Deuticke, Vienna.
60 Whitaker, S. (1977). Simultaneous heat mass and momentum transfer in porous media: a theory of drying, in Advances in Heat Transfer, J. P. Hartnett and T. F. Irvine (Eds), 13, 119–203, Academic Press, New York.
61 Zaremba, S. (1903a). Le principe des mouvements relatifs et les equations de la mécanique physique. Reponse a M. Natanson. Bull. Int. Acad. Sci. Cracovie., 614–621.
62 Zaremba, S. (1903b). Sur une generalisation de la theorie classique de la viscosite. Bull. Int. Acad. Sci. Cracovie., 380–403.
63 Zienkiewicz, O. C. (1982). Field equations for porous media under dynamic loads, numerical methods in geomechanics. Proceedings of the NATO Advanced Study Institute, University of Minho, Braga, Portugal, held at Vimeiro (24 August–4 September 1981). Boston: D. Reidel.
64 Zienkiewicz, O. C. and Shiomi, T. (1984). Dynamic behaviour of saturated porous media: the generalized Biot formulation and its numerical solution, Int. J. Num. Anal. Geomech., 8, 71–96.
65 Zienkiewicz, O. C., Chang, C. T. and Bettess, P. (1980). Drained, undrained, consolidating and dynamic behaviour assumptions in soils, Géotechnique, 30, 4, 385–395.
66 Zienkiewicz, O. C., Chan, A. H. C., Pastor M., Paul, D. K. and Shiomi, T. (1990a). Static and dynamic behaviour of geomaterials – a rational approach to quantitative solutions, Part I: Fully saturated problems, Proc. Roy. Soc. London, A429, 285–309.
67 Zienkiewicz, O. C., Xie, Y. M., Schrefler, B. A., Ledesma, A. and Bicanic, N. (1990b). Static and dynamic behaviour of soils: a rational approach to quantitative solutions, Part II: Semi‐ saturated problems, Proc. Roy. Soc. London, A429, 310–323.
68 Zienkiewicz, O. C., Taylor, R. L., and Zhu, J. Z. (2005). The Finite Element Method Set (Sixth Edition) Its Basis and Fundamentals, Elsevier.
3 Finite Element Discretization and Solution of the Governing Equations
3.1 The Procedure of Discretization by the Finite Element Method
The general procedures of the Finite Element discretization of equations are described in detail in various texts. Here we shall use throughout the notation and methodology introduced by Zienkiewicz et al. (2013) which is the seventh edition of the first text for the finite element method published in 1967.
In the application to the problems governed by the equations of the previous chapter, we shall typically be solving partial differential equations which can be written as
(3.1)
where A and B are matrices of constants and L is an operator involving spatial differentials such as ∂/∂x, ∂/∂y, etc., which can be, and frequently are, nonlinear.
The dot notation implies time differentiation so that
(3.2)
In all of the above, Φ is a vector of dependent variables (say representing the displacements u and the pressure p).
The finite element solution of the problem will always proceed in the following pattern:
1 The unknown functions Φ are “discretized” or approximated by a finite set of parameters and shape function Nk which are specified in spatial dimensions. Thus we shall write(3.3) where the shape functions are specified in terms of the spatial coordinates, i.e.(3.4a) or (3.4b) are usually the values of the unknown function at some discrete spatial points called nodes which remain as variables in time.
2 Inserting the value of the approximating function into the differential equations, we obtain a residual which is not identically equal to zero but for which we can write a set of weighted residual equations in the form(3.5) which on integration will always reduce to a form(3.6) where M, C, and P are matrices or vectors corresponding in size to the full set of numerical parameters . A very suitable choice for the weighting function Wj is to take them being the same as the shape function Nj:(3.7) Indeed, this choice is optimal for accuracy in the so‐called self‐adjoint equations as shown in the basic texts and is known as the Galerkin process.
If time variation occurs, i.e. if the parameters
Usually, the parameters
Typical finite elements involving linear and quadratic interpolation are shown in Figure 3.1.
Figure 3.1 Some typical two‐dimensional elements for linear and quadratic interpolations from nodal values.
In Section 3.2,