Computational Geomechanics. Manuel Pastor

Читать онлайн книгу.

Computational Geomechanics - Manuel Pastor


Скачать книгу
bold u equals left-bracket u 1 comma u 2 comma u 3 right-bracket Superscript normal upper T

      Similarly, we shall use wi and vi or w and v to denote the velocities of water and air relative to the solid components. These velocities are calculated on the basis of dividing the appropriate flow by the total cross sectional area of the solid–fluid composite.

      As mentioned in Chapter 1, σij and sigma double-prime Subscript italic i j refer to the appropriate total and effective stresses, with σ and σ being the vectorial alternatives.

      Similarly, εij or ε refers to the strain components. Further, pa, pw, and p = χw pw + χa pa will stand for air and water pressure and the “effective” pressure defined in the effective stress concept in Equation (1.11) when two fluids are present.

      Sa and Sw are the relative degrees of saturation and ka and kw are the permeabilities for air and water flow.

      Other symbols will be added and defined in the text as the need arises.

      The derivation of the equations in this chapter follows a physical approach which establishes clearly the interactions involved in the manner presented by Zienkiewicz and Shiomi (1984), Zienkiewicz (1982), Zienkiewicz et al. (1990a, 1990b), etc. This is a slightly different approach from that used in the earlier presentations of Biot (1941, 1955, 1956a, 1956b, 1962) and Biot and Willis (1957) but we believe it is slightly easier to follow as it explores the physical meaning of each term.

      Later, it became fashionable to derive the equations in the form of the so‐called mixture theories (see Green and Adkin 1960; Green 1969; Bowen 1976). The equations derived were subsequently recast in varying forms. Here an important step forward was introduced by Morland (1972) who used extensively the concept of volume fractions. Derski (1978) introduced a different derivation of coupled equations and Kowalski (1979) compared various parameters occurring in Derski’s equations with those of Biot’s equations. A full discussion of the development of the theory is given in the paper by de Boer (1996).

      For completeness, we shall include such mixture derivations of the equations in Section 2.5. If correctly applied, the mixture theory establishes, of course, identical equations with appropriately chosen parameters and rheological relations.

      It seems that despite much sophistication of various sets of coupled equations, most authors limited their works to conventional, linear elastic, behavior of the solid. Indeed, de Boer and Kowalski (1983) found it necessary to develop a special plasticity theory for porous, saturated solids. In the equations of Zienkiewicz (1982) and Zienkiewicz et al. (1990a), any nonlinear behavior can be specified for the skeleton and, therefore, realistic models can be incorporated. Indeed we shall find that such models are essential if practical conclusions are to be drawn from this work.

      

      2.2.1 Equilibrium and Mass Balance Relationship (u, w, and p)

      We recall first the effective stress and constitutive relationships as defined in Equation (1.16) of Chapter 1 which we repeat below.

      or

      (2.1b)bold sigma double-prime equals bold sigma plus alpha bold m p

      This effective stress is conveniently used as it can be directly established from the total strains developed.

      However, it should be remembered here that this stress definition was derived in Chapter 1 as a corollary of using the effective stress defined as below:

      (2.2a)sigma prime Subscript italic i j Baseline equals sigma Subscript italic i j Baseline plus delta Subscript italic i j Baseline p

      or

      (2.2b)bold sigma prime equals bold sigma plus bold m p

      which is responsible for the major part of the deformation and certainly for failure.

      In soils, the difference between the two effective stresses is negligible as α ≈ 1. However, for such materials as concrete or rock, the value of α in the first expression can be as low as 0.5 but experiments on tensile strength show that the second definition of effective stress is there much more closely applicable as shown by Leliavsky (1947), Serafim (1954), etc.

      For soil mechanics problems, to which we will devote most of the examples, α = 1 will be assumed. Constitutive relationships will still, however, be written in the general form using an incremental definition

      (2.3a)normal d sigma Subscript italic i j Superscript double-prime Baseline equals upper D Subscript italic ijkl Baseline left-parenthesis normal d epsilon Subscript italic k l Baseline minus normal d epsilon Subscript italic k l Superscript 0 Baseline right-parenthesis

      (2.3b)normal d bold sigma Superscript double-prime Baseline equals bold upper D left-parenthesis normal d bold epsilon minus normal d bold epsilon Superscript 0 Baseline right-parenthesis

      The vectorial notation used here follows that corresponding to stress components given in (1.1). We thus define the strains as

      (2.4)normal d bold epsilon Superscript normal upper T Baseline equals left-parenthesis normal d epsilon Subscript x Baseline comma normal d epsilon Subscript y Baseline comma normal d epsilon Subscript z Baseline comma normal d gamma Subscript italic x y Baseline comma normal d gamma Subscript italic y z Baseline comma normal d gamma Subscript italic z x Baseline right-parenthesis Superscript normal upper T

      In the above, D is the “tangent matrix” and dε0 is the increment of the thermal or similar autogenous strain and of the grain compression bold m ModifyingAbove p With ampersand c period dotab semicolon slash 3 upper K Subscript s. The latter is generally neglected in soil problems.

      If large strains are encountered, this definition needs to be modified and we must write