Computational Geomechanics. Manuel Pastor

Computational Geomechanics

and of the current time t

(2.55) normal x Subscript normal i Superscript normal pi Baseline equals normal x Subscript normal i Superscript normal pi Baseline left-parenthesis normal upper X 1 Superscript normal pi Baseline comma normal upper X 2 Superscript normal pi Baseline comma normal upper X 3 Superscript normal pi Baseline comma normal t right-parenthesis

To keep this mapping continuous and bijective at all times, the determinant of the Jacobian of this transformation must not equal zero and must be strictly positive, since it is equal to the determinant of the deformation gradient tensor F^π

(2.56) normal upper F Subscript i j Superscript normal pi Baseline equals normal x Subscript normal i comma normal j Superscript normal pi Baseline equals normal upper R Subscript i k Superscript normal pi Baseline normal upper U Subscript k j Superscript normal pi Baseline equals normal upper V Subscript i k Superscript normal pi Baseline normal upper R Subscript k j Superscript normal pi

where U^π is the right stretch tensor, V^π the left stretch tensor, and the skew‐symmetric tensor R^π gives the rigid body rotation. Differentiation with respect to the appropriate coordinates of the reference or actual configuration is respectively denoted by comma or slash, i.e.

(2.57) StartFraction partial-differential f Over partial-differential normal upper X Subscript normal i Baseline EndFraction equals normal f Subscript comma normal i Baseline or StartFraction partial-differential f Over partial-differential normal x Subscript normal i Baseline EndFraction equals normal f Subscript slash normal i Baseline

Because of the non‐singularity of the Lagrangian relationship (2.55), its inverse can be written and the Eulerian or spatial description of motion follows

(2.58) upper X Subscript i Superscript pi Baseline equals upper X Subscript i Superscript pi Baseline left-parenthesis upper X 1 Superscript pi Baseline comma upper X 2 Superscript pi Baseline comma upper X 3 Superscript pi Baseline comma t right-parenthesis

The material time derivative of any differentiable function f^π(X, t) given in its spatial description and referred to a moving particle of the π phase is

(2.59) StartFraction upper D Overscript pi Endscripts f Superscript pi Baseline Over italic upper D t EndFraction equals StartFraction italic partial-differential f Superscript pi Baseline Over italic partial-differential t EndFraction plus f Subscript slash normal i Superscript pi Baseline v Subscript normal i Superscript pi

If superscript α is used for StartFraction upper D Overscript alpha Endscripts Over italic upper D t EndFraction , the time derivative is taken moving with the α phase.

2.5.2 Microscopic Balance Equations

In the hybrid mixture theories, the microscopic situation of any π phase is first described by the classical equations of continuum mechanics. At the interfaces to other constituents, the material properties and thermodynamic quantities may present step discontinuities. As throughout the book, the effects of the interfaces are here not taken into account explicitly. These are introduced, e.g. in Schrefler (2002) and Gray and Schrefler (2001, 2007).

For a thermodynamic property ψ, the conservation equation within the π phase may be written as

(2.60) StartFraction partial-differential left-parenthesis rho psi right-parenthesis Over partial-differential t EndFraction plus rho psi ModifyingAbove r With ampersand c period dotab semicolon Subscript normal i slash normal i Baseline minus i Subscript normal i slash normal i Baseline minus rho g Subscript i Baseline equals rho upper G Subscript i

where ModifyingAbove bold r With ampersand c period dotab semicolon is the local value of the velocity field of the π phase in a fixed point in space, i is the flux vector associated with Ψ, g the external supply of Ψ, and G is the net production of Ψ. The relevant thermodynamic properties Ψ are mass, momentum, energy, and entropy. The values assumed by i, g, and G are given in Table 2.2 (Hassanizadeh and Gray 1980, 1990; Schrefler 1995). The constituents are assumed to be microscopically nonpolar; hence, the angular momentum balance equation has been omitted. This equation shows, however, that the stress tensor is symmetric.

Table 2.2 Thermodynamic properties for the microscopic mass balance equations.

Sources: Adapted from Hassanizadeh and Gray (1980, 1990), and Schrefler (1995).

Quantity	ψ	i	g	G
Mass	1	0	0	0
Momentum		t _m	g	0
Energy				0
Entropy	Λ	Φ	S	φ

2.5.3 Macroscopic Balance Equations

For isothermal conditions, as here assumed, the macroscopic

Скачать книгу